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ABSTRACT. Solutions are presented for the general number distribution function in 
cascade theory and for the kth moments of the distribution. "The Janossy G-equation 
is solved. ‘The methods of solution developed are applicable to the similar problem in the 
-electron—photon case. 


§1. INTRODUCTION 


HE fluctuation problem for nucleonic showers is of the same complexity as 
| that found in electron—photon shower theory. In both cases one is eventu- 
ally faced with the task of solving one or more partial integro-differential 
equations. Janossy (1949) in an effort to solve the fluctuation problem was 
successful in setting up the fundamental equation for a special case of the general 
number distribution function. His so-called G-equation was of great value, 
for from it one could derive the partial integro-differential equations for higher 
moments of the distribution. Recently (Messel 1951 a, Janossy and Messel 1951) 
the second moments of the distribution were found both for homogeneous nuclear 
matter and a finite absorber. However, little headway was made in solving for 
the distribution function itself or for its moments higher than the second. 

In this paper we shall give the solutions for the general number distribution 
function in homogeneous nuclear matter and a number of subsidiary functions 
arising therefrom. Furthermore, we shall present solutions for the nth moments 
as well as a formal solution of the Janossy G-equation. 

Throughout the work which follows we shall make use of one assumption 
only with respect to the cross section for nucleon-nucleon collisions, e.g. the 
cross section for the production of secondary and recoil nucleons of energy 
E, and E, by a primary nucleon of energy Ey is a homogeneous function of the 
primary and secondary energies, thus : 


w(Ey; Ey, By) dh, dk, =w(Ey/ Eo, Ey/Ey) Ey? dE, dEy.  ...... (1) 


Otherwise the function w is perfectly general. 
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§2. THE DIFFUSION EQUATION FOR THE NUMBER DISTRIBUTION 
FUNCTION AND ITS SOLUTION 

The most general distribution function in nucleon cascade theory is denoted 
by ¢nlery » -'+ €ni dey... dey. It defines the differential probability of 
finding the m nucleons with energies ¢,, e, + de,(e, = E,/Eo) in the ordere,, . . ., €,, at 
the depth x in homogeneous nuclear matter, due to a single incident primary 
nucleon of energy Ey. From the function 4, all other distribution functions of 
importance in the theory may be obtained. The function 


Ee gi-3 eg Hs byl egy wa egy M) SN leg ena) eee (2) 
perm 


expresses the probability of finding any m particles with energies «,, «,+de, at 
the depth x. Its equation-has been given previously in a number of papers. 
(cf. Messel 1951 b) 


7) 
an (Er, cannes a) = =A eo 89 


+2 APS fei Levy Gig ese )eo(egten te jeje sen ratte. (3) 
kzLJ 0 
where the primes signify that neither ¢,, nor <, are included among the m — 2 energy 
ratios «,. It may be shown (Messel 1951 b) that the dependence on the 
depth x of the F,, may be split off, thus 
Poles tenia €As 2) had BVA (Sp sep, ths Sax hy yin oe ee (4) 


where 
1 (eo) (oo) 
Moe al al Files e938) dep deg SN chee (5) 


The function f,(x) expresses the probability of finding m particles of arbitrary 
energies at the a x. Also X,f,(*)=1. From (3) and (4) it is seen that 


Ag(eiy <0) &y) = = iz, Sb. Ais (Gore es Sacleiene eee (6) 
Next we define the m-fold Mellin transform of A,, as 
Gee wees ik 7 S ay?! «wags tnd ches ee. Vile ae (7) 
Apply the n-fold transform to (6) and find 
Mise) oes: sz, zs WSs 51) MAUS} 522s See eee ee (8) 
where W545) = is ie PRET CATA 4) ed ae ee (9) 


Our initial condition is determined, since there can only be one waucleon of 
energy EK, at the depth x=0, thus A,(e)=8(1—e¢) and M,(s)=1. The solution 
of (8) is now immediate and may be given as 


heeet MAS, +» +) Sp) = (2 = Ly Ae (10) 
A.(S1) 0005 Sp) = py sa 5) XP Wa) --- eo Wis xs 2.00 
axbsécHkd aAzbyH... HY FZ 


Ee Pauenes (11) may be given in the aria form 


An(5y; ceey Ga) = Il Day WS», 5s. Sytat sey Ga))5 
v=n—1 
where X’ signifies summation over all combinations of the v+1 variables s,, 


seey Say S Sse ee ae 
taken two at a time. » Svy Syt > Snp 
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The factor 2" in (10) arises due to the symmetry of the W function. The total 
number of terms in (11) arising from the (x — 1)-fold summation is  ! (n—1) !/2"-1.. 
A few words should be added regarding the above summations. The superscripts 
beside each summation refer to the order in which it is most convenient to carry 
out the summation. The first summation simply yields W functions for all 
possible combinations of the m variables s,i=1, . . . ,n taken two at a time. 
The second summation then gives W functions for all possible combinations of 
the n—1 variables s,+5,,5,', . . . ,8,-2 taken two at a time, similarly the third 
summation gives W functions for all possible combinations of the m —2 variables. 
SetSq) 5", . . . ,S;/-3 taken two at a time, and so on for the remainder of the 
summations. We give in the Appendix the function M,, found using (11) for 
the cases »=2 and 4. Using (5), (10), (11) and applying an m-fold inverse Mellin 
transform we find for the solution of the diffusion eqn. (3), 


POA Ci eet ae) = (oi oat) ate) Level Site ok Salat eae aes (12) 
where 
, 1 Uy +i00 Un tio ee onthe F “ 
= eS —(S1 (Sn " 
iD. (ni) Mae ae an €y ese Sail nee Once (13) 


is the inverse Mellin transform operator, operating on the integrand A, in (12). 

From (12) it is a straightforward procedure to obtain the formal expressions 
for other distribution functions of interest in the theory. For instance, the 
probability O(c,7,,,;«) of finding m, nucleons with energies greater than 
«E, and n, particles with energies less than <E) is obtained from (12) by integrating 
over n, of the variables <,z=1, . . . ,m from « to 1 and over the remaining 
n, variables from 0 toe. Thus 


@(€, 2, Ng; x)= 


(2-2/4! my! (n—1)!)fq(x) er i 


el 
JOY ¢€ 


L 
bs | (PNG Na hee TE 
So far we have been concerned with homogeneous nuclear matter only. However, 
on the assumption that nuclei are spherical in shape, we may apply our results 
directly to nuclei by simply replacing the factor f,(x) in the above expressions 
by R,(D4) where 
RCD Ne | Det (Dxveaate we wad (kine (15) 
0 


and D, is the average number of collisions suffered by a nucleon in making a 
diametrical passage through a nucleus of atomic weight A. 


§3. THE DIFFUSION EQUATION AND SOLUTION FOR Cy (&,.- ++; €n3 ) 


Let C,(e1, . - - ,€n3*)de, . . . de, be the differential probability of finding 
one particle in each of the energy intervals ¢,,«,+de, and any number of particles 
with arbitrary energies at a depth x in homogeneous nuclear matter, due to a 
primary nucleon of unit energy. The m-fold Mellin transform of C,, is given by 


roo 


Va(Siy «229 5u5 8)= | a Cpe CC (eqn tye, 5 A) Gey ts Gey, | see (LO) 


70 
31-2 
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The relation between F’,, and C,, is expressed by 


Cee Sek) 
Bey =| See (As ON MN rk hy) SPAM sear! i a Alb te ay (17) 
0 . 


Since the function F’,, has been determined in (12), the solution for C,, is given 
by (17). However, owing to the complicated nature of the result for F,, it is 
difficult to obtain directly from (17) explicit values for C,, for different values 
of n. 

Explicit results for V,, are required for finding the nth moments of the number 
distribution function as well as in all problems dealing with the development of 
the cascade in a finite absorber. In previous work (Janossy and Messe! 1951) 
we developed:a method for obtaining V,, without knowing the solution of the 
diffusion equation for F,. ‘The method used was far from straightforward, and 
it would have been difficult to obtain explicit results for values of m equal to or 
greater than 3. We now present a new method for obtaining V,,. The diffusion 
equation for C,, is given and its solution obtained in the form of a simple recurrence 
relation. One may proceed in either of two ways to obtain the diffusion equation. 
The first method is to use (17) and the diffusion equation for /,,; the second is to 
set up the equation directly.* 

Consider the 7 particles in the energy intervals ¢,, ¢;+de, 7=1, ..., m, at 
the depth x. At some previous depth x—6 the last change in number must have 
taken place in such a manner that we obtain the m particles in the energy 
intervals de,. ‘There exist only two possibilities: 

(a) At the depth x«—0, n—1 of the particles may have attained their final 
energies, then one of the particles undergoes a collision giving rise to the 
nth particle. Once this collision has occurred no further change may take place. 
The probability for this process is given by 
Bem" dh | Clay's oy eh ay €3 BO) We fee de 
i=1 J0 
with 


We) = | {w(e,€,) + w(e1,€)} de, = | Z(c, €;) de, =2 | mle te) detrnesnet (18) 
0 0 0 
because of the symmetry of the cross section w(¢, €9). 

(6) Similarly for the second case where we have n—2 particles in their final 
energy intervals at the depth x—@. One of the particles undergoes a collision 
giving rise to the (w—1)th and nth particles. The probability for this process is 


> CO 
pent’ dé | Caley’) 1) Gnomes X— O)wle;/e, €/e)e* de. 
0 


JS#l 


‘Thus we find 


Cle ney Ens x) ies | " = Cer. eestets Epes ey x —0)W(e,/e)e} de 
=F 1, : 1 Cn-aler > wr Eien) e; x—O)w(ez/e, €,/e)e? acl 1 eee Re (19) 
4 


We believe Janossy gave in an unpublished manuscript the same equation as we give. He was 
not successful, however, in obtaining its solution. 
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This may be transformed by differentiation into the form 


aq Cale sy Om C= nC (e), ss. gt) 
+ Br Oye hs on gs ex) teele | de 
0 l=1 : 
2 i DCR (e] ) ores Enngs Co eel eye, erelemtiaes ©. sues: (20) 
0 J#Al 


‘To solve (20), take the n-fold Mellin transform defined by (16) and find 


ra) n 
ag V n(Si> cree Sy ahh) a a 2 Se natctn, Sate) eae $y) Sane Ae guna) 
ic l= 
APR WSs, 81) Hea Sig te ee Sate SFE STs Bee eee (21) 
J#l 
where W*(s) = | etW(e\ de Gua ie ceeds (22) 
0 


If we write «(s) = 1—W*(s) then the solution of (21) becomes 


115 CARR 8 ER WAL. 01 {-» x als) | exp {6 x als) x Ws, s;) 
i=l 0 i=1 J#l 
NON ONS, arg ho ose oy a de. aan (23) 


From our initial condition we obtain for the case n=1, V,(s; «) =exp {—«(s)x}. 
Explicit results for higher values of n may now be obtained by using (23). The 
integration though lengthy for high values of m is elementary. We give in the 
Appendix the expressions for V, and V;. The general expression for V,, can 
also be given but because of its length we omit it here. t 

Applying an inverse n-fold Mellin transform to (23) we find for our solution 


Or Gs 
Pe wrete ee Chat TeV Sts ueg Sat) tee Meee (24) 
which is necessarily equivalent to the solution given by (17). ‘The above result 


D4 
may be extended to nuclei by applying the averaging operator | 2% dx). P to it: 
0 


§4. THE G-EQUATION, ITS SOLUTION AND ALTERNATIVE 
EQUATIONS FOR Fnr(&,.--,€n3 x) AND Cnr(é,,..-, €n3 *) 


Janossy (1949) set up the following diffusion equation for O(e, m, m2; x) 
defined by us in §2, 


“£ oOo ro 
D(e, Ny, N25 x) = | e~ *—9) dg s | | m1) (: : ny, No! ; a) 
0 0 1 


0 Ny! +N” =Ny 
Nol +Ng” =No 


x® (: My", Mo” 3 a) Dieta ides dino we (25) 
€9 


+ The explicit expression for C,, and Vy is given in full by Messel and Potts (1952) : see their 
eqn. (30), where the recurrence relation for Vy is completely solved. 
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In a similar manner we have set up alternative diffusion equations for F,, and 
C,, which are practically identical in form with that of (25). 


Xv 
FUE; Be, Bas 2)= | en @—9) dp fs ie Fels Bes vay ee) 
h 0 aes 
x F,,(Eql, Ey", + ++y Ey," 5 9)t0( Eq! [Epp Eq”|Eo)Ey® dE! GEy", ....+- (26) 


where the summation is carried out over all possible ways of obtaining 1, +m, =n. 
Note that whenever either m, or m.=0 then the product F,, F,,, necessarily 
vanishes. The primes on the &, indicate that none of the EZ,’ are included among 
the £,". 


Cy(Eos By 2. Enxe)= iene “ ails afoe Cy (Ege TE tessa) 
X Cy,(Eo" Ey", «+» Eng; 0)20(Ep'|Ey, Ey" |Ep)Eg-? dEy’ dEy!, «0. (27) 


where the summation is carried out over all possible ways of obtaining n, +”, =n. 
In this instance when either 7, or 2, =0 the product C,,,C,,, does not vanish. 

It is a peculiar circumstance that both eqns. (26) and (27) may be solved 
directly by means of a n-fold Mellin transform and solutions obtained which may 
be shown to be equivalent to those already given by us in (12), (17) and (23) 
respectively, but not eqn. (25), even though it is only a special case of (26). The 
solution of (25) for O(c, m,, 2, ; x) is, of course, given by (14), but this was obtained 
by solving for the general probability distribution function F, first. 

The difficulty of solving (25) for ®(e,,,7,;x) led Janossy to introduce 
a generating function defined by 


Ge Us ay Ole, Ne a) 0 (28) 
NV 


in which U is written for u, and u,, N for n,, m, and U*> for -u,™u,™. 
Equation (25) could then be put in the seemingly simplified form 


rs) ora 
— G(e, U; x) +Gle, U; x)= | | ‘ (< ie v) G (: ,U; 3) wW(e1, €) de, dey 
Ox 0/0 €y Le 
Cece (29) 
. yee Bias [ a, for ale 
with the initial condition Ge, u,, v2; 0) = 1 eee (30) 
(u for e>1. 
The solution of the G-equation is of no interest in itself; however, for the 
sake of Pee es we give its solution below: 


1 
Ge, U; x) = E—— f,(x) Ue hs lire ily Pa Mids soci na eee 


N oa Tos 
obtained from (10) and (28). 


§5. THE RTH MOMENTS OF THE DISTRIBUTION FUNCTION 


The ‘kth moment of the distribution function giving the probability of 
finding n particles with energies greater than «HZ and an arbitrary number of 
particles with energies less than «Ky at a depth x is defined by 


o"G ; ; (oe) (ee) 
ee =X & n(n "eet er n, —k+1)O(e,m, Ne ; #) = Se, x). 
u u=1 m=0 n,=0 , 
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A formal solution for S*(c,«) follows immediately from the above since 
(ce, 2,,,; x) is given by (14), hence . 


fee) (oe) & 8 rl 1 ; 
S¥ex)= ES fu(2)/ng! (mF)! | mal i elk TE MME es \der rae 
n N2=0 Ons One € 


It is obvious that to obtain explicit results for a given k, from (33), is by no means 
an easy task. In order to circumvent this difficulty we give below the diffusion 
equation obeyed by the Ath moments and methods of solution which are more 
convenient for obtaining explicit results. 


From (29), (30) and (32) we obtain 


0 = e cl att es k € a = 
5, S62) + S462) | SE 2-2 (=) 92(£ 2) S*#(£ , x) we, 5) de, der, 
Ox 0 1 ey 


0 1=0 €5 


where k’ =k/2 when k is even and (k—1)/2 when k is odd. The exponent of the 
2 gives the value 4 when 2/=& and 1 at all other times. The solution of (34) for 
k=1 may be obtained by means of a single Mellin transform, thus: 


Ale; 2s = (peas €9exp4 ins) }ids/ Sen) ed lee (35) 


In a similar manner, using the solution (35) obtained for k=1, the solution for 
k=2 may be obtained. The procedure may then be repeated for k=3, ete. 
There is no need, however, to use this lengthy method, since the solution of 
(34) may be obtained directly. Assume that the solution is given by 


1 81’ +i0 8, +400 ds ds 
aS *(e, x) = ol aril ey aon os) Zo (5. 5655 $45 %)—... 
(272) "J 57 tc 3, —tw Si 


Substituting (36) into (34) yields the differential equation satisfied by the 2, : 
which can be shown to be identical with that satisfied by the n-fold Mellin 
transform of (27) when a given term is not repeated for all combinations of the 
k variables, i.e. for example, instead of 

{W(5, 52) W(sy + $2, 83) + W(s1, 83) W(sy + 53, 82) + W(S2, 83) W(S2 + 835 81)} 
we write 3W(s,, 5.) W(s, +55, 53). The solution for S*(e, x) is hence given by 
1 (aie (ee ds, ds, 


— Biers) (ee, Se os 
(2riy* Bee ce rae 


S*(e, x)= 


8-10 


§6. DISCUSSION 

The most striking feature of the solution of the fluctuation problem is this: 
any of the solutions for distribution functions or moments dealing with z particles 
invariably lead to an n-fold inverse transform. Once the form of the solution is 
known it is a relatively simple matter to find the exact solution. Furthermore, 
it is easier to solve for the most general distribution function itself, rather than 
for special cases of it. ‘The methods of solution given in this paper are directly 
applicable to the electron—photon cascade. We shall, in the near future, present 
the solutions for the distribution functions in this case as well. 
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We have not attempted to discuss the numerical evaluation of any of the 
expressions given in our work. This we hope to do shortly. For small numbers. 
of particles the evaluation appears straightforward. 
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> APPENDIX 
We give below the expressions for M,(s,, 5.) and M,(s,, 53, 53, $4). 


M,(5, Sy) a 2W(s,, Sy). 


93 
MiA(sy, $9, $3, $4) = 31 {W(s1, Sg) W(s1+Se, 83) W (sy +S2+Ss, 54) + W(5y, 52) W (sy Se, 84) W (54+ Sat Say Sa) 


+ W(s1, 83) W(s1 +53, $2) W(sy+S2+Ss, $4) + W(5, 83)W (81 +Ss, 54) W (81+ $34 Sa) Sa) 

+ W (Sy, 84)W(Sy+54, $2)W (51+ 52+ S4, $3) + W(53, 54) W (51 +Sq, 53) W (51+ 53+ Sa, Se) 
+ W (52, 53) W(s2+Ss, $1) W(s1+-S2+ Ss, $4) + W (Se, 83) W (S2+Ss, 54) W (S2+S3+ Sa; 51) 

+ W (Se, 84) W (So+ Sa, $1)W(s1+S2+ Sa, 53) + W (Se, Sa) W (S2+Sa, 53) W (So+S3+ Sa, 51) 

+ W (53, 84) W(s3+ Sa; $1) W(s1+534+Sa, Se) + W(Ss, 54) W(S3+-Sa; S2) W (So+S3+ Sa, 81) 
+2W (Ss, S2)W (Ss, 84)W (S152, 83+ $4) +2W (Sy, 53) W (59, 84) W (51+ Sg, S254) 

+2W (sy, S4)W (Sa, 83)W (Si +54, So+Ss)}. 


The expressions for V2(s,, 5.; «) and V3(s,, 52, 3; x) may be obtained either 
from (23) or from (27) after an n-fold Mellin transform has been taken. It 
will be found that in practice it is easier to obtain explicit results from the 
latter case. 

—2W(s4, 59) 
au($1 + Sg) — ($1) — (Sp) 
4W (s1, 82) W (S14 Sa, 83) nh opt a eh 
a($y) + (Sg) — oe($1 +S) (Sy + Sy+- $3) —a(sy +52) — ag) 

sat) Se eer er er 
($y + S2+S3) —a($1) —aX(S2) — a(S) 
4W (51, 83)W(s,+5s, $2) See et eset ater sea 
a(S) + (83) — ($1 + Sg) (Sy + $9-+S3) — a(S, + $3) — a(S9) 


_ £sP [ —{a(s1) + ox(se) + a(s3)}x] — exp[ —a(s;+S_+ “al 


V4(Sy, S25 x)= {exp[ —{a(s,) ++ ou( $2) }x°] — exp[—a(s;+s2)x]} 


V3(S1, Se, $33 %)= 


a(S, -+ $2+ $3) —ax(sy) — a(S) — a(S) 
4W (so, 53) W (s2+-53; 51) {Sek ee —exp[—a(s,+5_+5g)x] 
(Sg) -+-o(s3) — a(Sa+ $3) a(Sy+-S2+$3)—a(S2+ 53) —a(sy) 
ae exp[—{a(s1)-+ a(s2)-+ a(s3)}x] —exp[—a(sy+ 5.4 “vel 
(Sy + $9+53) —x(S1) — (Sg) —a(sg) : 
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ABSTRACT. ‘The complete analytical solution of the fluctuation problem for nucleon 
cascades developing in a finite absorber is given. Expressions are obtained for the general 
number distribution function and for the nth moments. 


§1. INTRODUCTION 

N a previous paper (Messel 1952, hereafter referred to as I) an analytical 
if solution of the fluctuation problem for nucleon cascades developing in 

homogeneous nuclear matter was presented. The methods of solution 
developed in I were of general applicability and in the present work we apply 
them to solve the similar fluctuation problem in a finite absorber. 

In §§2 and 3 of this paper we give the diffusion equations and their analytic 
solution for the number distribution function which expresses the probability 
of finding a stated number of particles with given energies at a given depth. 
From this distribution function one can in principle obtain all other number 
distribution functions of interest in the theory. An example is then discussed in 
detail for the case of four particles. In §$5—7 the diffusion equations for the 
mth moments and their solution are given. ‘The actual evaluation of the 
expressions obtained is not attempted in the present paper but is left over to 
subsequent work, where the problem of evaluating the function appearing in 
the solution of the fluctuation problem both for nucleon and electron—photon 
cascade theory will be discussed. 

Throughout our work we shall retain the same assumptions regarding the: 
cross section for nucleon—nucleon collisions as in I. 


PART I—THE DISTRIBUTION FUNCTION 
§2. DIFFUSION EQUATIONS 

Let H,(Ey; £y,..., £,; 9) dH, ... dE, be the differential probability that 
after a depth @ (measured in units of the interaction mean free path) a primary 
nucleon of energy F, has given rise to m nucleons with energies in the ranges 
E, E,+dE,,j=1,...min any order; and let Q,(E); E,,..., H,)dE, ...dE,d@ 
be the probability that a nucleon of energy Ey collides with a nucleus in travelling 
a distance dé and gives rise to n nucleons with energies in the ranges E,, E; + dE,,. 

j=l1,..., min any order. Assuming nuclei to be spherical in shape 


D 
Oe(iigy Boh 22H \de { © F( Bess Hig das. Bh is'ee)20 dx[D 42 sina: (1) 
0 
where, as in I, F,(E); E,,..., £,; «)dE,...dE, is defined as the differential 


n? 
probability that after a distance x in homogeneous nuclear matter a primary 
nucleon of energy Ey has given rise to m nucleons with energies in the ranges 
E,, E;+dE;,j=1,...,nin any order. D4 is the diameter of a nucleus of atomic 


number A. The function F,, has been determined in I and shown to be a 


‘474 H. Messel and R. B. Potts 


homogeneous function of the primary and secondary energies. Hence it folloms 
that O,, and therefore H,, are homogeneous functions as well. 

From a consideration of all possible first nucleon—nucleus collisions follows 
the diffusion aes for H,, in the form 


Oy 1 i) (oe) } 
H,,+H,= = =| dE i a dE wQ(Eo; Exo sey Eo) 
00 lk ee) 0 
» 
pad ia | Fg (Ejo; Ex, sey Ejng)s 0); eee (2) 


nv) j=1 
‘where &,,») signifies summation over all compositions of £,,...,H,, into 


v groups Ey,,...; yaa} ---3 Egy - - +» Eyn@ such that Xf_,n(j) = in all possible 7| 


ways, n(j)>0. The order within each group of n(j) terms is immaterial. There 
are just 2!/n(1)!...n(v)! ways of choosing the v groups n(j), 7=1,..., v from 
these, yokes. 

An alternative form of the diffusion equation may be obtained by considering 
all possible last collisions; in this case 


56 Ey ns ly He lee Eee ee 


rs) fe) 
H,,+nH,,= = DD OG; 
=1 (k)/0 
where 2X, signifies summation over all compositions of the £,,..., Z, into 


two groups £y’,..., Ey’; Eya1,..., ,'. It is this form of the diffusion equation 
which will be used in finding the solution for H.,. 


§3. SOLUTION OF THE DIFFUSION EQUATION 


The partial integro-differential equations (3) can be transformed to a recur- 
rence relation by taking an -fold Mellin transform with respect to the energy 


variables and a single Laplace transform with respect to the depth dependence. 
Let 


dS CE, RRR ays) | 
ss | Be | dE, | d0(E,|E,)® . . . (E,/Eq)’« €-?°H, (Ey; Ex, «> En} 9) 
0 ~ 0 ¥ 0 


BN ee ee (4) 
N,(s1, eS) Sn) = IF dE, L208 |, dE,,(E,/Ey)* aN (£,,/Eo)"O,(Ey; EF ak OS 5) E,); 
then eqn. (3) may \be-transformed into 2.4, 0) ws ee ee (3) 
(p+n)R,= = mt NGS ys => Se) Rae ball Sey = ee ee Ee Oa eee 
Sep ae (6) 
From the initial condition H,(E); FE; 0=0)=8(E)—E) and (3), (4) and (5) it 
may be shown that 
RS; Pa VPA). | eee ee ee (7) 
where h=1—2{1-(14+Dy)exp(-—D,)VWDe. 22... (8) 


The term on the right-hand side of (6) containing R,, for the case k=1 may be 
taken over to the left-hand side yielding 


(p+nh)R,,= 2B Nilsy', +++ 4 Se) RicnsilSprn «009 Sp 9 St eee. «ee | Deeds 
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The solution of the recurrence relation expressed by (9) may be given in the 
following compact form: 


1 $ : 
Ray = prh 2 2 NV ni+(Saj—van 62+) Sagy Saar + +++ + Sneath +[qG)+ LAE, 


h e j=m com 


eas (10) 
where 
(a) =, signifies the sum over the 2”! compositions of n; c is the composition 
n(1), n(2),..., m(m). Xi yn(j)=n, m is the number of nucleon—nucleus 


collisions, n(j)+1 gives the number of nucleons arising from the jth nucleon— 
nucleus collision. 


(6) g(j)=n(1)+ ...+n(7), q(0)=0 and g(m) =n. 
(c) =X signifies summation over all combinations C%#}t} of the g(j)+1 


symbols 5), 59, .. . , Sq. Soggak «+ - +5n41 taken n(y)+1 at a time. 

It will be found most expedient to carry out the summation », first, then 
the product [I]}_,, and finally to sum from the right. 

hie 2" * een of n describe the different ways in which a nucleon 
can, by collision with nuclei, give rise to the +1 nucleons. If from the 
jth nucleon—nucleus collision n(j)+1 nucleons emerge and there are m collisions, 
then &”,n(j)=n; hence the composition n(1),n(2), ...,m(m) describes a 
eee in ae the primary nucleon collides with a nucleus giving rise to 
n(1)+1 nucleons; one of these nucleons then collides with a nucleus giving rise 
to n(2)+1 nucleons, etc., until the mth collision in which n(m)+1 nucleons 
emerge. ‘The total number of nucleons is then 2+1. 

It is evident from (10) that the p-dependent factor of R,,,, is 


TT” 91 /{p +[9(j) + 11h} 


and that this factor is the same for each composition c. The case Ry, 1.e. n=3, 
will be discussed in detail in the following section. 
The solution for H,, is now immediate: 


Dib Oo ORR od Se pel TE oN DPN Sal Ce a Rex st Praha) Hanes ME eg (11) 


where I,, is the inverse Laplace—Mellin transform operator : 
nee Uy, +400 A+io E 8i+1 E $4 
I= oan Br" { cea asap |, eeu (3) (2) en, 
On Ut u,—t0 A—iw Ey E,, 


From the distribution function H,, whose solution is given by (11) all other 
number distribution functions of interest in the theory may easily be obtained. 
For instance the probability of finding m, particles with energies greater than F 
and nz, particles with energies less than E at a depth @ due to a primary of energy Ey 
may be obtained by integrating over m, of the variables of H,, from E to o and 
the remaining 7, from 0 to £; hence 


1 (oa) o rH EE 
(Ey; E, m, ny; 6) = ——— | i) | el 2 CT pei) aan (13) 
Ny: Nog:JE E/0 0 
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The physically interesting quantity is, however, the probability of finding 


n, particles with energies greater than E and an arbitrary number of particles — 


with energies less than . at a depth 0, which is given by 
. E 
O(Ey, E,n,;6)= > sal, ae es all fide ode = hae (14) 
0 


Ne =0 1, | 
The integrations in ies and (14) over the energy variables are easily 
performed since the function H,, is a homogeneous function of the primary 
and secondary energies as seen from the solution given in (11) and (12). 


§4. THE CASE R, (S1, Sa, Ss, Sa; Pp) 


It is evident that the solution for H,, given by (10), (11) and (12) is by no 1 


means elementary. In view of this we have thought it advisable to discuss one 


case in full, the case R,. 
First it is necessary to list the compositions of »=3; there are 2?=4 in all: 


Composition m ny)+1 Form 
c(i) 15151 3 ORD) as fig. 1 (a), (6) 
2) ey 2 26 as fig. 2. 
cS) all 2 352 as fig. 3. 
c(4) =3 1 4 as fig. 4. 


The values of m in the second column give the number of nucleon—nucleus 
collisions, the values of n(j)+ 1 give the numbers of nucleons emerging from the 
m nucleon—nucleus collisions, and the diagrams represent possible ways in which 
the traverse can occur. These will be discussed further below. 
For the composition c(1) the corresponding term in expression (10) is 
1 
wn a x 
(PEIN P+ DAH SAD + AM) og 860 80 fe Nala b 80 +89) 
si No($15 $2 + S3 +54) 
C2 


- SEHOTGTINGTD = x ey 54){No(S2, $3 + $4) No(S1, 52 +53 +54) 

i N2(51, 83 +54) No(Sey 5) + 53 + 54) 

+ No(S2, 51) No(S3 +54, So +5;)}, «2.00. (15) 
where U4 signifies summation over all combinations of s,, 53, 83, S, taken two 
at a time. When this last summation is carried out the result becomes entirely 


analogous to the expression for M,(s,, s., Ss, $4) given for homogeneous nuclear 


matter in the appendix of I. This is to be expected since the composition 
c(1)=1, 1, 1 corresponds to a traverse in which each nucleon—nucleus collision 
gives rise to two nucleons just as in nucleon—nucleon collisions. 

It is interesting to notice that the terms in (15) after the last summation 
has been performed can be represented by ‘trees’ of the type given in fig. 1. 
The terms are constructed from the trees by reading from bottom to top and 
adding the s variables at each collision. ‘The two trees given in fig. 1 correspond 
to the terms 


N2(53, 54) No(52, 53 + 54) No(S1, S2+53 +54) and No(s3, 54) No(s, 52) No(S3 + Sq, 51 +59). 


The remaining terms in (15) may be obtained from fig. 1 by means of a simple 
permutation of the variables. 
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The p-dependent factor in (15) may also be read off from the tree, 
fig. 1(a). Reading from top to bottom, there are at first one, two, then three 
and finally four branches at each level; the corresponding p-dependent factor is 
1/(p+h)(p+2h)(p+3h)(p+4h). These factors are always determined by trees 
of the type in fig. 1 (a) (and not 1()). It is known that the factors depend only 
on the composition, and there is a tree of this type for every composition. 


5, +52 
Sgt Sa 
$i S233 54 

(B) 
Fig. 1. 
5p* 53+ Sa 
oa 
S, Se sew as GT he US 8S, ies Pe eke, ? 

Fig. 2. Bivees. Fig. 4. 


For the composition c(2) the corresponding term in (10) is 
INET ING TH 2 N3($25 S35 $4) a N4(51, 82 +53 +54) 
Sj oleae eee 
— (p+h)(p+2h)(p+4hy ° 

+ N,(S3, S4, 51) No(S2, $3 +54 +51) + No(S45 $15 $2)No(S3, $4 + 51 + 59) 
NAS; Sota) ah Sy, Sytoee Te Sa)t- Pet) Ne  yis. (16) 

‘The tree from which these terms may be obtained is given in fig. 2. 
For the composition c(3) the term is 


N3(Sp, $3, S4)No(S1, 52 +53 +54) 


oa WEDET EAE ee y i 
(D+Ap + 3A(—+4N) os Na(Say $4) 2% Na(S1y S25 83+ 84) 


1 ; 
(p+ Ip + Mp + ah) 
+ N9(Sq) $1). N3(S2, 83) $4 +51) + No(S1, S2)Noa(Sa5 Say 51 + 52) 
+ No(s, 53)No(Sa, $1, $2 +53) + No(S1, 53)No(S2) Say 51 +55) 
PUVa(Sas Sa)ValSus Sar Setsa)f, ett tng (17) 
corresponding to the tree given in fig. 3. 
Finally the composition c(4) yields the term 

(1/(D+A)(P+4h)}-NalS1s S25 $3 Sa)y nes (18) 

and the tree representing it is given in fig. 4. 


The solution for R, is now obtained by adding the terms given by (15), (16), 
(17) and (18). 


3 


Sg $4)Na(S1, Se 53 +54) 


——" 
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PART II—THE MOMENTS 
§5. DIFFUSION EQUATION FOR THE MOMENTS 

The nth factorial moment of the distribution function ® found in (14) is 

defined to be er 
T (£9; £; 6)= beany (EN ERE OVA 4th sgt aes (19) 

Since ® is known, we have in principle the solution for the T,. ‘The solution 
for T,, as expressed by (19) is not, however, very susceptible to explicit evaluation. 
A more suitable expression for the moments is found by solving their diffusion 
equation, which will be given below. 

The starting point is the G-equation corresponding to the function ® (see 
Messel 1951): 


fa} foe) 1 co [o-e) 7. 
OtG= S| dy... | dlyoQ,(Ey; ry -++> Ero) 1 G(Eyy; Bu; 9), 
(Pett UE) 0 j=1 


where GUEy EF; 050) = Os EO Ro). ee eee (21) 
k=0 


Since we have an analytical solution for the ® expressed by (14), the solution 
for the G-equation is immediately obtained from (14) and (21). It follows 
from (19) and (21) that 

Lone OD) =(@2G cn) ee eee (22) 


By successive differentiation of (20) and making use of the relation (22) and the 
fact that G(E); FE, uw=1; 6)=1, we find the diffusion equation satisfied by T,,: 
0 


ive} 1 co ioe) 
56 Uae be Set [. dE)... I. dE» OA(Eo; Exo, - - - » Ero) 


n! a 
x 2 Ney i Tolls £30), eeece (23) 
where & signifies summation over all compositions of n into r groups n(j), 
G21, on (mays be- zero): 

The solution of eqn. (23) for the cases m=1 and 2 has previously been 
obtained (Messel 1951). One could obtain the higher moments by successive 
substitutions into eqn. (23). This proves to be a very lengthy task for values 
of n>3. Equation (23) may however be solved in a manner analogous to that 
given in I for homogeneous nuclear matter. In order to do so, we must first 
consider the solution of a subsidiary diffusion equation, given below. 


§6. DIFFUSION EQUATION FOR J, 
Let J (Lo; Be os. ea ae ees dE, be the probability of finding 
n particles with energies in the range E,, E;+dE; at a depth 0, and any number 
of particles with arbitrary energies, due to a primary of energy Ey. As for the 
case of the H,,, two separate diffusion equations may be set up. 
By considering all possible first collisions we find 


rs) : eg 1 (ce) r co 
Hintia= BS { aE yy | dEO (EO Bae Bo) 


0 


x » J i (Ejo5 Ey, eee y Ejny)3 6), aieeteie: 6 (24) 


nr) j=1 
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where &,,,) signifies summation over all compositions of E,,...,E£,, into 
r groups. Fy, ..., Fina} +++} By). ++; Epa, Such that 2,_,2(j)=2 and where: 
n(j) may be zero. There are just n!/n(1)!...n(r)! ways of dividing the 
E,, ..., E,, into the r groups. 

An alternative form of the diffusion equation is obtained by considering alk 
possible last collisions; in this case we find 


ra co 1 00 
J,+nJ,= > BB email, aa. fs. d 
06 b=1 r=k (kt) (r—)! a 1 a iat 


x 0.(n; Ey’ ee | Ey ? rap eR) End n—naa(Eo3 Es DWC) Ee 3 0), 


Soh: (25) 
where 2 ;) signifies the summation over all compositions of E,,..., E,, into 
two groups £,’,..., E,, and Ey 441... Ey > 

Let us now define Ay 
Bette Perey | CAR yy tab Egy )2e de) 8 82 2 (26) 
0 


where from I 


Che ee | 1 hela | -LSppled da UD gid Doheny eve.) 
= 0 0 ‘ 


hence using (1) 


co i foe) r 0 
SE ay ee | ee | _ Ent Qnra Eo; Ex +++» Ents) 


The function D,, expresses the probability that a primary nucleon of energy Eo, 
on colliding with a nucleus will give rise to m nucleons with energies in the 
intervals dE;,7=1,...,m, and any number of nucleons with arbitrary energies. 
In I the solution for C,, was left in the form of a recurrence relation. It may 
be shown that the explicit solution for C,, is 


CR ey ee (29) 


1 


Ay 1 
where V AY ae | Sn ) > ey a a I = 
we 1 +1 p) pta(sy+ eels +Sya9) jan comb 


WSs flat + + Sai {Pta(sy) +... Falsyp)ta(Sejaat »-» +5n41)}4 


and the notation is the same as in (10) with the composition 
n(1)=n(2) = - =n(n)= =1. We refer the readers to I for the definitions of « and W. 
Now since C,, is known, the value of D,, may be obtained by perlonaine 
the x integration as in (26). 
Using the relation (28) for D,, eqn. (25) may be reduced to the form 


0 “2 bs / / / / 
Ne fee eS [ Dy BEE i. op Bn baa 
00 k=1 () 10 
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The solution for J,, may now be obtained from (31) by applying a single 
Laplace and an n-fold Mellin transform and solving the resultant equation in a 
manner similar to that used to find H,. We find 


Dee ee a Oe eee (32) 
1 1h 
Vals > 5 Sas D) = Fh ee = Bigs 
i Sakyeay eine + ss Sac qtpau ce acoaete ne 
x {PtA(S) +... FA Sap) +MSa pit «++ FSn4 Fs veers (33) 
where B,, is the -fold Mellin transform of D,, and 
h(s)=1—2[1 — {1 + Dya(s)}exp{— Dya(s)}]/{D4a(s)}*. ...... (34) 


The terms which occur in (33) can be constructed from trees just as for (10). 
The function-B, occurs in exactly the same way as the function N, but the 
p-dependent factors are now different for each term in Xo». They can, however, 
be read off from the tree. For example, fig. 2 corresponds to the term 


B5(Sq, 83) 54)Bo(S1, 52 +53 + S4){[P + A(sy + 59 + 53+ 54)] 
x [p+ h(s1) + A(S. + 83+ $4) IP + A(S1) + Ase) + (Ss) + A(Sa) 
in the expression for Y,. The p-dependent factors are written down from the 


tree by reading from top to bottom and at each level choosing as the variables 
of h the values of s on the branches. 


$7. THE SOLU TTPONSS OR ein 
We are now in position to give the solution for T,,. Assume the solution to be 


rUy +400 PU, +4 00 —(Sit... +S) 
TBs Bs O= off (By 


(252)8 U,—t0 Uy — i Ey 
ds ds. x 
& Lal Saini es SetO) cae ee ee (35) 
Sy Sn 
It remains for us to determine Z,(s,,...,5,;9). Substitute the value of T,, 


given in (35) into eqn. (23). The resultant differential equation for Z,, can then 
be shown to be zdentical with that satisfied by the n-fold Mellin transform of 
eqn. (24) for J,, when a given term is not repeated for all combinations of the 


n-variables s;,7=1,...,n. Hence the solution for 7, is given by 
1 Us +000 Uy, +t / AY\—(si+... +8) 
T(E); E; 0—)==— in = 
“A 0 E 0) ae Lee: @ 
= ds ds 
Y De eee 
x Sis ) 5nd ) Si S, 5) (36) 
= 1 *A+¢0 
WiCTES = Viel Si) a ceee 8 ut) ani | CPE (Si, aie gist D) CPaaat mee (37) 
Tt) pj Hm 


‘This completes the solution of the fluctuation problem in nucleon cascade theory. 
There only remains the explicit evaluation of the different analytical expressions 
arising above. We hope to consider this problem in some detail in the near future. 
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On the Scattering of Fast Electrons by Nuclei 


By L. R. B. ELTON 
Wheatstone Department of Physics, King’s College, London 


Communicated by H. S. W. Massey; MS. received 7th February 1952 


ABSTRACT. It is shown that, contrary to expectations from the non-relativistic case, 
the phase-shift C,, between the partial waves of electrons scattered by a point nucleus and 
by an extended nucleus respectively is positive for sufficiently large n, and that this is a 
spin effect and not merely a relativistic one. An approximate formula for C,, is evolved 
and its usefulness is discussed. 


wis UN ERO DUGLLON 

ECENT calculations, by Elton (1950, to be referred to as 1) and Acheson 
Rs on the one hand and by Parzen (1950b) on the other, of the 
scattering of high energy electrons by a point and extended nucleus 
respectively have shown up an apparent disagreement, in that Parzen finds that 
the. phase shift ¢,, (his 6,) between the partial waves for the two cases tends 
to a as m tends to ©, while Acheson and the author assurned that it tends to zero 
monotonically. Since we both find that ¢)<0, Parzen concludes that tan ¢,,<0 
for small m, but is greater than 0 and tending to +0 for large n, while Acheson 
and the author assumed that tan ¢,,<0 for all m and tends to —0 monotonically 
for large n. It will be shown now that neither view is correct, in that ¢,>0 and 
tan €,>+0asn+o. It follows from this result that phase shifts ¢,, have to be 
calculated up to values of m for which ¢,,>0 and negligible, since a negligibly 
small negative value of ¢, may be followed by a larger positive one. ‘The 

calculations in I were therefore taken a stage further and ¢, was calculated. 

The notation of I will be followed throughout in this paper. 


SJ PRO OR TE Aa sey 0 

As Parzen has pointed out (1950b, eqn. (17)), @,>0 as no, where @,, is the 
phase difference between the partial waves corresponding to the regular and 
irregular solutions for the point nucleus, and hence it follows from I(12) that 
tanl,-+0. To show that ¢,-+0 we prove that for large m the functions G,,(7) 
and %#,(r) which are the regular solutions of the radial equation for a point 
nucleus and an extended nucleus respectively, have no zeros inside the nuclear 
radius, and that therefore on the nuclear radius they cannot be out of phase by 


as much as 7. 
Both Y, and %, satisfy the differential equation | (8) 
Ce] |) ee ee ON eis (1) 
where 
MRO Vg ee Ve etic le ye he. Le 
r Be Pa—“r 4 40% 2a.’ | 
and «=(W—V+mce?)jhe. Therefore, provided f,(r)>0, e7">0 as long as 
of >0, 1.e.0%’ increases as long ase’ >0; and since from the power series expansion 
of o¥ it is clear that for sufficiently small 7,” andc”’ have the same sign and can 
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alt) = — Re + 
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therefore be taken positive, it follows that c/ is positive and increasing as long 
as f,(r)>0. Now for large n and 740, but small, the dominant term of f,(7) 
is clearly n(n+1)/r2, and as r increases, all other terms of f,(7) remain small 
compared with n(m+1)/r? except possibly the first and third. Hence /,(7) 
vanishes first when 7 is approximately given by (1—2WV/h?c*k*)k?r? =n(n + 1), 
ie. when kr=O(n). Hence Y, and #,, do not vanish until kr =O(), which 
proves the proposition. 


§3. PROOF THAT tan &= +0 
It is next shown that for sufficiently large m, tan ¢,>0. In this proof affixes 
P and E denote quantities referring to a point nucleus and to an extended nucleus 


respectively. 
Equation (1) can be re-written for an extended nucleus 
d?of n(n +1) 
Se S- ce ete U(r) fof = V ,(r)of fo Peitarerete (2) 


where U,(r)=f,(r)+k?—n(n+1)/r? and V,(r)=U,"(r)—U,*(r). Thus 
V,(r)=0 for r>R, the nuclear radius. We choose a solution %, which 
satisfies the following conditions: 


A ,(0)=0,- # ,,~const. exp{i(kr+y log Zkr)}, y=Ze*/hv.-.. 2: (3) 
Let G,, and ¥,, be the solutions of the equation 
d?of n(n + 1) 
7 t & Pages ur) | of sQyrs ts ti ee (4) 


which satisfy the following conditions: 
G(0)=0, GY,~p? sin (kr+y log 2kr—4n7+7,"), 


JF ,—~pPR exp {i(kr+y log 2kr—4nw+7,F)t, ss eae (5) 
where »p =(W+me?)/ic. Then following Mott and Massey (1949) we have for all r 
dG, Che agen 
In pod te ==] cee eee (6) 
andhence’  #', =. | IV,H,dr-F, I Gy ge See (7) 
oa) 0 


To this solution we must add Y,,, since when V,(r)=0, %,, must reduce to G,,. 
Therefore 


A ,~e*? sin (kr +-y log 2kr — nz +7,°) 
—u?k- exp {i(kr + y log 2kr—4n74+7,P)} i CoH Sart ee (8) 
Now ts 
A ,—~e-\? sin (kr + log 2kr — dna +7,,") 
=p"? exp(—7€,,)[sin (kr +y log 2kr —4n7+7,,°) 
+ {exp (2i¢,,) — 1} exp {i(kr+y log 2kr—4n7+7,?)}], ...... (16) 


where,,=7,," —1,'- Comparing (8) and (10) we have, aftera little rearrangement, 


; ee ees 2 rr 
sin £, = ra GV Gt Ce PER (11) 
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Now for large n, ¢,, is small and hence to first order %,, can be replaced by Y, 
inthe integral. Also V,(r)=0forr>R. We thus obtain the approximate formula. 


wh 
eee al GeV EE “dete dip cba meenoens (12) 
RJ 
Therefore, for large n, the sign of ¢,, is opposite to that of V,(r). But 


1 n+1 / , 
V (0) = page Van Vn) 2W—Vy—Vy) + "2" (@ - 2) 


ap aH 


3 (on? ap”? 1 fay’ ap” 
ealcasteet) me edezal 
and since all terms are O(r-?) near the origin, it is clear that for an attractive 
potential ¢,, will have the same sign as apay’—apep’ for sufficiently large n. 
Now ap=u+A/r, where A=yv/c, and we can write «,=+Ag(r), where for 
any physically likely charge distribution g(r) decreases and rg(r) increases. 
monotonically with 7, and 0<g(r)<r-, all for r<R. Then for r<R, 


ge =0, re 420) Stee | er <1 )r?. 
Thus the force at any point is less than the force due to the corresponding point 


nucleus. (Note that the converse of this, 1.e. that if |g’| <r-®, then rg(r) increases 
monotonically with 7, is not true.) ‘Therefore 


py’ — pep’ =(u+A/r)Ag’ + (u+Ag)A/r? 
==Al (gl: Ljr*)-P ar (2 2/7) 
0, 


using the above inequalities. Thus ¢,>0 for sufficiently large z. 

This rather surprising result is peculiar to the Dirac equation, since for the 
Klein—Gordon equation, V,(r)=(Vg—Vp)(2W-—Vy—Vp)/h?e?, and for the 
Schrédinger equation, V,(r)=2m(V_—V p)/h?, and in both these cases ¢,,<0 
for all m since Vy >V> for all r. Thus the effect is not a merely relativistic one, 
but one due to the spin of the electron. 


§4. USE OF THE APPROXIMATE FORMULA 

The accuracy of (12) has been checked by comparing values of ¢,, obtained 
from it with the exact values, for electrons of energy 40mc? being scattered 
by gold nuclei. ‘The values are given in the table and it is clear that even for ¢, 
the approximate formula gives good results under these conditions (B and C 
refer to uniform and shell nucleus respectively). Lastly ¢,° was calculated and 
it was found that ¢,°= —0-00002. Thus ¢,° is still negative, but it is clear from 
the trend of ¢,, that all further phase shifts can be neglected. 


Phase shifts ¢,,%, ¢,,° obtained by exact calculations and from the 


MF n 
approximate formula (12) 
Phase shift C8 fo" £8 €,° 
Exact —0-196 —0-271 —0-002 — (0-004 
Approximate —0-231 — 0-321 —0-002 —0-004 


The approximate formula will be particularly useful in calculating the effect 
of different models, once the phase for the spherical shell model has been 
calculated exactly. This is easiest to calculate since %,,° is a Bessel function. 


32-2 
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Further, since in that case %,,° is known, the phase shifts for other models can 
be obtained from the formulae 


R 
Lone F | G(UB—Up)G, dr snes (13) # 
0 
and 
; bL jee A ~ ¢ 
ee | HO By Lr a. ate (14) 
0 


where the affix E refers to any model for an extended nucleus, other than the 
shell model, and the correct result must clearly lie between the results obtained 
from (13) and (14). This method is likely to give more accurate results than the 
variation principle evolved by Parzen (1950a) for the same purpose, since 
Parzen in effect replaces Y,(kr) or %,,°(kr) by the less accurate Bessel function 
J,,(kr) and neglects squares of the potential. A numerical check has been carried 
out, by calculating £,8 —¢,° by Parzen’s method and from formula (13), which 
confirms the superior accuracy of the present method. 
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ABSTRACT. 'The gamma-rays emitted in the decay of *Re have been investigated 
by means of absorption and coincidence techniques and by the use of a thin lens beta-ray 
spectrometer. ‘The gamma-ray intensities (expressed as a percentage of the total number 
of beta-disintegrations) have been measured and a probable decay scheme is proposed. 


§1. INTRODUCTION 

LOW neutron irradiation of rhenium produces the two radioactive isotopes 

S *8°Re and !**Re which decay by beta-emission to the osmium isotopes 186Os 

and '**Os with half-lives of 90 hours and 18 hours respectively. 18%Re has 

previously been investigated by Miller and Curtiss (1946), Goodman and Pool 

(1947), Cork, ShrefHer and Fowler (1948) and Beach, Peacock and Wilkinson 

(1949), but the only gamma-ray intensity measurements have been those of 

Miller and Curtiss who gave rough estimates of the relative intensities of the 
gamma-rays. 


* Now at the Department of Natural Philosophy, Glasgow University. 
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§2. THE DETERMINATION OF GAMMA-RAY INTENSITIES 

The secondary electron spectrum ejected by the gamma-rays from a 
thick lead radiator (150mgcm-?) was measured in the spectrometer and is 
shown in fig. 1. Photoelectron peaks corresponding to gamma-rays of 
energies 0-152+0-001 Mev (K and L levels), 0-476+0-005 Mev (K only), 
0-638 + 0-005 kev (K), and 0-933 + 0-008 Mev (K) were observed (the 0-152 Mev 
peak is not shown in the figure). These values are in good agreement with those 
obtained by Miller and Curtiss and by Beach et al. The 1:40 Mev gamma-ray 
reported by these authors was too weak to be detected in the present experiment. 


Ey = 0-476 Mev 


E’y = 0-638 Mev 


Counts (min-!) 


Ly = 0-933 Mev 


1000 - ta! 
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Fig. 1. Photoelectron spectrum of 1**Re y-rays. 


Using the method of Ellis and Aston (1930) as modified by Richardson (1950) 
the relative intensities of the 0-476, 0-638 and 0-933 Mev gamma-rays were 
estimated from the secondary electron spectrum to be in the ratio 1:2:1-5 with 
an accuracy within about 10%. Since the loss of electrons from the peak of 
the line by multiple scattering and anisotropic emission of the photoelectrons 
may be very serious at low energies, no attempt was made to compare the 
intensities of the hard gamma-rays with that of the 0-152 Mev radiation. 

The absolute values of the gamma-ray intensities were determined by taking 
an absorption curve of the gamma-radiation in lead, using a Geiger counter 
whose efficiency as a function of gamma-ray energy had previously been 
estimated. ‘The curve could be resolved into four components corresponding 
to gamma-rays of energies 0-15, 0-42, 0:96 and 1-32Mev. ‘These values are in 
satisfactory agreement with those obtained on the spectrometer. No component 
corresponding to the 0-638mMev gamma-ray could be resolved and it was. 
assumed that the ‘0-42Mev’ component contains the 0:-476Mev and the 
0-638 Mev gamma-rays since these have similar mass absorption coefficients. 
The relative intensities of the four resolved components were determined from 
the absorption curve, and the absolute intensities, expressed as a percentage of 
the total number of beta-disintegrations, were then estimated using a calibrated 
beta-counter in conjunction with the gamma-counter. The results are: 


’ Gamma-ray energy (Mev) 0-15 0-42 0-96 yey 
Gamma-ray intensity 70+ 10 Tae 2 5a 2 sae 72 
(Total B-ray intensity = 100) 
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These figures are of low accuracy because the continued subtraction of the 
higher energy components from the absorption curve leads to cumulative errors. 
They are combined with the more accurate relative intensity values obtained 
by means of the beta-ray spectrometer to give the following estimate of the 
absolute intensities of all the gamma-transitions: 


Gamma-ray energy (Mev) 0-152 0-476 0-638 0-933 1-40 
Gamma-ray intensity (0c LOW Sete Oat cia, 4:54+2 
(Total B-ray intensity = 100) 


§3. THE CONTINUOUS BETA-RAY SPECTRUM AND THE 
CONVERSION ELECTRON SPECTRUM 
A Kurie plot of the continuous beta-ray spectrum is shown in fig. 2. In the’ 
region above W=3m,c? it is substantially linear and leads to an end point energy 
of 2-07+0-02mev. However, because of the limited resolution of the 
spectrometer and the general high background due to the large gamma-ray 


"FE (plibin, p20, BOE 12 bs OTR eh da eae 
W (Total Energy Including Rest Energy) (v2,c? units) 


Fig. 2, Kurie plot of 18%Re continuous f-ray spectrum. 


intensity it is not certain whether this end-point energy refers to the transition 
to the ground state, the transition to the 0-152 Mev level in 188Os or to some 
mean of these two transitions. This difficulty does not affect the decay scheme 
put forward in §5, but does give an uncertainty to the exact disintegration energy 
(shown as the dotted line at ~2-1 ev in fig. 5). 

Below W=3m) ec? the beta-spectrum of 188Re was obtained by subtracting 
the contribution of the 18*Re spectrum (which was measured after the 188Re had 
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been allowed to decay for a period of 10 half-lives), and is therefore subject to 
cumulative errors. ‘The Kurie plot in this region is still substantially linear, 
but this cannot be taken to indicate that no beta-ray branching occurs since the 
plot is not of sufficient accuracy to show a branching of less than about 10%. 
In fact, gamma-ray intensity measurements and B-y coincidence experiments 
(see §4 below) show that such branching is present. 

The conversion electron spectrum was investigated in an attempt to determine 
the multipolarity of the 0-152 Mev gamma-ray. The spectrometer resolution 
(An/y at half-maximum of a conversion line, where 7=electron momentum) 
was 5%. ‘The Geiger counter used for detecting the electrons had a thin nylon 
window which transmitted electrons with energies greater than 10kev. 
Rhenium metal was precipitated from a solution of potassium perrhenate 
(KReO,) and deposited as a source of thickness 0:5 mgcm~?. The source backing 
consisted of a nylon film on to which a thin layer of aluminium had been 
evaporated, giving a total thickness of approximately 0-2 mgcm™. 

Three conversion lines were found at electron energies of 0-0779, 0-1388 and 
approximately 0-150 Mev, the last two being incompletely resolved. These 
correspond to K, L and M conversion lines due to a transition of energy 
0-152 + 0-0005 Mev, a value in good agreement with that obtained by Cork et al. 
The relative intensities of the K and L conversion lines were estimated both 
from the heights of the lines on an undivided momentum plot and from their 
areas on a divided momentum plot (i.e. a graph in which the ordinate P(7) is 
obtained by dividing the counting rate at a given momentum setting by the 
corresponding momentum 7). The mean value thus obtained for the ratio of 
the number of K conversion electrons to the number of L conversion electrons 
Nx/Ny, was 0-42+0-04, By comparing areas on the divided momentum plot 
the ratio of L conversion electrons to beta particles Ny,/N, was found to be 
0-076 +0-008. By the same method a value of 0-032 + 0-004 was obtained for 
the ratio N,/N, and from this value the ratio of K-conversion electrons to 
gamma-quanta N;/Ng was estimated to be 0-046 + 0-01. 

An empirical curve of the ratio Nx/N;, for electric quadrupole radiation 
by Goldhaber and Sunyar (1951) predicts a value of 0-5—0-6 for 48°Os which 
agrees well with the observed value. On the other hand according to the tables 
of Rose et al. (1951) the experimental value of Nx/No agrees best with the 
calculated value for electric dipole radiation and is less than the value for electric 
quadrupole radiation by a factor of 7: 


' 


Electric Magnetic Observed 
1 1 2 3 1 2 3 
Nx/No Op Vi e033 0°95 1-58 8:00 29-0 0:046+ 0-01 


It does not seem reasonable to attribute this discrepancy to experimental error 
in the determination of the intensity of the 0-152 Mev gamma-ray because intensity 
measurements of the 0-136 Mev gamma-ray of 1**Re, using the same experimental 
equipment, have given a value agreeing well with the intensity of the associated 
partial beta-spectrum measured on the spectrometer (Grant and Richmond 
1949; see also Steffen 1951). Scattering and absorption due to non-uniformity 
in the source undoubtedly reduced the number of K-conversion electrons 
observed, but careful calculation shows that this effect could not have reduced 


Coincidences (min-') 


Fig. 3. 
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the K-conversion peak by more than about 15% at most. It would therefore 
appear that there is a real and substantial disagreement between the theoretically 
predicted number of K-conversion electrons and the number actually emitted. 


§4. COINCIDENCE MEASUREMENTS 

Curves of the absorption in lead of y~y coincidences are shown in figs. 3 and 4. 
These were taken using two photomultipliers (E.M.I. type 6262) in conjunction 
with a coincidence circuit of resolving time 2 x 10-8 second. The curve of fig. 3 
was obtained by placing equal thicknesses of lead absorber between the source 
and each multiplier and measuring the coincidence counting rate as a function 
of absorber thickness. Curve A of fig. 4 was obtained by placing different 
thicknesses of lead absorber between the source and one multiplier only. A lead 
absorber of thickness 1-8 g cm~? was then placed between the source and the other 
multiplier and the experiment repeated to give curve B. In all cases— 
aluminium absorbers were placed adjacent to the source to absorb the 
beta-particles. 


Coincidences (min~') 
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Absorption‘jof y-y coincidences (absorbers _ Fig. 4.. Absorption of y-y coincidences (absorbers on 
on both sides of source). one side of source). 


Observations of B-y coincidences were also made, the beta-particles being 
excluded from one multiplier by means of an aluminium absorber. Gamma-ray 
absorption measurements showed that beta-particles were in coincidence with 
the 0-152Mev gamma-ray and also with the harder gamma-rays. Absorption 
of the beta-particles in aluminium showed that particles of energies greater than 
1-5 Mev were in coincidence only with the 0-152 Mev gamma-ray. 

Coincidence measurements made after the decay of the 188Re showed that 
*®6Re radiations could not have affected the above results. 


§5. DECAY SCHEME 
The five gamma-rays may most simply be fitted to the decay scheme shown 
in fig. 5. ‘The deductions which may be made from the experimental results 
are summarized below and it will be seen that they are all in accordance with this 
scheme. 

Intensity Measurements 


(i) The 0-152 Mev gamma-ray has an intensity four times as great as the total 
intensity of the four harder gamma-rays, and must therefore come from the 
lowest excited state. (ii) Since the maximum beta-ray energy is about 2-1 Mev, 
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energy considerations show that the 1:40 Mev and 0-933 Mev gamma-rays cannot 
be in coincidence. Coincidences between the 1-40 Mev gamma-ray and the 
0-638 and 0-476Mevy gamma-rays would correspond to beta-transitions of 
energies about 0-1 and 0-2 Mev respectively. In either case the intensity of the 
1-40 Mev gamma-ray would be very much less than the observed value of 4:5%. 
The 1-40 Mev gamma-ray can therefore be in coincidence only with the 0-152 Mev 
gamma-ray. 
Coincidence Measurements 

In the coincidence absorption curves only three components may be 
differentiated : (a) the 0-152 Mev gamma-ray; (b) a component containing either 
the 0-476 Mev or the 0-638 Mev gamma-ray or both; (c) a component containing 
either the 0-933 Mev or the 1-40 Mev gamma-ray or both. 

(i) The curve of fig. 3 shows only components (a) and (6). This is in 
agreement with the proposed decay scheme since, when the 0-152 Mev gamma-ray 
has been removed, only coincidences between the 0-933 Mev gamma-ray and the 
two gamma-rays of component (0) will be recorded and each of these coincidences 


I88Rp 188Q) 5 
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Energy Levels (MeV) 


Fig. 5. Decay scheme of 1*°Re. 


will have the same absorption coefficient as the lower energy member of 
component (b). (ii) The curve A of fig. 4 can be resolved to show all three 
components. Therefore all three components are involved in coincidences. 
(iii) In curve B of fig. 4 the 1-8 gcm * absorber reduces the 0-152 Mev gamma-ray 
to 4% of its original intensity without appreciably affecting the other 
gamma-rays. Hence this curve shows the presence of coincidences between 
the 0-933 Mev gamma-ray and either the 0-476 Mev or the 0-638 Mev gamma-ray 
or both. (iv) The parallelism of curves A and B shows that the 0-152 Mev 
gamma-ray is in coincidence with both (6) and (c). 

While none of the experimental evidence contradicts the proposed scheme, 
the ambiguity attached to components (6) and (c) makes it possible to construct 
more complicated schemes (containing five energy levels) which would give the 
same experimental results. 
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ABSTRACT. "The Landau (1941) picture of liquid helium II as a modified Debye 
model is critically compared with the London—Tisza picture of liquid helium II as a 
condensed Bose-Einstein assembly in an attempt to find some common features. "Two 
such common features are found, first that the criterion for Bose-Einstein condensation 
in a one-particle model is (0S/0N)pg,y <0 or g=(0F/¢N)y,p >0, a condition that can also 
be satisfied by the Debye model below a certain temperature. Secondly, an assembly 
satisfying such a condition, is in principle, capable of being in thermodynamic 
equilibrium with a weakly adsorbed film or with free atoms of small energy. It is also 
suggested that such atoms may make possible the observed phenomena of superfluid 
flow because their wave-functions are abnormally sensitive to small applied external fields 
of force, and some evidence in support of this is presented. 

A numerical study of the Landau—Debye model shows that, in accordance with 
observed facts on helium II, its energy spectrum is probably modified considerably in 
channels of less than 10~* cm width. It is also concluded that the Debye model is probably 
not a workable approximation near the A-point, but that it is correct at sufficiently low 
temperatures. 


Si INTRODUCTION 

HE behaviour of liquid helium can be quantitatively described by means of 
| the ‘two-fluid’ picture, which is the assumption that liquid helium behaves 

as a solution of two fluids, each with its own density and stream velocity, one 

of them, the superfluid, being devoid of viscosity and unable to carry any entropy, 
while the other component behaves like an ordinary liquid, with a finite viscosity 
that can be measured by an oscillating disc type of experiment (Andronikashvili 
1948). It.is also necessary to postulate that the ‘superfluid’ effects disappear 
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progressively as the velocities increase, and Gorter and Mellink (1949) have shown 
that many of the experimental facts can be interpreted by the assumption of a 
mutual friction between the two fluids and there is also evidence of a critical 
velocity of superfluid flow (Kapitza 1941, Atkins 1951, Bowers and Mendelssohn 
1950). Gorter and Mellink also find some evidence of a decrease in the effective 
viscosity of the normal fluid for very narrow channels (~10~*cm) which they 
suggest may be a mean free path effect. Also, for channels narrower than this 
approximate width, the critical velocity of superfluid flow begins to increase with 
diminishing width of channel, this critical velocity being roughly inversely 
proportional to the width. 

‘Two physical models which lead to the two-fluid theory as a consequence have 
been proposed. F. London (1939) has suggested that the perfect gas model is a 
satisfactory first approximation, and that the appearance of superfluidity is 
connected with the condensation phenomenon that occurs in a perfect gas obeying 
Bose-Einstein statistics. Landau (1941) takes the Debye model of a solid as a 
first approximation, supplementing it by postulating the presence of ‘ rotons’ as 
well as ordinary sound quanta or ‘phonons’. ‘The ‘roton’ is supposed to be a 
quantized state of rotational motion of the liquid. It turns out that the total 
effective mass of these excitations, defined in a special way to be described below, is, 
at very low temperatures, less than the actual mass of the assembly, and in this way 
it is also possible to arrive at a two-fluid picture. The predictions of the two 
theories only differ appreciably at temperatures below 1° and in this region it 
has been shown by Atkins and Osborne (1950) that the measured velocity of pro- 
pagation of ‘second sound’ or temperature waves is in much closer accord 
with Landau’s theory. 

In this paper we shall study a simplified form of the quantum-mechanical 
many-body problem, and we shall show how, as the strength of the interactions is 
progressively reduced, we pass over continuously from a slightly modified form of 
Landau’s model to the perfect gas model. At the same time, interpretations will 
be given of the anomalies associated with flow through very narrow slits noticed 
by Gorter and Mellink (1949), of the influence of zero-point energy on the pro- 
perties of liquid helium, pointed out by Simon (1934) and of the peculiar behaviour 
of the viscosity of the normal fluid as a function of temperature (Andronikashvili 
1948). Finally, we shall point out a possible reason for the occurrence of super- 
fluidity or, more precisely, flow of liquid helium without appreciable transport of 
entropy under very small pressure heads. 


§2. THE RELATIVE SUCCESSES OF THE TWO MODELS 

We now give a list of the experimental facts (a) that do not point decisively in 
either direction, (4) that seem to confirm the Landau type model, (c) that seem to 
confirm the Bose-Einstein type model. Under heading (a) we may first of all 
mention the equilibrium properties such as the equation of state and the specific 
heat curve. Neither of these show anything particularly unusual except for the 
small negative expansion coefficient and absence of a triple point. ‘The anomaly 
in specific heat is of a type very commonly encountered. As Simon (1934) has 
pointed out, there is evidence of an abnormally high zero-point energy (of the order 
of 70 cal/mole) which has the effect of making the lattice spacing abnormally high. 
A theory of Debye type gives a zero-point energy of 9NR@/8 which is of the right 
order of magnitude. ‘Temperley (1947) has pointed out that a generalization 


492 H. N. V. Temperley 


of the perfect gas model, in which each atom is supposed to be moving in a spatially 
periodic field representing its interactions with the other atoms, also possesses a 
zero-point energy and may have a coefficient of expansion of either sign. 

A second set of experimental facts under this heading seems to be those that can 
be interpreted by the two-fluid theory, without enabling one to calculate a definite 
value of the ratio p,/p,. Among these we may mention the thermo-mechanical 
effect, the experiments on heat transport, discussed by Gorter and Mellink (1949) 
and the experiments on flow in capillaries discussed by Atkins (1951). 

A third set of experimental facts, the oscillating disc type of experiment and 
measurements of the velocity c, of second sound above 1°x do lead to definite 
values of the ratio p,/p, (Gorter, Mellink and Castelejn 1950) but the required 
values of this ratio are consistent with either theory. 

Under heading () come the experiments of Atkins and Osborne (1950) on the 
velocity c, of second sound below 1°K. ‘These confirm entirely the results of 
earlier work which indicated a rapid rise in c, below 1° and point to values of c, in 
close agreement with Landau’s theory. As has been shown (Temperley 1951 a, 
Dingle 1951), such results can only be obtained from a theory of London—Tisza 
type if one abandons the implicit assumption that p,/p=N,/N where Np is the 
occupation of the lowest state. A second piece of evidence is the steep rise in the 
viscosity of the normal fluid with falling temperature in this region, which is, as. 
we shall see, characteristic of a ‘solid-like’ rather than a ‘ gas-like’ model. 

Under heading (c) comes the occurrence of a transition point. According to- 
Landau (1941) the transition to HeI should occur when the effective density p, as- 
defined by him becomes equal to the actual density p, but his theory in its present 
form gives no explanation at all of why this equalization should be accompanied by a 
co-operative transition. Such a transition is, however, a natural consequence of 
the Bose-Einstein type of theory (F. London 1939). A second piece of evidence 
that a ‘gas-like’ model may be more appropriate near the transition point is the fact 
that the viscosity is rising with increasing temperature in this region (as well as for 
He lI). 

What might seem to be a third piece of evidence is that the formation of thick 
and mobile adsorbed films below the transition temperature receives a natural. 
explanation (Temperley 1949, Band 1951), if one supposes that the adsorbed 
film is in equilibrium with a condensed Bose-Einstein assembly. The Landau 
theory, in its present form, gives no such explanation, but we shall see that if 
properly developed the Landau theory could explain this fact; thus this piece of 
evidence really belongs in category (a). 

‘Taken all together, the evidence appears to indicate that the Debye-type model 
is a good first approximation at very low temperatures, but that a ‘ gas-like’ model 
may be more appropriate near the transition point. The behaviour of the viscosity 
of the normal fluid, which apparently passes through a flat minimum at about 
1-7°k, also suggests that more than one mechanism is responsible for it. 

The Debye model will now be examined quantitatively and it will be shown 
that it probably is a good approximation at very low temperatures. The manner 
in which it passes over into a ‘ gas-like’ model as the interactions are weakened or as 
the temperature rises will also be traced. A study of its entropy, as a function of 
the number of atoms present, shows that the model is theoretically capable of a 
‘condensation effect’ which has not been noted before. This seems to provide a 
connecting link between the two types of model, and to give another instance of the 
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already known fact, see §5, that the Bose-Einstein condensation phenomenon 
is not confined to ‘ one-particle’ models but that analogues can exist in the many- 
body model even with very strong interactions. On this basis, the formation of 
thick adsorbed films can be understood on the basis of the Landau model (1941) as 
well as on the London—Tisza model (F. London 1939). 


§3. THE DEBYE TYPE MODEL 

It is known (see, for example, Peierls 1929) that the Debye model is practically 
equivalent to the assumption that the potential energy of a crystal lattice is a 
quadratic function of the distances of the atoms from their equilibrium 
positions, or from one another. Peierls uses the linear lattice to illustrate the 
departure of the modes of an actual atomic lattice from the ideal Debye spectrum. 
We may write V =3mw,?X,(x,—%,,,)? and obtain for the classical equations of 
motion ay 


dt 


By making the substitution €,=%,.x,exp(2rs71/N) Peierls finds for the possible 
modes of oscillation 


Bea (Lee ine a ety) — Oia at at) 2) AUS (1) 


NE VENTE CORD a ee Rie (2) 

These differ from the ideal Debye spectrum in that the possible frequencies are 
proportional to sin (s7z/J) instead of to s, which means that one automatically gets 
N normal modes without having to assume an arbitrary cut-off of the spectrum as 
Debye does. If we write the quantum-mechanical version of eqn. (1) we get 

we en FER PES (a, ep PAO. (3) 
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To avoid end effects, we may suppose fh atoms arranged in a ring of total length / 
and that the last atom interacts with the first, the interaction between them being 
proportional to (x, —*y+/)? rather than to (x,—xy)?. ‘This secures that the total 
energy is a minimum when the atoms are spaced equally round the ring, instead 
of when they are all bunched together. Symmetry between the N atoms can then 
be restored if we write y, = x1, V2 =, —l/N, y3=x,—2l/N...etc. This makes the 
potential energy again a quadratic function of the coordinates, and the terms 
involving w,7/? simply ae a change of the SHS zero. Equation (3) becomes 
separable if we write C,=%,y,cosrs0, S,=2,¥y,sinrs@ (where 0=27 |). It 
becomes 


Nh [32 2? Ima? . 5 (8 mi Went 
| ptm (sc3 + aga) +E Pe Wupoge a (Fic. Sabai 


A typical solution of eqn. (4) is 
Oras Ff Gy S, 
b= Mexp(— ate) (Si (=), oe ceee (5) 


$ s 


‘E 


where the H,, H,, are Hermite polynomials, not the same for all values of s, the 
order of each of these polynomials indicating the number of sound quanta or 
‘phonons’ associated with the corresponding variable. ‘he energy associated 
with the variable C, or S, is a half-odd multiple of hw,, where w, is given by eqn. (2), 
and r, is a distance given by (Nh/47mw, sin $86)", and is a measure of the range 
through which C, or S, can depart from the zero values corresponding to uniform 
spacing of atoms round the ring before the wave function ¢s becomes negligible. 
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We can obtain a relation between the constant w, specifying the strength of the 
interactions and the measured velocity of sound c, at low frequencies by an 
application of the correspondence principle. As Peierls (1929) remarks, in a 


problem of this kind the wavelength cannot be strictly defined, and its place is 


taken by the parameter s, which specifies the phase difference between neighbour- 
ing atoms associated with phonons involving the variable C, or S,. For small 
values of s this phase difference is practically the same as that associated with a 
sound of wavelength //s, so that we may deduce the relation between c, and a» 
qa d= 5 Ftsin (SS) = ae. . saa (6) 
If we wish to preserve the acoustic concept of a wavelength //s for large values of s, 
it is necessary to admit that the corresponding velocity of sound depends on s. 
The changes necessary if we work with an atomic lattice instead of a continuum are 
an upper limit to the possible frequencies, given according to eqn. (2) by w/z, and 
a dispersion in the velocity of sound. Substituting for w, in the relation for 7, we 
get 

r oN Ala me, sy)" <.  &  e e (7) 
which enables us to estimate the distance y, by which a single molecule is shifted 
from its equilibrium position as a result of the zero-point motion associated with a 
quantized sound wave of wavelength //s. For small values of s, the corresponding 
classical sound wave would be represented by a space variation of the type 
V,=YC0S(rs), Yo being the amplitude. Comparing this with the expression 
C,=,¥,cos (rs?) and the condition that C, may attain values of the order of 7,, 


we conclude that yor 2r,|Na(hl)Nrmesy?. ae (8) 


A result very similar to (8) follows if we suppose that the energy $hv, is shared among 
N molecules. Taking the case s=1 and inserting numerical values 
l/N =3-6 x 10-8§cm, c=2-5 x 10*cm/sec and m=6 x 10-**g appropriate for liquid 
helium, we find yy~10-8cm; thatisto say that the zero-point motion has amplitudes. 
comparable with the atomic spacing, a result that has been pointed out by several 
writers. 


§4. THE LIMITATIONS OF THE DEBYE MODEL 

We have to consider first the changes that have to be made in the one-- 
dimensional Debye model to make it apply to a three-dimensional lattice, secondly 
the manner in which it is likely to break down as the interactions are reduced or the 
temperature raised, thirdly the modifications necessary to allow for the symmetry 
properties of the helium atoms, and lastly the effects of cubic and higher terms in 
the energy. 

(i) The Change to the Three-Dimensional Case 

This question has been considered by Peierls (1929) and, in more detail, by 
Blackman (1935). As Peierls shows, the introduction of interactions between 
other than nearest neighbours in the linear model introduces no new feature of 
physical interest while, in the two and three dimensional cases, the introduction of 
‘ cross-bracing’ interactions between second nearest neighbours implies the possi- 
bility of propagation of transverse or shear waves. ‘Temperley (1951 b) has. 
suggested that Landau’s ‘rotons’ might well be replaced by quantized shear waves 
which do represent rotational motion, and that the available evidence is consistent 
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with their velocity of propagation being comparable with, perhaps even somewhat 
greater than, the ordinary velocity of sound. The introduction of these extra 
interactions makes the energy spectrum more complicated but as long as we confine 
ourselves to terms in the interaction energy V that are quadratic in the coordinates, 
the equation can still be separated. One can always find the orthogonal trans- 
formation of coordinates that reduces such a function V to a sum of squares, and 
such a transformation changes ©, A,? into the same operator in the new coordinates, 
Associated with each of these coordinates, we can define an ‘effective velocity of 
sound’ by a process analogous to that leading to eqn. (6). The correspondence 
principle shows that the continuum treatment must be accurate at sufficiently low 
frequencies but even here we may have to allow for both longitudinal and transverse 
oscillations. We thus arrive at the ordinary Debye model with zero-point 
energy $/w, associated with each normal coordinate, the longitudinal and transverse 
velocities of sound being in general different, and both varying with frequency. 
The dispersion is a consequence of the fact that we are dealing with a discontinuous 
lattice, and is already present in the one-dimensional model. 

The physical picture of a phonon arrived at above differs in some respects from 
the rather vague picture of a ‘wave packet of sound’ (localized in space 
and travelling with the group velocity), that most writers seem to have in mind. 
Our assembly is to be described by a superposition of states ¢ (t,, u,...) given by 
(5), where the #’s and w’s take all possible values, but, in practice, a complexion 
in which ¢, departs appreciably from its statistical mean value 1/{exp(hw,/RT) — 1} 
is very unlikely, so that the assembly is probably described sufficiently accurately 
by a pure wave function in which the ?’s and w’s have values near to these mean 
values. If so, the emission or absorption of a phonon of frequency w,/27 corre- 
sponds simply to a change in ¢, oru, of unity, but such a process affects all the atoms 
on an equal footing, and cannot therefore be associated with any definite part of the 
assembly. We cannot meet this difficulty by describing the assembly by a mixed 
rather than by a pure wave function without blurring the concept of emission or 
absorption of a phonon (we cannot get a clear wave-mechanical picture of the 
absorption of a quantum of radiation, unless we begin with an atom in a definitely 
specified state). We can, however, use Peierls’ device (1929) of dividing up our 
assembly into L sub-assemblies of N/L atoms, and then assign a set of wave 
functions to each sub-assembly. ‘The disappearance of a phonon from one 
sub-assembly and its appearance in another can then be given a definite meaning, 
but it remains true that each phonon must be associated with a large number of 
atoms and therefore that it cannot be precisely localized. 

Another important point is that the energies and frequencies of phonons are 
quite unaffected by a uniform external field such as gravity. For the addition of a 
potential energy mg(y,+2+...) to eqn. (4) affects only the equation involving 
Cy=(¥,+32+.-.) which governs the motion of the centre of gravity of the 
assembly, and none of the equations involving C, or S, (s40). (A non-uniform 
field, involving something like X,y," would introduce stresses into the assembly, 
and this could, of course, affect the frequency spectrum.) We can also show that 
the energy of the phonons is unaffected by a Galilean transformation of coordinates. 
For, if we put y,’=y,+ vt, t'=t, we have 0/dy,’=0/d0y, 0/0t=0/dt' + vu,d/dy, 
so that the change in energy of the whole assembly(h/277)(0/dt' —d/dt) is 
(vh/27i)»,,0/0y,’ =(vh/2ni)0/0Cy, which energy change is all attributable to the 
change in velocity of the centre of gravity of all the atoms, and there is mo energy 
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change associated with the phonons themselves. Wave functions involving C, or 
S, actually correspond classically to stationary sound waves, but the above result 
remains true even if we use the quantum analogue of a travelling sound wave 
(see for example Wilson 1936, chap. 6). Landau’s (1941) derivation of the effective — 
mass of an assembly (see p. 503) seems to be based on the assumption that phonons 
behave like material particles and that the energy of each one is changed by pV by a 
Galilean transformation of axes. However the change in energy of any assembly 
of material particles, such as a Debye lattice, must , as we have seen, be proportional 
to their total momentum, that is, to the momentum of the centre of gravity, so 
that Landau’s derivation needs revision. Nevertheless, it is possible to obtain 
results, precisely equivalent to Landau’s for the effective mass of an assembly by a 
somewhat different argument, which will be found in the Appendix. In this 
revised derivation the ‘ momentum of a phonon’ is defined simply as the momentum 
that would be imparted to a material particle if the phonon disappeared, and nothing 
is assumed about any change of energy of an individual phonon as a result of a 
Galilean transformation. 


(ii) The Effect of Raising the Temperature 

At extremely low temperatures, the energy associated with the normal 
coordinate C, is $iw, where w, is given by eqn. (2) and, as we have seen, the 
corresponding migration of an atom from its equilibrium position associated with 
one particular frequency is of the order of 10-§cm. If, however, we raise the 
temperature, the probable energy associated with C, is hw,/{exp(hw,/RT) — 1} =RT 
forkT>hw, and the corresponding migration of an atom is multiplied by a factor 
of the order of (RT/hw,)"”. ‘This factor will be largest when s=1 and, as we have 
seen, this value of s already gives the largest migrations at absolute zero. Thus, 
according to the Debye theory, at 1°x the thermal vibrations of the lattice are 
already displacing individual atoms through distances of the order of 
(RT/hw,)'? x 10-§=10->cm, and the conception of liquid helium as a crystal 
lattice with the atoms bound in the neighbourhood of equilibrium points is 
apparently becoming far-fetched. Nevertheless, at least one other instance is 
known of a bound quantum-mechanical system in which both the potential and 
kinetic energies are large compared with the actual binding energy, namely the 
deuteron. The characteristic feature of this system is that the proton and neutron 
are very likely to be at a distance apart which is significantly larger than the range of 
action of the force between them, but nevertheless they continue to form a bound 
state. Furthermore, the behaviour of such a system is not very sensitive to the 
precise form of the interaction as a function of the distance, but only to various 
weighted averages of this function. In an analogous way, it seems possible to 
understand why in liquid helium the crude assumption of a square law form of 
interaction (which is certainly not valid for an atom that has migrated as much as 
one lattice distance) can nevertheless give results in reasonable agreement with 
experiment. 

‘To follow the changes that take place as the temperature rises or as the inter- 
actions are weakened, we should have to assume somewhat different boundary 
conditions in eqn. (4). So far, we have implicitly assumed those appropriate to 
the ordinary simple harmonic oscillator problem, namely that the wave functions 
must represent bound states, and must thus vanish sufficiently fast as any of the 
C.>o. If the interactions were weaker, we should have to take account of the 
effect of the walls, and the problem would become similar to that of simple harmonic 
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oscillators enclosed between rigid walls, which has been treated by Auluck and 
Kothari (1945). The effect on the energy spectrum would become noticeable at 
“ordinary’ helium temperatures when the walls were at distances of the order of 
10-°cm which may be the explanation of the fact, noted by Gorter and Mellink 
(1949), that the calculated viscosity of the normal fluid apparently falls off sharply 
for slits of and below this width, and of the increase in the critical velocity of super- 
fluid flow for such slits. The calculated figure of 10->cm for the migration of a 
typical atom associated with non-zero values of C, or S, seems at first sight rather 
small, both the above effects already being quite noticeable for slits of 10-4cm width. 
We must remember, however, that this calculated figure represents the migration 
of an atom associated with one frequency only, whereas there will be a displacement 
of an atom associated with each of the normal modes. These displacements have 
to be added up in a random manner because, in a solution such as (5), the phases of 
C,and S, are unspecified. In this way it seems easily possible to account for the 
discrepancy of a factor of 10 or so. 


(iii) The Symmetry Properties of the Atoms 

So far. we have taken no account of the fact that our theory is to be invariant 
if we interchange one or more pairs ofatoms. The interaction energy in an equation 
such as (3) ought to be a symmetrical function of the coordinates, whereas we have 
neglected all interactions except those between nearest neighbours. This could 
be formally remedied by replacing the interaction function by one such as 
x ,&,V4(x,—x,)? (which would actually lead to a model of Einstein rather than 
Debye type, all the normal modes having the same frequency). Such a model, 
however, implies that the interaction between distant atoms is far greater than that 
between nearest neighbours, and it is, for that reason, unphysical. What we 
really require, to represent a solid, is an interaction energy that is always large 
except for configurations corresponding to a nearly perfect crystal lattice. 

For such a model, it is shown by Rushbrooke (1949) that the indistinguish- 
ability of the atoms does not matter, and that we can get the correct 
partition function even if we make the physically inadmissible assumption that the 
atoms are arranged on the lattice in a fixed numerical order. In quantum- 
mechanical terms, we obtain the same set of energy levels whether we use wave 
functions such as (5), or the correct ones obtained from these by the symmetrizing 
process of permuting the atoms and superposing the resulting wave functions. 
As we weaken the interactions, exchanges of places between the atoms become more 
_ and more frequent, and it becomes a poorer approximation to use wave functions 
such as (5). Finally, when the interactions become very weak, we are back to the 
perfect gas model, for which the symmetry properties of the atoms are very 
important. Thus, we reach Wergeland’s conclusion (1947) that we have to 
consider the symmetry properties of the atoms when, and only when, the inter- 
actions are weak. 


(iv) The Effect of the Cubic Terms in the Energy 

Peierls (1929) showed that, even in an ideal lattice at very low temperatures, it is 
quite essential to consider the effect of the cubic terms in the interaction energy, 
otherwise the thermal conductivity of the lattice diverges. ‘This remains true 
even if we postulate the existence of regions in the lattice where the effective value. 
of the velocity of sound is different from that in the surrounding regions. 
Scattering proportional to v* occurs, but the thermal conductivity still diverges 
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because of an ‘infra-red catastrophe’. A similar result is easily seen to be true for 
viscosity. Taking both the irregularities and the cubic terms into account, 
Peierls (1929) considered three-phonon ‘ Umklapp’ processes and concluded that 
the conductivity at very low temperatures should behave like 1/7. Landau and 
Khalatnikov (1949) considered four-phonon processes without ‘Umklapp’ and 
concluded that the viscosity due to phonons only should behave like T~°, and they 
also found it essential to take dispersion into account. 

This type of variation of the transport coefficients with temperature seems to be 
characteristic of the many-body problem. In contrast we have the results of 
calculations of ‘ gas-like’ viscosity and thermal conductivity, when only two-body 
collisions are considered but the cross sections are calculated quantum mechani- 
cally (see, for example, Buckingham, Hamilton and Massey 1941). These give a 
viscosity which generally increases with rising temperature except near absolute 
zero, where the results depend greatly on the precise form of the interaction. The 
processes considered by Peierls (1929) (three-phonon Umklapp) and Landau and 
Khalatnikov (1949) (four-phonon) always involve the emission or absorption of 
small numbers of ‘ phonons’—processes which involve the re-arrangement of a 
very large number of atoms (theoretically all the atoms in the assembly are affected). 
Conversely, a process such as the collision or interchange of two atoms only, would 
be described in the Debye picture as a process involving all the normal coordinates. 
C, and S,, and therefore also affecting a very large number of phonons. The 
relation between C,, and x, is reminiscent of that between a variable and its Fourier 
transform. A process involving a single C, involves all the x,’s and vice versa. 

It therefore seems fair to regard the ‘two-body’ and ‘ Debye’ contributions to 
the transport coefficients as fundamentally distinct from one another, and to 
suppose that they can, in first approximation, simply be added together. If this is 
so it seems that, in liquid helium II, the two types of process are of comparable 
importance. The ‘ gas-like’ two-body calculations by themselves give a viscosity 
of the order of magnitude observed in the liquid. The ‘Debye’ contribution 
to the viscosity is very difficult to estimate numerically, and one can predict only 
its probable variation with temperature. To sum up, if we wish to use a model of 
Debye type, we have to admit the possibility of variations of the velocity of sound 
c, with density and temperature (these variations are experimentally observed). 
In addition we have to allow for the atomicity of the crystal lattice, which, as we 
have seen, implies a variation of c, with frequency. It is also likely that the 
velocity of propagation of transverse sound waves is not a constant. Thus, near 
the transition temperature, where the variation of c,; with temperature begins to. 
become rapid (Atkins and Chase 1951) a model of Landau—Debye type would 
contain so many adjustable constants as to be almost useless, quite apart from the 
fact that the large migrations of the atoms perhaps make it a poor approximation in 
any case. 


§5--A GENERALIZATION OF THE MECHANISM OF 
BOSE-EINSTEIN CONDENSATION 
Haselgrove and 'Temperley (1952) have investigated the statistical mechanics of 
Bose-Einstein condensation from the standpoint of the analytic theory of partition 
of numbers, one of their conclusions being that a sufficient condition for the 
occurrence of ‘condensation into the lowest level’ is probably that the thermo- 
dynamic probability W or py (Z), the number of ways in which the energy E can be 
shared among N atoms, is amaximum for variable N and fixed E. By Boltzmann’s 


Modes of Liquid Helium II 499 


hypothesis, we have S = const + k log W so that the criterion for the appearance of a 
condensation phenomenon is 


os oF 
(ar), : <0 or tsa) 4 10): weeees (9) 


Consider first a model of the same type as the perfect gas, in which the energy 
levels refer to the individual atoms, the so-called ‘ one-particle’ type of model. At 
ordinary temperatures one can always gain entropy by adding another atom to the 
assembly at constant volume, even if no energy is added with it. This follows 
from the principle of equipartition of energy, valid for ‘classical’ assemblies. As 
the total energy EF is reduced there will, however, come a stage when the number of 
ways of partitioning this energy into N non-zero parts among N atoms no longer 
increases with N. For N greater than this critical value, one can actually gain 
entropy by assigning some of the atoms no energy at all—that is to say, the most 
probable state will be one in which a finite fraction of the atoms are in the very 
lowest energy level of the assembly—the energy of which level we can always 
choose to be zero. 

Haselgrove and Temperley (1952) have shown that this phenomenon probably 
occurs for a very wide class of ‘ one-particle’ models and that it is almost certainly 
identical with Bose-Einstein condensation. For the perfect gas it can be formally 
proved that the known formulae for the energy and entropy lead to the well-known 
formula for the transition temperature 


N h2 Bet oO aie 
> (scaer) Ba LO ita a tae (10) 


as a result of imposing the condition (9) (0S/AN), » <0. The formula (10) is. 
usually derived from the consideration that the occupation of the lowest level must 
be finite, but condition (9) can be applied to any assembly, whereas considerations 
depending upon occupation numbers can usually be applied only to models of 
‘one-particle’ type. 

We therefore postulate condition (9) as the criterion for what we may call 
generalized Bose—Einstein condensation. Its consequence, that particles can be 
removed from the assembly or relegated toa state of zero energy without decreasing 
the total entropy, can also be exhibited in another way, which may make its physical 
meaning clearer. Condition (9) is equivalent to 

BORON rem 20; (ial wf Shin Sadakes (11) 
where Fis the Helmholtz free energy, which states that if a particle can be reversibly 
removed from such an assembly under isothermal conditions and put in a state of 
zero free energy no mechanical work need be done. Naturally, the thermo- 
dynamic treatment gives no information on the possibility of such a transfer, it 
merely draws attention to the existence of two types of assembly. Ina ‘classical’ 
assembly, the introduction of a further particle without energy would bring the 
principle of equipartition of energy into operation with a resulting increase of 
entropy. In a ‘condensed’ assembly the energy is not shared and no such 
entropy increase occurs. ‘I'emperley (1949) has called attention to the fact that if 
an adsorbing surface is in equilibrium with a condensed Bose—Einstein assembly, 
an adsorbed film much thicker than usual may result. ‘The reason for this is. 
that it is thermodynamically profitable for the film to build up as long as the energy 
of an adsorbed atom in the top layer is less than that of the lowest state of the Bose— 
Einstein assembly, whereas, in ordinary conditions, one needs an adsorption 
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energy of several times RT to compensate for the relatively small number of ‘film’ 
states. (The same point has recently been made by Band (1951) without any 
acknowledgment.) 

We shall now show that the Landau-Debye model of liquid helium also 
possesses the property (9) or (11), and, in fact, that there is a temperature below 
which the signs of inequality hold. This, as we shall see, results in the ‘ osmotic 
driving force’ that is required to explain second sound and the thermo-mechanical 
effects. The perfect gas model itself is a limiting case and the signs of equality 
hold below the transition temperature, and the predictions of the thermo- 
mechanical effect and velocity of second sound follow only as limiting cases of 
the theory to be outlined below. In eqn. (13), for example, g is zero below the 
transition temperature, so that (0g/0P)7 and (dg/0T)p are both zero. 

For the Debye model at low temperatures, we have for the Helmholtz free 
energy F=9 NkO@/8 — 4Nx*kT(T/©)3 where N is the number of atoms present and 
© the Debye temperature. By successive approximation we find that eqn. (11) 
is satisfied for temperatures below T, 0-750 if we neglect any variation of © with 
density and any change in the zero of energy such as the term V4 in eqn. (4) and 
thus, at a temperature of this order, the chemical potential g changes sign, becoming 
positive below it. Such a change of sign has never before been thought to have 
any particular physical significance, in view of the fact that energies in thermo- 
dynamics are often measured from an arbitrary zero. It does, however, seem to 
play a part in the theory of liquid helium. In actual helium, the energy Ep at 
absolute zero is probably negative (F. London 1939), the shift in the energy—zero 
due to the attractive forces being a little greater than the zero-point energy. 

Tisza {1947) has proposed as a model of liquid helium II a Debye-like lattice in 
which the shear modes have undergone Bose-Einstein condensation, but he has 
given no actual proof of the possibility of such an effect. However, another 
many-body model, in which Bose-Einstein condensation would also occur, 
results from Bloch’s ‘spin-wave’ treatment of ferromagnetism at low temperatures 
(1930). It has been remarked by many writers that this model is mathematically 
equivalent to the perfect gas at sufficiently low temperatures, so that it shows a 
transition of the sametype. (‘This is a rigorous consequence of the assumptions on 
which the model is based; whether these are physically valid near the predicted 
transition temperature is another matter.) Thus, it is by no means a new sug- 
gestion that Bose-Einstein condensation may have analogues in the many-body 
type of model, and the possibility has been rigorously proved in at least one instance 
besides the one now being examined. 

Consider a Debye lattice containing N atoms, and suppose it to be in thermo- 
dynamic equilibrium with an adsorbing surface containing M sites, each of 
energy. For equilibrium we may suppose that m atoms have left the lattice and 


are adsorbed on the sites. ‘The entropy of m atoms on M sites is k log ea) so 
that the free energy of the m adsorbed atoms is i 


Fas =me —RT log = : 


Adding to this the free energy of N —m atoms in the Debye lattice, and differenti- 
ating with respect to m, the condition for equilibrium is 


kT log (57 +e-g=0 ie ioe ee (12) 
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If, say, half the adsorbed sites are occupied, the first term in (12) will vanish, and if 
(¢F/ON)y, p org is negative, then « must also be negative in order to get appreciable 
adsorption, whereas if it were possible to have an adsorption sites with nearly zero 
energy, they would begin to fill up as soon as g became positive. Expressed 

differently, an adsorbed film in equilibrium with a Debye lattice may be expected, 

below the temperature estimated above, to build up layer by layer to a thickness. 
such that the adsorbing energy of an atom in the top layer is negligible. 

This picture is in many respects similar to that of an adsorbed film in 
equilibrium with a condensed Bose-Einstein assembly discussed by 'Temperley 
(1949). In contrast with the perfect gas, the simple Debye model does not 
predict any definite transition temperature below which the film should suddenly 
begin to thicken, but simply suggests a slow build-up in thickness as the tempera- 
ture falls. ‘The Debye model is, however, quite capable of being modified to give 
a co-operative effect and make g a non-analytic function of temperature. There is 
clear evidence (Atkins and Chase 1951) that the approach of the transition 
temperature, from above or below, is heralded by a rapid fall in the measured 
velocity of sound c,. If this situation is to be properly handled by a Debye-like 
’ theory, we shall have phonon energies that depend on temperature and the model 
must then be handled by the special technique described by Rushbrooke (1940). 
_ We do not attempt this at present, as there are few data on the form of the energy 
spectrum, while the particular feature of the interaction between helium atoms. 
that is responsible for this rapid variation of c, with temperature has not yet been 
identified. A similar remark applies to the effect of the variation of cy or wy with 
density. There is a corresponding variation of the zero-point energy with 
density, which, as many writers have pointed out, can be described as a ‘ zero-point 
pressure’ which is responsible for ‘ blowing out’ the liquid to an abnormally low 
density, but we are not yet in a position to predict this variation in detail from the 
interaction function. 


§6. THE FOUNTAIN EFFECT 
Consider now two vessels of liquid helium connected only by an adsorbed film, 

which we may suppose to be at a temperature that is the mean of that in the two 
vessels. Equation (12) shows us that the condition for the helium in either vessel 
to be in equilibrium with the film is g=e,). where €,j). is the adsorption energy 
of a layer that is half occupied (M=2m). (This condition is independent of the 
temperature of the film itself.) Both vessels will then be in equilibrium with the 
film when g is equal in both vessels, which means that the pressure and temperature 
differences are related by 

AP (@g/T)p 

AT. (n/a), sof) 8 He ts elena bare g er cro) 
an equation that was originally derived by H. London (1939) by a thermodynamic 
cycle, assuming the existence of the film, which is capable of transferring matter 
without any appreciable energy transfer. 


§7. SECOND SOUND 
(i) The Fourth Propagation Equation 
The hydrodynamics of second sound as formulated by Landau (1941) involves, 
besides the two ordinary equations of hydrodynamics, the assumptions that all of 
the entropy is carried by the normal fluid (the third propagation equation) and. 
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that the acceleration of the superfluid is proportional to the space-gradient of g 
(the fourth propagation equation). The third equation can be accepted as 
reasonable—it is in accord with the experimental fact (Kapitza 1941) that no entropy 
transfer occurs when helium flows through a very narrow channel. It also follows 
as a natural consequence of Landau’s definition of p,, see below (p. 503). The 
fourth equation is based on the thermodynamic result that, if one atom of helium 
is removed from a given region at constant volume and entropy, the mechanical 
work done is (@E/0N)y,5=(0F/ON)y,7=g but the difficulty has been to explain 
why such a process can apparently be carried out reversibly only in liquid helium 
so that one can get an acceleration of superfluid proportional to gradient g. ‘The 
appeal to the analogy of osmotic pressure does not really help to solve this funda- 
mental difficulty, because, although ordinary osmotic pressure is governed by a 
similar thermodynamic theory, the flow which occurs in response to a difference of 
osmotic pressure is invariably of ‘diffusion’ type and is thus irreversible, the 
velocity of flow being proportional to gradient g. 

The discussion we have just given seems to provide a clue to the answer. We 
have seen that helium atoms of very small energy, or forming part of a weakly 
adsorbed film, can be in thermodynamic equilibrium with a condensed Bose— 
Einstein assembly, or with a Debye lattice below the temperature at which g =0, 
and the movement of such atoms from place to place could then give rise to the 
transfer of matter without any appreciable transfer of energy, the very process 
envisaged by all theories of second sound and the fountain effect. If, however, 
the temperature is raised so that g becomes negative such atoms of low energy can 
no longer be in equilibrium with the assembly, instead the assembly shares some 
of its energy with them so that they become incorporated in it, and the process then 
described becomesimpossible. The idea of special states of low energy, ‘ z-states’, 
was put forward by Mendelssohn (1945). Until now, the only known mechanism 
by which such low-energy states can be occupied below a certain temperature but 
practically empty above it is the Bose-Einstein condensation. However, the 
associated phenomenon of the entropy becoming a maximum as a function of NV 
is applicable to many body models as well, and it is put forward here as the natural 
generalization to such models of the Bose-Einstein condensation effect. Clearly, it 
is not by itself a complete explanation of superfluidity—besides accounting for the 
-occupation of the ‘z-states’ it is still necessary to show that the atoms in them can 
flow from place to place under a very small pressure gradient. Evidently the 
Rollin film, in which the top layers of atoms are weakly adsorbed and probably 
therefore mobile, constitutes one type of ‘z-state’, but the existence of second sound 
‘seems to require the existence of ‘z-states’ in the bulk liquid as well. I hope to 
make a suggestion about their nature in another paper. (The tempting hypothesis 
that the superfluid effects in liquid helium are all due to a ‘sub-surface film’ 
broadly similar to the Rollin film is untenable on various grounds, the most 
important of which are the fantastically high velocities of film flow that would be 
required to explain the existence of ‘second sound’ of quite moderate frequencies 
and the observed fact that the measured velocity of second sound is independent 
of the shape of the vessel (D. V. Osborne, private communication).) However, 
is quite unnecessary to postulate, for example, that at a temperature at which p,/p 
is 1%, 99% of the atoms in the assembly are in 3-states. ‘The fallacy of supposing 
that p,/p is necessarily equal to (N— N,)/N where N, is the number of atoms in the 


2-states, has been pointed out by Dingle(1951). It is true for the perfect gas model, 
but is not true generally. 
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The fact that the measured chemical potential of liquid HeII is decidedly 
negative ~ 10—15 cal/mole shows that the above model is still too specialized. 
It is hoped to consider this point in another paper. 


(11) Landaw’s Concept of Effective Mass 


Landau (1941) showed that it is possible to assign an ‘effective mass’ to an 
assembly of phonons by the following process. One imagines the assembly in 
thermodynamic equilibrium with a vessel the walls of which are moving with a 
small tangential velocity v. Relative to the walls, a difference of energy is 
postulated between a phonon moving in the same sense as v and an exactly 
similar one moving in the opposite sense, the difference in energy being assumed to 
be 2p,v, where p, is the x-component of momentum of the phonon. Thus if 
the assembly of phonons is in equilibrium with the moving walls, there is a net 
momentum of the whole assembly. Averaging over all angles and energies, 
Landau shows that the net momentum is proportional to v for v small, and he 
thus defines an ‘effective mass’ of the assembly which he identifies with the mass 
apparently ‘carried round’ by a disc oscillating in liquid helium. This concep- 
tion of ‘effective mass’ is reminiscent of the effective mass of an assembly of 
electrons in metals, but is definitely not equivalent to it. Landau defines effective 
mass as (mean momentum/Vvelocity), whereas the effective mass of an assembly of 
electrons is defined, see for example Seitz (1940) as 


Total force acting on electrons 
Rate of change of electron flux’ 


The latter definition is inapplicable to phonons, which are unaffected by any linear 
external field acting on the atoms. See p. 495. 

So far, Landau’s theory may be regarded as an attempt to predict theoretically 
the apparent mass of liquid ‘ carried round’ by an oscillating disc. His subsequent 
identification of the difference between this mass and the ordinary inertial mass of 
the liquid helium itself as the ‘effective mass of superfluid’ is much less easy to 
understand physically. Consider, for example, a toroidal vessel of liquid helium 
subjected successively to the processes of pushing the whole liquid forward by 
means of a piston, and then dragging the normal fluid back an equal distance by 
setting the walls in motion. Landau would interpret the final effect of these two 
processes as equivalent to the motion of superfluid alone, and it is easily shown that 
the entire entropy of an assembly of Debye waves is carried along with the walls 
if they are in thermal equilibrium with the assembly (Dingle 1951). Landau has 
therefore shown that this process of superfluid transfer would occur, if at all, 
without any transfer of entropy, and should thus be thermodynamically reversible, 
but he has not given any physical picture at all of how such flow is to take place. 


(iii) A Suggestion about the Possible Nature of Superfluid Flow 

We have suggested above that superfluid flow takes place through the medium 
of z-states which can, in principle, occur in any assembly. Since we do not know 
the nature of z-states associated with a Debye-like model, but have simply shown 
that if they occur, they will be occupied only below a certain temperature, we go 
back to the perfect gas model, the nature of whose ‘ z-state’ is precisely known. 
‘The experimental fact that has to be explained is not a ‘persistent velocity’, 
analogous to the persistent current in superconductors, because no such effect 
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has yet been observed, but the fact that superfluid flow, in films or bulk liquid, can 
take place, at a measurable rate, under extremely small pressure heads, of the order 
of a fraction of a millimetre of liquid helium or less (Daunt and Mendelssohn 
1946, Bowers and Mendelssohn 1950). In this connection, it is of interest to call” 
attention to a feature of the condensed Bose-Einstein gas that does not seem to have 
been explicitly noticed before, namely its extreme sensitivity to an applied external 
field such as gravity. It is fully realized that, in practice, this picture must be 
considerably modified by atomic interactions, but it seems not unreasonable to 
think that it has some relation to the real situation. 

The energy levels of a material particle in a constant external field were found 
by Breit (1928). If the field is along the z-direction, we have to solve the equation 

anaes 
SE + (B— megs + myx) =0 SiN ep (14) 

the choice of the zero of z being arbitrary, and the factors of # involving x and y 
being the same as in the perfect gas. The relevant solutions of this equation are 
oscillatory when the kinetic energy is positive, but fall away rapidly to zero as soon 
as mg(z—2%) >. Insertion of numerical values shows that, for helium atoms, 
the characteristic distance is only of the order of 10-*cm. Thus, if we choose the 
origin of z at the bottom of the container which we may suppose to be of dimensions 
of the order of 1 cm, the wave function for the lowest state will vanish at z=0, 
attain its maximum for z~10~4cm and tend rapidly to zero for larger z, the corre- 
sponding factor of % being very nearly an Airy function. The next highest wave 
function will pass through two maxima as a function of z, and so on. 

The interesting point is that, for particles of the mass of helium atoms, the 
lowest states are very much affected by a field of the order of the gravitational 
field. Without the field, the lowest state would have a factor sin (72/h) where / is. 
the height of the container, whereas the presence of gravity has resulted in this. 
state being of negligible amplitude except near the very bottom of the container. 
For an uncondensed gas, the probability of occupation of each state is very small,. 
and the result of averaging over all states would lead simply to the ordinary 
sedimentation effect that occurs in any gas under gravity, but in a perfect gas below 
the Bose-Einstein condensation temperature many of the atoms enter the very 
lowest state. Naturally, any interactions between the atoms would modify the: 
form of the lowest wave function very conrsideably, but it seems most unlikely 
that they would remove this effect altogether. In the perfect gas, if it were possible 
to reduce the effect of gravity by a large factor, the effect of such a reduced field on 
the lowest states would still be quite noticeable. From this it seems at least. 
plausible to argue that, in actual liquid helium, the ‘z-states’ may be abnormally 
sensitive to an applied external field and one may expect a similar result for a 
hydrostatic pressure gradient, although it does not follow a priori that in a quantum. 
assembly a pressure gradient and an applied field can simply be added together,,. 
as is possible in classical hydrodynamics. 


§8. HELIUM 3 
It is natural to ask whether superfluid effects are to be expected at sufficiently 
low temperatures. It is probable that, at sufficiently low temperatures, He, can be 
represented by a Debye model. 'Tomonaga (1950) has shown that this is possible 
for any assembly of Fermi particles whose interactions obey some not very onerous 
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conditions, but we remain completely ignorant of the order of magnitude of the 
Debye temperature. (It should be possible to estimate this when data on the 
compressibility of He, become available.) In the meantime, it can certainly be 
said that the absence of superfluid effects in He, above 1° k can be understood if it 
should turn out that it has a Debye temperature significantly smaller than that of 
He, 

We have seen that, in the region of temperatures where a Debye model is a 
satisfactory approximation, the type of statistics becomes unimportant. The 
connection between the Debye © and a possible transition temperature can be 
illustrated by considering the simplest Debye model, longitudinal sound only 
without dispersion. For this we should have for the transition temperature the 
equations 
(a) according to Landau’s criterion p,/p 


es kr) {(4) leading to T,210°x; 


(6) according to the criterion (11) (@F/AN)p »=0 
T, 0-750 ~25°xK, 


if the shift in the zero of energy due to the attractive forces is ignored, see eqn. (4). 
The introduction of dispersion, of variation of c, with density and of transverse 
sound quanta or rotons would modify these numerical values very considerably but 
would not destroy the essential connection between 7, and the Debye © just 
referred to. 


$9. THE PROBABLE FORM OF THE SPECTRUM OF HE, AT 
LOW ENERGIES 


The transition temperature of He, is so much lower than these ‘order of 
magnitude’ estimates that it is quite clear that the Debye approximation is poor in 
the transition region. Landau (1947) assumed arbitrarily that the roton part of the 
spectrum could be described by E(p) =A + (p —pp)?/2u but Temperley (1951 a) has 
pointed out that such a spectrum would imply the unphysical assumption of an 
‘infinite’ velocity of propagation of transverse sound of low frequency and also 
that the assumption that this velocity is zero or small compared with c, would 
be inconsistent with the observed approximate agreement of the low temperature 
values of c, and c,/1/3 (Atkins and Osborne 1950, Atkins and Chase 1951). It 
therefore seems advisable to assume, provisionally, that the longitudinal and 
transverse velocities of sound are comparable, as in an ordinary crystal lattice. 

For their calculations of viscosity, Landau and Khalatnikov (1949) were mainly 
interested in the low energy end of the spectrum, and they supposed that Landau’s 
detailed model which agrees reasonably with the data for both c, and specific heat 
could, for this purpose, be replaced by the simpler spectrum 

ED) SOA Di 2) ee ee RE oy las (15) 
They derived gz ce x LOS Or? CIs eceOmiail al a 10 ceri lista (16): 
from Landau’s (1947) assumed values of A, and py and they then use this value of 
y to derive values of viscosity in reasonable agreement with experiment. 
It is of interest to compare this value of y with various other estimates : 


(a) Blackman’s (1935) work on the specific heat of an ideal crystal lattice suggests 
that significant departures from Debye’s law may occur down to temperatures as 
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low as 0:01. For a Debye temperature of 30°K (appropriate for liquid helium) 
the two terms in (15) would be comparable at an energy corresponding to 0-1 0 if 
yao 10" gem seen ee eee (17) 

(b) The assumption of an energy spectrum such as (15) enables both the specific 
heat Cy and the velocity c, of second sound to be calculated in closed form. We — 
can calculate the partition function, energy, entropy and specific heat directly, by 
integrating over all possible magnitudes and directions of the momentum p. For 
T <Q the upper limit of p can be put equal to infinity. The velocity c, can be 
derived from the formula, Landau (1941), 

(2a ae ; 
pn OS/OT 

Ps/p, being obtained from Landau’s formula (derived in the Appendix) for the 
effective mass-of an assembly in motion 


Px yb exp (c,/RT)/{exp (e|RT) — 1} 


p Nm, “RT 
where p, is the x-component of the momentum of rth state. We find easily 
Cy  967k(RT)? c(6) 2 
et Mm i+; (4) 150y cores (a YA ees (18) 
Pemet ves Re) kay 
c= + [1 7 36007 (= isin: Pecivee al) eas ae (19) 


the second terms in the brackets representing the corrections due to dispersion of 
first sound, resulting from the assumed spectrum (15). Equations (18) and (19) 
give us two further estimates of y as follows. 

Hull et al. (1951) measured the specific heat down to 0-6°K and represented 
their data by means of the formula C = 0-024 T®? cal/g deg. Comparing this with 
the expected result for a Debye model with c, ~2-4 x 10*cm/sec they conjectured 
that a change-over to the Debye law should occur at a little below 0-6° kK. Kramers 
(1951) has recently reported measurements at still lower temperatures, confirming 
this conjecture. If the two terms in (18) are comparable at 7 =0-6°, this would 
ECan ENE yes X10 e-* cmesect: (20) 

From the data of Atkins and Osborne (1951) we see that c, has fallen 
appreciably below its limiting low-temperature value at 7 ~0-4°, and, if we choose 
y of the same order as that given by (20), we see that it also makes the two terms in 
(19) comparable at just about this value of T=0-4°. Obviously a more detailed 
spectrum than (15) and the retention of more powers of T in (19) are required to 
describe the behaviour of c, at higher temperatures when the two terms in 
expression (19) nearly cancel. 

We conclude, therefore, that a value for the constant y certainly no larger than 
that suggested by Blackman’s (1935) work on crystal lattices will describe both the 
observed low temperature behaviour of both c, and Cy. It would, no doubt, be 
possible, by introducing further terms into (15), and admitting also the existence 
of a similar spectrum of transverse sound quanta, to obtain complete agreement 
with all the data, but, as we have seen, we should have a great many adjustable 
constants, and a Debye model is probably a poor approximation near the transition 
point in any case. 


Modes of Liquid Helium II 507 


‘This value of y is distinctly greater than that given by Landau and Khalatnikov 
(1949), which, if their mechanism of viscosity is the dominant one, would mean 
that we cannot describe the low-temperature behaviour of all three quantities, 
Cy, ¢ and y in terms of a single parameter y. However, it may well be that 
processes of Peierls’ (1929) ‘Umklapp’ type (three-phonon process), would give a 
‘sufficient contribution, and the variation of 7 with temperature would be similar to 
that found by Landau and Khalatnikov (1949). The numerical constant is, how- 
ever, extremely difficult to estimate and would, in any case, require knowledge of 
the whole frequency spectrum. It does, however, seem appropriate to consider 
three-phonon processes as a possible mechanism before we invoke four-phonon 
processes. 

It is sometimes suggested that the large attentuation of first sound at low 
temperatures (Atkins and Chase 1951) is evidence against the validity of the phonon 
model near absolute zero. However, since the frequency used (of the order of 
10’ c/s) is much less than the phonon frequency corresponding to a temperature of 
1°x (of the order of 101° c/s) it follows that what has been demonstrated is not that 
these latter frequencies correspond to non-stationary states, but that transfer of 
energy between such states is becoming difficult. On the simple Debye model such 
transfers are impossible, so that once again we meet the necessity of making the 
model realistic by introducing a parameter such as y, which means that the Debye 
modes are not quite stationary states. On this basis Khalatnikov (1950) was able 
to predict satisfactorily the observed trend of attenuation as a function of 
temperature. 

We now seem to be in a position to state with fair confidence that the properties 
of helium II can be explained if it can be shown that the quantum-mechanical 
many-body problem has two consequences : (i) We have already arrived at the 
idea that helium II can be described by a Debye model in statistical equilibrium 
with some states of nearly zeroenergy. Sucha picture is fully in keeping with the 
abnormal weakness of the atomic interactions and the large lattice spacing. (11) 
We must still ask whether this picture leads to superfluidity and to the apparent 
presence of two distinct densities. 

It is clear that, if a Debye model is a valid approximation at all, we must admit 
to the existence within a given volume of two quite distinct momenta, even in an 
ordinary liquid. 

1. The mean momentum P, or pz, associated with the Debye waves. ‘This 
momentum must be assigned physical reality, because, by absorption of 
the phonons, it can be transferred to material particles, as is pointed out in the 
Appendix. It is there shown further that, by considering the Debye waves to be 
in statistical equilibrium witha moving body, we can arrive at separate definitions of 
p, and 4, the effective density and stream velocity of the Debye waves, equivalent to 
Landau’s definitions (1941). 

2. The total momentum of the centre of gravity of the Debye lattice, described 
for instance by the coordinate Cy in eqn. (4), together with the momentum, if any, 
of the atoms in states of nearly zero energy. ‘This total momentum we call pV, 
where p is the total density of the assembly. 

Comparison with the two-fluid picture shows that p, must be identified with the 
density p, of the normal fluid in order to obtain the correct variation of c, with 
temperature. Further, we must identify v, with v,, because of the experimental 
fact that the entropy is carried by the normal fluid and Dingle’s result (1951) that 
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the entropy of a Debye assembly is transported with the velocity 7). Since sound 
waves never transport any mass, we must identify the rate of transport of mass, 
PnUn +Psv, With pV so that we have v,=(pV —py%,)/(p—py)- On this basis, we 
can derive equations in first approximation equivalent to those of the two-fluid. 
theory, using the postulates of conservation of mass, of rate of change of material 
momentum being proportional to gradient of hydrostatic pressure, of conservation 
of entropy (transported with velocity v,) and finally of rate of change of momentum 
of Debye waves being proportional to gradient of sound pressure. Sound 
pressure is defined in a manner analogous to radiation pressure (Ward and Wilks. 
1951). 

Superfluid effects manifest themselves, not through wv, being different from. 
zero but through », being different from v,, that is, through v, being different from. 
V. Now on the simple Debye model, there is no coupling between any of the 
C, and C), but we have seen that we have to admit such coupling between the C,’s. 
in an actual liquid in order to explain the finite viscosity and thermal conductivity. 
Furthermore the absence of superfluid effects in ordinary liquids indicates a 
reasonably strong coupling between C, and the remaining C,’s. In the Appendix, 
it is shown that Debye waves will be in statistical equilibrium with a body only ifits. 
velocity of translation is equal to vu, and we may apply this result, in particular, 
_ to the motion of the centre of gravity of the lattice. We must suppose that in an — 
ordinary liquid the coupling (due to the cubic and higher terms in the interaction 
energy) is strong enough to maintain statistical equilibrium between C, and the 
other C,’s. In liquid helium II we can explain a difference between v, and V in one 
of two ways: (a) the presence of atoms in z-states means that V may no longer be 
quite equal to the velocity of the centre of gravity of the lattice, which velocity is. 
still constrained by the coupling to be equal to uv). (6) The coupling of C, to the 
other C,’s suddenly becomes too weak to maintain statistical equilibrium. 

On present knowledge we cannot exclude either of these possibilities, but (a) 
seems far more plausible than (d) in view of the sudden onset of superfluidity and 
the known existence of a thick film. 


CONCLUSIONS 

We conclude that the Debye model is a satisfactory approximation to He, at 
sufficiently low temperatures, and that its range of application can probably be 
extended somewhat by taking account of the probable presence of transverse waves, 
and the variations of the effective propagation velocities with density, temperature 
and frequency. It is not yet possible to say whether it will ultimately prove to be 
a useful approximation in the region of the transition temperature. 

The Bose-Einstein gas-like model has been very valuable in gaining physical 
insight into the behaviour of liquid helium II and it seems very likely that two of its 
features, the condensation phenomenon, and the sensitivity of the lowest states to 
an applied field of force have analogues in the many-body model, even when the 
interactions are very strong, and that both these analogues will find a place in the 
theory of liquid helium IT. 

Many of the properties of helium II follow from a model in which a Debye 
lattice is in statistical equilibrium with atoms in states of small energy and entropy. 
This model would have a transition temperature above which the latter states 
would be unoccupied. 
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APPENDIX 


Revised Calculation of the Effective Mass of a Moving Assembly 


We re-apply the methods of statistical mechanics to an assembly in which we 
specify linear momentum as well as energy. (The extension to angular 
momentum is immediate.) We suppose that each level in the assembly may be 
occupied any number of times, and that to the th level we may assign an energy e, 
and an x-momentum p,, the latter quantity being the change in momentum 
produced at a wall if the occupation number of the 7th level changes by unity as the 
result of a collision. (Since the collision of excitations such as those described by 
eqns. (4) and (5) with a wall does not yet seem to have been considered explicitly we 
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can, strictly speaking, only infer the relation «=c|p| for phonons from the 
correspondence principle. We shall however assume it generally. With this 
definition of the momentum of phonons, momentum is conserved when electrons 
collide with a crystal lattice. * 

The number of ways in which one can realize a state of total energy H and total 
x-momentum P is the coefficient of y’z” in the generating function II,(1—y?rz*r)-1. 
If the total number of excitations were fixed (e.g. if the levels referred to individual 
particles) we could introduce a third selector variable. Then, applying the saddle- 
point method, we have for the entropy S 


S~kX, log (1—¢”6%)4—-kPlogd—RElog@, _....... (A 1) 
where y =¢, z=8 is the position of the saddle-point. By general theory we may 
infer that 9 =exp(—1/RT). For the free energy we have 


F=kT%, log{1—¢”exp(—<,/RT)}+RTPlog¢d, ....... (A2) 
where ¢ is determined by 
P, as 
Pas je ea RS | (A 3) 


~ 1 —gPexp(—e/RT) 


We now assume that, to every level with a given p, corresponds a level with the same 
energy but an x-momentum —p,. This follows quite generally from symmetry 
considerations. Now put¢=e%. We then have 


P, exp («,/kT) sink (Bp,) er 


om - 1+exp(2e,/RT)—2exp(e,/RT)cosh(Bp,)’ 9-7” 


kT 
F= — dlog{l+exp(2e,/RT) —2exp(<,/RT)cosh(Bp,)} +RTP,B. 


2 
For small P, f is also small, and these expressions reduce to 
2 


, {exp (e,/RT)— 1}" 
RT .. exp(e,/RT) 
F=RT Xlog {1 - —e,/RT — Bp? & ———* 
We now determine the value of f so that this assembly shall be in statistical equili- 
brium with a body of mass M, moving with a velocity v, in the x-direction. The 
total momentum of assembly and body is P+ Mv which is constant, while the 


total free energy is 


+0 (6%). 


TB. _ exp(c,/RT) 


kT d log {1 —exp(—e,/RT)}+ s 5 +$.My,. 


2 + {exp(e,/RT)—1} 


B has to be chosen so that F is a minimum, subject to the momentum being 
conserved, which requires that 8B =v,/RT. The mean momentum of the assembly 
1 y_PPexp(e/RT) 
kT , {exp(e,/RT)—1} 
which gives an ‘effective mass’ 

1, pPexp(e/kT) 

kT, {exp(e,/RT) —1}? 


* See, for example, Wilson (1936), chap. 6. 


is then 
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which is equivalent to Landau’s expression —%,p,2(0/de,){exp(e,/RT) — 1}-! and 
we can complete the calculation for the phonon model by putting «=c|p|, 
P,= |p| cos é@ and integrating over all angles and all positive values of |p|. 

This treatment still assumes that a phonon has momentum that is transferable 
to material particles, but avoids Landau’s dubious assumption that the energy of 
such a phonon is different for different observers. The treatment on p. 495 indeed 
seems to show that this assumption is incorrect. 

It is of interest to recall an old result of Brillouin’s (1922) who concluded that, 
if an assembly consisting of Debye waves were sheared with uniform velocity, the 
sound waves would be refracted in such a way as to tend to increase the momentum 
content of the faster-moving regions of the assembly. This shows that a physical 
mechanism, by which a moving assembly can approach the state of statistical 
equilibrium calculated in (A4), actually exists. We also notice that 8, and 
therefore the mean momentum P, both depend uniquely on v,.. This means that 
an assembly of Debye waves cannot be in statistical equilibrium with two bodies, 
e.g. a wall and the centre of gravity of the crystal lattice, unless they are both 
moving with the same velocity. 


The Thermal Conductivity of Liquid Helium I 


By R. BOWERS* 
The Clarendon Laboratory, Oxford 


Communicated by K. Mendelssohn; MS. received 14th September 1951 


ABSTRACT. The thermal conductivity of liquid helium has been measured between 
2-2 and 4° k. It is found to increase with temperature over the whole range and shows 
no sign of anomalous behaviour just above the A-point. 

This variation is compared with that of the viscosity and entropy in the same 
temperature range. 


St INTRO DUG ELON 

N an earlier paper (Bowers and Mendelssohn 1950) we discussed the fact 

that the fall in entropy associated with the lambda transition begins at a 

slightly higher temperature than the actual transition point of 2-19° kK. 
It was thought that this preliminary entropy drop might manifest itself in the 
transport properties of liquid helium I within about } degree of the A-point. 
In the earlier work, a fall in the viscosity in this temperature region was indeed 
observed. The present paper describes the measurement of the thermal 
conductivity of liquid helium I.t 

As was the case with the viscosity determinations, earlier work below 2:19°k 
(e.g. Keesom, Saris and Mayer 1940, Allen and Ganz 1939) gave no indication 
that high heat flows would persist above the A-point. However, the type of 
heat flow observed by these authors in helium II was of a very anomalous kind 
which cannot be characterized by a thermal conductivity in the normal sense. 


* Now at the Institute for the Study of Metals, Chicago. 

+ A short preliminary note of the results of this work has been published earlier (Bowers and 
Mendelssohn 1951). More recently, a determination of the heat conductivity of He I has been 
reported by Grenier (1951), whose results are in satisfactory agreement with ours. 
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Our work on the viscosity suggested that even though the anomalous type of 
transport vanished at the d-point, there still remains the possibility that the 
normal heat conductivity could change markedly within } degree of the -point. 
The only previous determination of the heat conductivity of liquid helium I 
was that of Keesom and Keesom (1936). They measured a single point at 3-3°K 
and obtained the value 6 x 10-5 cal/deg. cm sec, which is of the order of the thermal 
conductivity of a gas at room temperature. However, Pellam and Squire (1947) 
have suggested that there might be a considerable variation of the thermal 
conductivity of helium I with temperature, the conductivity increasing with 
decreasing temperature. This suggestion arose from their study of ultrasonic 
attenuation in liquid helium. 


§2. METHOD 


The measurement was carried out on a vertical cylinder of the liquid. Heat 
was supplied to the top of the cylinder and was conducted down to a heat sink 
at the bottom of the cylinder which could be maintained at a constant 
temperature. ‘Two thermometers were fixed at different heights within the 
cylinder, mounted parallel to its base. For zero side loss of heat from the cylinder, 
the thermal conductivity K may be calculated from QO=KA(6,—6,)/L where 
Q=?R=the rate of production of heat in the heating plate, 6, and 6, are the 
equilibrium readings of the upper and lower thermometers for the heat flow Q, 
A is the area of cross section of the cylinder, and L is the distance between the 
thermometers. 

The actual form of the conductivity cell is shown in fig. 1 (a). A is a copper 
base on which were mounted the two resistance thermometers C and D. The 
heater B and the resistance thermometers were mounted parallel to the top of 
the copper base, being supported by three tubes of german silver of } mm diameter 
and §,mm wall thickness. The heater was a coil of Eureka wire of about 110 ohms 
and the heat generated in it was distributed uniformly over the circular cross 
section by the plate of copper to which it was attached. The thermometers, 
which were required to measure the temperature in a plane of circular cross 
section, consisted of leaded-brass wire wound over a former of amber. About 
20 cm of resistance wire was required in order to obtain the necessary sensitivity 
and this wire was wound on to the former in the manner shown in fig. 1 (8). 
It was unavoidable that the wire, which ideally should have been entirely 
accommodated in one plane, tended to bow on the former. However, since the 
bowing was symmetrical about the plane of the former, no appreciable error was 
introduced and the recorded temperature could be taken as equal to that in the 
annular plane. The leads for the thermometers and for the heater (three leads 
for each thermometer and two for the heater) passed out of the cell through 
tube E and then through eight ‘Covar’ glass-to-metal seals soldered into the 
copper base. ‘The heater and thermometers were enclosed by an inverted glass 
Dewar vessel G which ensured that all the heat generated in the heater flowed 
down to the heat sink, i.e. it prevented side loss. 

The conductivity cell was suspended in the experimental space of a small 
scale Linde liquefier (Daunt and Mendelssohn 1948). During an experiment, 
the apparatus was surrounded by liquid helium, and liquid helium could be 
condensed into the measuring cell through tubes E and H. 

Measurements of the thermal conductivity of liquid helium I present a 
number of peculiar difficulties which result from a combination of high specific 


| 
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heat, very low heat conductivity and an extremely high fluid mobility. Since 
a detailed discussion of these difficulties has been presented elsewhere (Bowers 
1951), the experimental problems involved need only be mentioned briefly in 
this paper. It was found that even small vibrations in the laboratory were 
sufficient to produce convection currents inside the cell which made accurate 
determinations of equilibrium temperatures difficult. ‘The experiments had 
therefore to be carried out at night, but even under these conditions, vibration 
represented a serious cause of random experimental error. A further difficulty 
was introduced by the fact that the heat conductivity of any suitable container 
material, even glass, is very much larger than that of liquid helium I (ratio of 
3-5 to 1 in the case of ordinary soft glass at 3°k). It is therefore essential to use 
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—— Copper 


Wood's 
Alloy Seal 
ee 
Gy, 
Fig. 1 (a). The conductivity cell. Fig. 1 (6). Construction of thermometer. 


a cylinder of helium of the maximum radius possible and the thinnest container 
wall in order to prevent an excessive quantity of heat flowing down the container. 
A limitation on the size of apparatus which can be placed in our liquefier-cryostat 
led to the fact that about 50%, of the heat from the heating plate flowed down the 
inner glass wall of the container Dewar. ‘Thus a correction of this magnitude 
had to be applied to the measured conductivity of the cell in order to obtain the 
conductivity of the liquid helium itself. 

The experimental procedure adopted for the measurements was as follows. 
As helium was being liquefied into the cryostat containing the conductivity cell, 
a pressure of about 100cm Hg was maintained in the liquid in the conductivity 
cell. During the filling of the cryostat, the cell was found to fill to about two-thirds 


‘of its total capacity with liquid. The cell was completely filled by pumping the 


outer bath below the A-point, when the condensation became very rapid and all 
temperature gradients in the liquid disappeared. ‘The liquid in the cell was 
then brought to the temperature required for a determination of the conductivity 
or for a calibration of the thermometers by warming the outer bath to this 
temperature. The temperature of the liquid in the cell followed that of the 
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outer bath quite rapidly because the principal process of heat flow involved was. 
that of convection. Due to a combination of low viscosity and large expansion 
coefficient in liquid helium I, thermal convection currents can equalize 
temperatures rapidly. 

' As the liquid in the cell warmed up from below the A-point, the pressure 
exerted on it was maintained at a steady value of 90cmHg. ‘This pressure, 
which always exceeded the vapour pressure, was maintained throughout the 
experiment so as to prevent any evaporation of liquid inside the cell. Moreover, 
during a determination of the conductivity, the cell was regularly inspected for 
the appearance of bubbles. 

The resistance of the thermometers was calibrated for a number of different 
temperatures determined from the vapour pressure of the liquid in the outer 
bath. The sensitivity of the thermometers changed from 3Q/deg.k at 4°K to. 
4Q/deg. K at 2:2°x. Individual calibration readings fell on the smoothed 
calibration curve to within about ;4,°K. This evidently represented the accuracy 
with which we can determine the absolute temperature of the liquid helium — 
within the cell. Although the main shape of the calibration curves was. 
determined in these initial experiments, calibration points were subsequently 
taken before each determination of the thermal conductivity but no change in 
the calibrations was found. 

For making a measurement of the heat conductivity, the following procedure 
was adopted. Having filled the cell below the A-point the outer bath was brought 
to a temperature slightly below that at which the heat conductivity was required. 
When conditions had become steady, a current was passed through the heating 
plate of the cell. ‘The heat required can be divided into two parts. Firstly, 
enough heat is required to establish a convenient temperature gradient 
(0-1 or 0-2°k/cm) along the cell and secondly, a constant heat current is required 
to maintain this gradient during the measurement. The first quantity of heat is. 
determined by the thermal capacity of the liquid in the cell. The second quantity 
of heat is to be supplied asa steady heat current and is determined by the thermal 
conductivity of the liquid in the cell. With our arrangement, it was not convenient 
to supply the heat necessary for establishing the temperature gradient using the 
same steady current required to maintain it. The thermal capacity of the liquid 
is so large and its heat conductivity so small that the steady flow heat current 
would have taken four or five hours to supply the total heat necessary for 
equilibrium. It was clearly more convenient to put in the heat required for the 
establishment of the temperature gradient in a short time, using a large current 
through the heater for a few minutes and then reducing the current to the value 
required to maintain the gradient of temperature. Using this technique, 
equilibrium could be reached in about 2 hours and readings were continued for 
a further two hours in order to establish accurately the temperature difference 
between the thermometers. 

$3. RESULTS 

A series of observations was made at temperatures between 2:2°x and 4°K. 
The results of these measurements are shown in fig. 2. In this figure we have 
plotted an effective conductivity K, of the cell filled with helium between the two 
thermometers. ‘The values of K, together with the temperature differences used 
are given in the table. K, is not the conductivity of the helium because it 
contains a large contribution due to conduction down the inner wall of the Dewar 
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and a small contribution due to conduction down the german silver supports 
_ and the Eureka wires. The correction necessary to allow for this residual 
conductivity is also shown in fig. 2. This correction was calculated from 
measurements of the conductivity of a specimen of glass of the same batch as 
that of the Dewar vessel; it has been found that the conductivity does not vary 


Details of Heat Conductivity Measurements 


rag T Is K 
(ma) (° k) (° k) (x 10-5 cal/deg. cm sec) 
0-50 0-047 2-780 8:20 4-25 
0-70 0-104 2-824 7-44 3-42 
0-75 0-111 2-719 7:98 4-13 
1-01 0-183 2-786 8:79 4-82 
0-90 0-166 27573 7:69 4-09 
0-80 0-171 2-425 5:91 2-54 
diel'5 0-216 Soilsles 9-65 5215 
0-65 0-105 2-292 6:38 3-23 
0-88 0-161 2°353 7-63 4-38 
0-91 0-202 2-505 6-45 3-00 
0-75 0-120 2°320 7-40 4-20 
0-86 0-177 22352 6:59 3-34 
0-70 0-131 2:450 5-90 2°50 
0-70 0-121 2-281 6°39 3-26 
0-55 0-076 2251 6:05 3-00 
1-25 0-221 3-655 TO 5-65 
Ta 0-165 3-941 ey, 6:90 
1-05 0-167 3-097 10-4 5:95 
O75 0-121 2-450 (35 3:91 
0-97 0-189 2°515 . 8-30 4-80 
0-87 0-184 2-444 6°52 Su 2 
hat? 0-216 3-380 9-15 4-22 
1-24 0-196 3-973 12-4 6°57 


i=steady current through heater (resistance =109Q), AT=resulting temperature difference, 
T=mean temperature of determination. 
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Fig. 2. The variation of effective conductivity with temperature. 
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appreciably from one sample to another.* ‘The evaluation of the residual heat 
flow by a direct method involving evacuation of the measuring cell was not 
possible because the heater and thermometers were not connected directly to 
the glass. The use of a gas at low pressures was also not practicable because of 
the long equilibrium times involved in this measurement. However, it should 
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Fig. 3. The variation of heta conductivity of helium I with temperature. 


be remembered that our main interest in this work was the variation of the 
conductivity of helium I with temperature, particularly just above the A-point, 
and not its exact absolute magnitude. Thus a precise evaluation of this correction 
is not of the greatest importance. 

The final curve of the heat conductivity for liquid helium I is shown in fig. 3. 


§4. DISCUSSION OF RESULTS 

The scatter of points in fig. 2 is rather large; the mean scatter of points from 
the straight line drawn is6%. ‘There does not appear to be any correlation between 
the deviation of any point from the line and the temperature gradient used in 
its measurement. We have no evidence that would suggest that the heat flow 
was non-classical, so that we regard the scatter of points as purely experimental 
error. Probably the most important single contribution to this was the 
disturbance due to vibration mentioned earlier. During the determinations, 
fluctuations in temperature were observed on the cell thermometers which were 
of the order of 10-20% of the temperature difference between them. However, 
the random error introduced into the final result was reduced by averaging the 
observed fluctuations. ‘These fluctuations observed in experiments carried out 
at night were an improvement by at least a factor 2 over the disturbance observed 
during the day. It was this fact which led us to investigate mechanical vibration 
as the source of the disturbance (Bowers 1951). 

Our observations suggest that the heat conductivity of liquid helium I varies 
with temperature, falling from 7x 10-°cal/deg.cmsec at 4°K to about 
3x10 cal/deg.cmsec at 2:2°xK. The relation between conductivity and 
temperature is roughly linear. Extrapolating the straight line drawn through the 
points of fig. 3 to zero conductivity, it is found to cut the temperature axis at 0-75° k 
The absolute value of our results is in good agreement with the single value ane 
to Keesom and Keesom. The results of our experiment are not precise enough 
to allow a detailed discussion of the variation of the thermal conductivity of 


* We wish to thank Dr. R. Berman for determining the conductivity of the glass required for 
evaluating this correction. 
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liquid helium I with temperature. It should be recalled, however, that we were 
primarily interested in the question of whether the conductivity changes in a 
way which would suggest a preparation for the superfluid state as the A-point 
was approached. Our results do not show any marked anomaly near the 
A-point, although the scatter is too large to rule out changes in the conductivity 
of the order of 10%. In no case did we observe very high values of the 
conductivity just above the A-point, so that the effective conductivity of the 
liquid helium must change very sharply on passing through the )-point. 

In order to get some further qualitative information on this sudden change 
at the A-point, a short experiment was carried out which involved heating the 
cell through the A-point. During this experiment the temperature difference 
between the two thermometers was measured as the liquid helium was slowly 
heated from below the A-point. It was observed that a temperature difference 
between the thermometers began to appear within a few thousandths of a degree 
above the A-point. This is further evidence for our view that the change in 
effective conductivity is very rapid at this point. 


§5. CONCLUSION 

In order to analyse the results of this experiment we shall consider to what 
degree they are consistent with the earlier observations on the viscosity of 
liquid helium I and with the known shape of the entropy curve (Keesom and 
Keesom 1935). The viscosity showed a preliminary fall just above the A-point 
in the same temperature region in which the entropy showed its first anomalous. 
decrease. So that both the viscosity and the entropy just above the A-point 
indicate some preparation for the superfluid state. Our earlier work (Bowers. 
and Mendelssohn 1949, 1950) showed that the observed fall in the viscosity really 
represented a fall in the true viscosity of helium I, there being no sign of a flow 
with the anomalous character as is observed below the A-point. ‘Thus the 
preparation for the superfluid state took the form of a change in the well-defined 
viscosity and could not be explained by the inception of a small degree of 
superfluidity. The heat conductivity, however, does not appear to show a 
similar preparation for the superfluid state, but such an observation is not 
necessarily inconsistent with the variation of viscosity and entropy in this. 
temperature region. 

According to the classical kinetic theory of gases, the thermal conductivity K, 
the viscosity 7 and the specific heat at constant volume Cy, are related by the 
relation K=const. 7C,. Although helium I is a liquid, it does have a number 
of gas-like properties. Due to its large zero-point energy, its density is quite 
low (Simon 1934). Moreover, its viscosity increases as the temperature increases, 
a property characteristic of gases. 

The dependence of K on the product 7C, may provide an explanation for 
the insensitivity of the thermal conductivity to the approach of the superfluid 
state. As the temperature is reduced below 2-6°k, the viscosity decreases while 
the specific heat increases, so that the product 7C, may not change appreciably 
in this temperature region. A more quantitative analysis of the product 7C, is 
hardly justified, considering the crude gas model underlying it. It is, however, 
worth mentioning that if the mean free path A is calculated from the gas relation 
K=4pcaC, then a value of the order of 10~-7cm is obtained. This is of the order 
of the interatomic distance in the liquid as calculated by London (1936). 
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Heat Transport in Lead-Bismuth Alloys 
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Clarendon Laboratory, Oxford 


Communicated by F. E. Simon; MS. received 3rd March 1952 


ABSTRACT. ‘The variation with temperature and with magnetic fieid of the heat 
conductivity of a series of lead—bismuth alloys has been measured. As the concentration 
of bismuth is increased the heat conductivity in the normal state falls below that in the 
superconducting state, contrary to the case in pure superconductors. Estimates have been 
made from the experimental data of the processes contributing to the heat conductivity, 


and of the resistances due to the various scattering mechanisms limiting these 
conductivities. 


Sel uN OD Ui Cr On 
RELIMINARY results of a series of measurements of the thermal 
conductivities of superconductors have already been published (Mendelssohn 
and Olsen 1950a,b,c). These results showed that processes other 
than the conduction by the electrons must be important components of the 
total thermal conductivity in the superconducting state. It is the object of this 


paper to present fuller experimental data and to attempt a quantitative analysis 
of the observed effects. 


§2. APPARATUS 

The measuring part of the apparatus is shown in fig. 1. The specimen S 
was placed horizontally and connected thermally to the expansion chamber C of 
an expansion liquefier (Simon 1932) by the copper connecting rod K. The 
temperature gradient produced by heating electrically with the heater W was 
measured by two gas thermometers T’ and T”. The pressures in these 
thermometers were transmitted to an external manometer system by the 
capillaries M. ‘The specimen and thermometers were surrounded by the 
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radiation shield R, whose interior communicated with the vacuuth space of the 
liquetier. Rees fag 1 

The temperature of the expansion chamber could be kept at any temperature 
between 20°K and about 10°xK by allowing hydrogen in a container (not shown 
in the diagram) attached to the expansion chamber to evaporate at an appropriate 
pressure. ‘Temperatures between this and the critical point of helium could 
be obtained by allowing the helium in the expansion chamber to expand at a 
suitable rate (MacDonald and Mendelssohn 1950). Lower temperatures were 
obtained by letting the helium boil under a reduced pressure. 

The thermometers were filled to a pressure of between 15 and 75cm Hg 
when the apparatus was at 20°K, and subsequent pressure differences between 
the two thermometers were measured by a differential oil manometer. Since 
the difference in height of the two oil levels could be read by eye to about 0-2 mm, 
changes in temperature differences of less than 0-001° could be read with the higher 
filling pressures in the thermometers. The temperature differences set up by 
the heater were usually made of the order of 0-1°, and there was consequently 
a possibility of obtaining the thermal conductivities with an accuracy within 1%. 


Fig. 1. Arrangement of specimen and thermometers. 


The conduction along the thermometer capillaries and supports was too 
small to be measured in the absence of a specimen, and has therefore been 
calculated. It is found to lie below 1% of the conduction of the specimen in all 
cases. Departure of the helium in the gas thermometers from the ideal gas laws 
causes a negligible error in the temperature difference, but a correction has 
been made to the absolute temperatures for this effect. ‘The dead volume of the 
capillaries introduces an error which is always difficult to correct for, but in our 
case this is less than 1%. One of the most dangerous errors in thermal 
conductivity measurements, because it may pass unnoticed, is that due to 
conduction through the residual gas in the vacuum surrounding the specimen. 
The use of a combined liquefier—cryostat has the advantage that this danger is 
completely eliminated, since a deterioration of the vacuum is immediately shown 
by an abnormal evaporation of helium from the cryostat. ‘This evaporation 
would become obvious when the conductivity of the residual gas reached about 
10-* watt/degcm which would correspond to an error of 5 x 10-4 watt/degcem 
in the measured conductivity of a specimen of usual size. 
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The specimens were made from pure lead obtained from Johnson, Matthey 
(Lab. No. 1932) to some of which 0-02, 0-1, 0-2, 0-5, 10 and 30% of bismuth (Lab. 
No. 4213) from the same source had been added. The shape of the specimens 
is seen in fig. 1. The two copper discs D, to which the thermometers were 
eventually attached with solder, were a part of the mould. They were wiped 
inside with lead to make them adhere to the specimen as it solidified. This 
avoided the falsification of magnetic field values by persistent currents set up 
in rings of solder. 


§3. VARIATION WITH TEMPERATURE 
(i) Experimental 


In fig. 2 the thermal conductivities of two specimens of pure lead are shown. 
These specimens were cast in vacuum direct from the ingots obtained from the 
makers. They may be expected to be polycrystalline and less pure than the triply 
recrystallized specimens of de Haas and Rademakers (1940, Rademakers 1949). 
This is also in agreement with the fact that we find a maximum in the conductivity _ 
which is three times lower than theirs. While their specimens were stated to 
be single crystals, both our specimens probably had a fairly small grain size.* 
The superconducting conductivity K, of our specimens departs sharply from 
the normal conductivity K, below the transition temperature 7, just as has 
been found by previous workers. However, while they found that it falls to a 
minimum and then rises and has a maximum at about 3:5°K, below which 
temperature K, falls steeply, no such maximum was observed in our specimens. 

When impurity is added (see fig. 3) the absolute value of the normal 
conductivity is decreased, and the maximum in the conductivity is forced to 
higher temperatures, as predicted by theory for the normal state (Makinson 1938, 
Sondheimer 1950). It is also seen that as the impurity increases the low 
temperature end of the superconducting curve is raised in relation to the normal 
curve. When we come to the specimen with 0-2% Bi the low temperature end 
of the superconducting curve actually lies above the normal curve, which it 
crosses at 4:8°K. At this impurity content the high temperature end of the 
curve still remains below K,, but when the impurity is increased to 0-5°% and 
10%, K, becomes greater than K, throughout the range of measurement. The 
thermal conductivity of an alloy containing 30% of bismuth is shown in fig. 4 
on an enlarged scale. 

It is clear that these alloys behave quite differently from our pure specimens, 
and from all published results for pure metals. On the other hand the behaviour 
described here is not a general feature of high impurity content alloys. Thus an 
alloy of lead with 30% tin has a thermal conductivity in the superconducting 
state which is lower than that in the normal state throughout the range of 
measurement (Mendelssohn and Olsen 1950 a). The results shown in the same 
publication for a pure sample of niobium should not, however, be taken to 
indicate a real deviation from the usual characteristics of a pure superconductor. 
The magnetic field available during the measurements was not large enough 
to destroy superconductivity completely at the lower temperatures, and the 


* The grain size was subsequently estimated as = mm. by etching a surface cut in the lead with 


a chisel. ‘This seems to be surprisingly small, and although this method of preparing a surface is 


generally recommended, the result obtained thus must probably be accepted with a good deal of 
reserve. 
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curve shown as normal does not in fact represent the thermal conductivity with 
superconductivity destroyed. The crossing of the two curves in this case is 
therefore of little significance. 

(ii) Discussion 

The thermal conductivity of a normal metal may be split up into two 

components, and we may conveniently write (Kénigsberger 1907) 

I OW Ct eer re (1) 
Here K is the total thermal conductivity, K, the electronic thermal conductivity, 
and K, is the thermal conductivity due to the lattice vibrations. We shall denote 
by K,, and K,, the electronic and by K,, and K,, the lattice thermal 
conductivities in the normal and superconducting states respectively. 
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When we come to examine K, and K, separately it is simpler to consider 
their reciprocals. This is because the total thermal conductivity is made up of 
two conductivities in parallel, each of which is limited by resistances acting 
as if in series. In the normal state we have for the electronic contribution 


(Makinson 1938) 1/K.. <AT?+B/T 2M 
Km Ie oe On ner ee. 
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where the first term AT? is due to the scattering of the electrons by lattice -_ 


vibrations with 4 a characteristic of the metal, and the second term is due to 
impurity scattering with B proportional to the residual electrical resistance. 

For the lattice resistivity 1/K,, which will be denoted by W, the main 
contributions are: W, due to the scattering by the electrons, W, due to 
scattering at the boundaries of the specimen and crystal grains, and Wp due to 
scattering by impurities. We thus have 

| 1/K,=W=Wyt+WytWp. nena (3) 

This simple sum may be used as a first approximation, but it cannot be relied 
upon implicitly. Klemens (1951) has pointed out that in dielectrics the 
corresponding expression becomes unreliable when two or more of the 
contributions to W are of the same order of magnitude. ‘The experiments of 
Hulm (1950) and Berman (1951a) on normal alloys confirm that for these 
W, is proportional to T~? as is to be expected theoretically. The expected 
linear dependence of Wy upon temperature, and the proportionality of W, to T* 
have been verified in experiments on dielectric crystals (e.g. Berman 1951 b). 
Berman has also shown clearly the inaccuracy of the additive resistance 
approximation in the region predicted by Klemens. 

When a substance becomes superconducting the electronic thermal 
conductivity will change from its normal value, and it is not impossible that the 
lattice thermal conductivity may do the same. It is therefore of interest to - 
examine the behaviour of both the ratios K,,/K,, and K,,/Kgp. 

The first of these two has been found for a number of metals by Hulm (1950). 
In a pure metal K,, is negligible compared with K,,, and if it is assumed that this 
is also the case with K,,, then K,/K, will provide a good estimate of K,,./Kgp. 
Hulm finds that in substances where the impurity scattering term is the 
dominant one in the normal electronic resistance, the ratio K,,/K., is fairly 
accurately given by the function suggested by Heisenberg (1948) 

Race) i: 
Koa ETT) Oe Se as 
We shall follow Hulm and denote this function by f(7/T,). At the lowest 
temperatures this simple formula does however need some modification 
(Goodman 1951). 

For substances where scattering by the lattice vibrations is dominant, Hulm 
has suggested the empirical relation 


Kog/Kon = (T/T)? =(T/T.) say. wees (5) 
As the temperature is lowered the behaviour of the ratio K,,/K,, may be expected 
finally to tend to that of eqn. (4) since the lattice vibrations will die out as T°. 
In many cases a metal will not fall definitely into either group, and we shall 
need an interpolation formula to describe these intermediate cases. A suggestion 
for such an interpolation formula for K,, has been made by Hulm (1950), but 
it does not appear to give satisfactory agreement with the published results. 
We find that the experimental results are described slightly better by the formula 
Ku _ a T/T Se( T/T.) +f(LIT,) ; 
en T/T.) +1 ae): sSavinoeGies (6) 
Where « is defined by «= AT,*/B so that it is the ratio of the two contributions 
to 1/K,,, in (2) at the transition temperature. The formula (6) thus represents 
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a mean of the two functions f(7/7,) and g(7T/T.), weighted according to the 
amounts of the two scattering processes present at a particular temperature. 

An examination of eqn. (2) shows that a plot of 7/K,, against T° will be 
a straight line, and suggests an obvious way in which « may be evaluated. This 
has been done for some published experimental results, and in fig. 5 observed 
values for K,/K, for specimens with «=0-1, 1 and 10 are compared with values 
of K,,/K,, calculated from eqn. (6). Agreement is seen to be fair in most cases 
at the higher temperatures. Hulm has shown how well the pure Heisenberg 
formula, to which (6) reduces for «=0, describes the experimental results for 
metal with high Debye @ where « is quite small even in very pure specimens. 
For the intermediate cases our calculation does at least give a better estimate 
of K,,/K,, than either f(T/T,) or g(T/T,) taken separately. At the lowest 
temperatures there is always a good deal of disagreement. 


@ Tin «=I (Rademakers 1949) 
w Mercury @=l 
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Fig. 5. Kes/Ken calculated from eqn. (6) compared with experimental values for Ks/Kn. 


If the theoretical curves in fig. 5 are compared with our results in impure 
substances in fig. 3, it is seen that there is a serious disagreement and that the 
observed K, is larger than predicted. It must however be remembered that the 
curves in fig. 5 only show the behaviour of K,,/K,,, and that the conductivity 
has a lattice component and may also contain other contributions which we 
will denote by X. We must therefore write 

Kit Kec es Ae yh Bee ei ers (7) 
It is readily seen that if K,, does really have the value predicted from the 
calculations above, and if we take it that K,, is not appreciably different from its 
normal value, then X will represent an important part of the total conductivity. 

It seems clear, of course, that the need for this extra term must be due to 
the fact that some important contribution has been omitted in calculating 
either K,, or K,,. In fact formula (7) is very similar to one proposed in an 
earlier paper (Mendelssohn and Olsen 1950a) in which the extra term was 
considered to be part of the electronic conductivity in the form of a circulation 
process similar to that in liquid helium II. It was then pointed out that while 
this process could be made to fit our observed data, numerical evaluation would 
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be unconvincing because of the difficulty of checking the assumptions made for 
this process. It is, of course, also possible that the high value of K, may be 
due to some increase in the lattice conductivity in the superconducting state. 
Indeed Hulm (1950) has pointed out that a reduced scattering of lattice waves 
by the superconducting electrons may cause an enhanced thermal conductivity 
of the lattice. Since the proportion of X or K,, in such substances as our alloy 
of lead with 10° bismuth is larger by an order of magnitude than found in any 
of Hulm’s results, explanation by an enhanced lattice conductivity did not 
appear convincing at first sight. However, the results on the graded series of 
lead—bismuth alloys described in the present paper offer an opportunity of 
estimating the postulated enhanced lattice conductivity in a more quantitative 
manner. 

‘The most convenient first step in our estimate of the magnitude of X will 
be to calculate X+K,,. Since X+K,,=K,—K,,, all we need do is to 
subtract K,, from our observed superconducting thermal conductivity. We 
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have used formula (6) to derive K,, from the normal electronic thermal 
conductivity K,,. The values of « required were found as described above and 
ranged from very nearly zero for the most impure specimens, to «=2 for the 
specimen containing 0-02% bismuth. In principle K,, may be derived from the 
electrical conductivity and a knowledge of the appropriate value of the 
Wiedemann—Franz constant. Figure 6 shows how we have used this procedure 
to obtain K,, for the specimen containing 10°% bismuth. 

This simple method of calculating K,, cannot be used with the purer 
specimens, because the electrical data are not sufficiently accurate, and because 
the Wiedemann—Franz constant is not reliably known for them. In order to 
calculate K,, we have therefore subtracted a calculated value of K,, from the 
total normal thermal conductivity. This value for K,, has been derived from the 
lattice conductivity of the 10% bismuth specimen between 7°K and 21°K. In 
this region the only important factors limiting K,, are electron and impurity 
scattering, and a knowledge of their theoretical dependence on temperature 
allows them to be separated. We thus obtain Wy, and Wp, for the impure 
specimen. In the purer specimens the impurity scattering will be negligible 
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compared with the scattering by electrons, and it has been assumed that the 
resistance from this cause, W,,, is the same in the pure and the impure specimens. 
The value thus obtained for Wy, in lead is 25 deg cm/watt at 10°x. This can be 
compared with the value of 50 deg cm/watt to be expected from the expression 
given by Makinson (1938). We would here remark how much better is the fit 
of this expression for the lattice thermal resistance in the presence of electrons, 
than is the fit of the expression given in the same publication for the electronic 
thermal resistance in the presence of lattice vibrations. 

In fig. 7 is shown the conductivity of the specimen containing 0:5°% Bi. 
K,, has been calculated by subtracting from the total normal thermal conductivity 
the value for K,, described above. 
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Fig. 8. Contribution to Ks other than Kes for specimens with various percentages of bismuth. 


In fig. 8 are given the differences between K,, and K, for the whole series of 
lead bismuth specimens with the exception of the two pure specimens for which 
K, — K,, is too small a fraction of K, for reliability, and the 30% Bi specimen which 
belongs to a different crystal system. ‘The fact that all the curves coincide at 
T =T,, is simply a consequence of the way in which we have calculated the lattice 
thermal conductivity in the normal state. All the curves rise with falling 
temperature, and appear to reach a maximum at about 0-47,. Assuming 
tentatively that X is zero and that K,—K,, is indeed to be interpreted as K,,, 
we have for more convenient discussion plotted in fig. 9(a) and (4), instead 
of K,,, W=1/K,, against T/T,. The resultant curves for all five specimens 
are seen to be very much alike. 

In fig. 9(b) we have also shown the terms W, and W) due to scattering at 
grain boundaries and impurities calculated for the 10% Bi specimen in the 
normal state. A comparison of the total observed W with these two contributions 


shows that they are of the same order of magnitude at the lowest temperature of 
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measurement. This would seem to indicate that the only other contribution 
to W, the scattering of the lattice vibrations by electrons, should be negligible 
at this temperature. Thus the high values of the observed thermal conductivity 
might indeed be explained by a rapid decrease of the scattering of the lattice 
vibrations by the electrons in the superconducting state. If we consider such 
a process possible it is of course tempting to try to derive the temperature 
variation of the scattering of phonons by electrons in the superconducting metal. 

Assuming this to be the same function for all the specimens we obtain 
Wag%(T/T,)”, where n is about 4. When we take into account the uncertain 
validity of eqn. (3) in the region of the minima in K, then the results are probably 
in agreement with this power law to within the limits of error, We should 
certainly beware of placing too much reliance upon a conclusion reached after 
such a long chain of uncertain calculation. The accuracy of the interpolation 


50 $0 


deg cm/ watt 
deg cm/watt 


(0) 


T/T. 


Fig. 9. 1/Kgs for specimens with various percentages of bismuth. Ws and Wp are calculated 
for specimen containing 10% bismuth. 


formula for K,,/K,,, in particular, has not been very soundly verified. If, 
however, we accept the above result, then we find 

Wr Mae.) (8) 
which we may compare with the fifth-power law for K,,/K,, suggested by 
Hulm (1950) for the case when the electrons are scattered only by lattice 
vibrations. 

The lattice resistivity W of the specimen containing 30% bismuth is shown 
in fig. 10. In considering this diagram it must be borne in mind that this alloy 
belongs to a different crystal system from that of the other lead—bismuth alloys 
examined. In addition we must remember that the amount of impurity scattering 
will be greater than in these specimens. With these allowances the behaviour 
appears in at least qualitative agreement with the explanations given in the 
previous section. 

(111) Comparison with other Results 


It is of interest to see if the picture given above and deliberately fitted to our 


results for the lead—bismuth alloys with bismuth contents below 10 °% will fit other 
results as well. 
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The previously published (Mendelssohn and Olsen 1950a) thermal 
conductivity results for niobium did not contain reliable values for the normal 
conductivity, but these may be deduced from electrical measurements. A value 
for Ky/ Ken derived from these data is shown in fig. 11, and the corresponding 
lattice resistivity is shown in inset. The agreement with the treatment given in 
the previous section is only partly satisfactory. An extrapolation of the lattice 
resistivity curve according to a 74 law would lead to a value of 1 000 deg cm/watt 
for the lattice resistivity at the transition temperature, and this appears to be 
approximately ten times greater than the value to be expected on the basis of 
Makinson’s (1938) theory. 

The fact that the alloy containing 30° tin does not behave as the bismuth 
alloys and does not have an abnormally high value for the thermal conductivity 
in the superconducting state naturally requires some explanation. There would 
appear to be two reasons: the electronic thermal conductivities of the alloys 
containing tin are considerably greater than of those containing an equal impurity 
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of bismuth which belongs to a different group in the periodic table; in addition 
to this, while the electrons are less affected by the presence of tin than of bismuth, 
one may expect that because of the greater difference in atomic mass the lattice 
vibrations will be more strongly scattered by tin than by bismuth impurity. 
The situation thus arises in this alloy, that although only a fairly large increase 
in K, would be noticeable, there is too great a resistance due to impurity 
scattering of the lattice vibrations to allow K,, to be increased much above K,,. 
To make these reasons completely convincing, however, it is clearly desirable 
to have much more information on the heat conduction in other superconductive 
alloy systems. 

An examination of the results obtained by Hulm (1950) for the lattice 
conductivity in superconducting tantalum, and in tin of high impurity content 
appears to indicate that perhaps the fourth-power variation with temperature 
of K,, which we have deduced for lead is not a general feature; or else, and 
perhaps more probably, that we have obtained too high a power for the 


temperature dependence of K,,, 
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§4. FIELD TRANSITIONS 


Further information on the conduction mechanism can be gained by 
measuring the thermal conductivities in the interval of temperature and magnetic 
fields in which the actual destruction of superconductivity takes place — 
(Mendelssohn and Pontius 1937). Such measurements have been carried out 
in longitudinal and transverse fields for the alloys of the lead—bismuth series. 

In fig. 12 are shown conductivities of our two purest specimens in a 
longitudinal field. It will be seen that just as the transitions measured by 
Mendelssohn and Pontius (1937) on lead, and by Hulm (1950) on tin, these are 
fairly sharp, and that any lack of sharpness may be accounted for by the 
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Fig. 12. Thermal conductivity in longitudinal magnetic field for pure lead specimens. 


demagnetization factor of the specimens, which are not infinitely long. The 
observed width of the transition at 4-6°K is somewhat greater than the width 
of the electrical transitions found by MacDonald and Mendelssohn (1949) in 
long specimens at the same temperature. At lower temperatures where these 
authors found very wide electrical transitions, the thermal conductivity changes 
as sharply from its superconducting to its normal value as at the higher 
temperatures. [or comparison we have also shown the shape of a transition of 
the electrical resistance p in the diagram for 2:7°x. It is plotted so that p=0 
coincides with K, and p=pyorma coincides with K,. The greater width of the 
electrical transition, and the fact that the transition in the thermal conductivity 
remains sharp just as the magnetic transition, would appear to lend support to 
the view that the wide electrical transitions observed in lead are due to the 
existence of only very small regions having high critical field. 

In fig. 13 we show the transition of a pure specimen in a transverse field at 
3:3°K. We see that the transition in conductivity is far from linear. If the 
same transition is plotted in terms of thermal resistance, on the other hand, it is 
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found that the variation is almost linear with field. A very simple calculation 
shows that this is an indication that the laminae forming the intermediate state 
(Landau 1943) lie at right angles to the axis of the cylindrical specimen. The 
argument leading to this result is as follows: If the fraction of the volume of the 
cylinder occupied by normal and superconducting material are a and b 


respectively, then the conductivity of the cylinder will be given for the case 
of longitudinal laminae by 


WSs (aks WOK.) K@s:3)] came Mirae, Tae. oat (9) 
but if the laminae lie perpendicular to the axis of the cylinder, then 
1/K =(a/K,,+6/K,)/(a+6). ‘ © .....- (10) 


It is known from magnetic experiments (Shoenberg 1948, etc.) that once $H, 
has been exceeded for a cylinder in a transverse field, then the amount of normal 


material increases approximately linearly with field until all the material is 
normal at H... 
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Fig. 13. Transverse field transition at 3-3° k for pure lead specimen B. 


An anomalous effect indicating a rather more complicated behaviour has 
been found in specimens with a small degree of impurity or strain (Mendelssohn 
and Olsen 1950a,c). Similar results have recently also been obtained on 
pure lead, particularly at low temperatures and under complicated geometrical 
conditions (Webber and Spohr 1951). 

In fig. 14 is shown the variation of thermal conductivity with field for our 
specimens of intermediate purity. At the highest temperatures there is a 
.transition which is very like that observed in pure specimens. ‘The departure 
from linearity is less obvious than in these, but this is simpiy a consequence of the 
much lower ratio of K, to K,, At asomewhat lower temperature there is a sharp 
drop in conductivity when $H, is exceeded, and the conductivity then recovers 
and rises to its value in the normal state. At a still lower temperature only the 
initial drop is found, and the recovery is almost non-existent. ‘The return curves 
in a decreasing field all lie below the increasing field curves. 

This anomalous drop in the conductivity between the completely normal 
and the completely superconducting state creates a difficulty when an explanation 
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is attempted. As we have already pointed out when announcing the first of these 
results, it is impossible to interpret the results simply as being due to a set of — 
successive laminae of alternating heat conductivity. A series of thermal 
resistances of which the individual components have either the normal or the 
superconductive resistance must clearly always result in an overall value 
intermediate between these two resistivities, as for example in the case in fig. 13. 
A heat conduction which is lower than in either of these states cannot be 
explained in this simple way, but requires some more complicated mechanism, 
such as for example the circulation process referred to previously. 

The only other explanation for the heat conduction results in superconductive 
alloys is that of the enhanced lattice conductivity, and it is therefore interesting 
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Fig. 14. Transverse field transitions for specimens containing 0-1, 0:2 and 0:5% bismuth. 


to see what additional assumption would have to be made for this alone to explain 
the experimental results. It is interesting to note that by calculating the lattice 
and electronic heat conductivities separately, and treating them as parallel 
throughout the length of the specimen, a minimum in the heat conductivity in 
the intermediate state can be produced. This is equivalent to assuming that the 
two heat currents flow independently without transfer from one to the other. 
From what has been said in the preceding section we should expect the overall K, 
to decrease rapidly as superconductivity is destroyed, while at the same time 
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K, increases. Little can be said, of course, concerning the way in which these 
changes depend on the field, but as fig. 15 shows a minimum can be produced 
by a suitable adjustment of the constants. 

The trouble with this picture is that it appears to contradict our concept 
of an intermediate state made up of a set of laminae each with their own value 
of the thermal conductivity. However, it serves as a useful guide to show what 
assumptions would have to be made to explain the phenomena by an enhanced 
lattice conduction. It seems clear that in order to obtain an approach to the 
case of parallel conduction processes (a concept which is incidentally similar to 
the circulation hypothesis) there must be little interchange of energy between the 
lattice and the electronic heat currents over a distance of the order of magnitude 
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Fig. 15. Expected variation of Ke and Kg during superconducting transition with zero energy 
transfer between the two processes. 


of the thickness of the laminae. ‘Thus we would require the periodicity of the _ 
laminar structure to be of the same order or smaller than the mean free path of 
the lattice vibrations or electrons. 

Summarizing the work described in this paper, we see no easy and simple 
way of explaining the variety of the observed heat conduction phenomena. The 
validity of such ideas as a circulation process orthe enhanced lattice conductivity 
must still rely on assumptions which can only be checked by further experiments. 
The present work shows that the second as well as the first of these processes. 
might be capable of explaining the phenomena, and that so far neither process. 
can be ruled out by the evidence provided. Beside the recording of the 
anomalous phenomena associated with heat flow in superconductive alloys, 
experiments of this kind may offer a useful new method for investigating the 
intermediate state. 
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Infra-Red Absorption by Metals at Low Temperatures 
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ABSTRACT. The infra-red absorptivities of several electropolished metal surfaces 
have been measured at liquid helium temperatures with a specially designed radiation 
‘cryostat with a black body at room temperature as the source and a differential helium 
gas thermometer as the detector. It is found that the absorptivity of tin at 1-95° k for the 
mean estimated wavelength 14 4 remains the same in the normal and superconducting 
states to within a probable error of 0:3%. The absorptivities of the normal metals, which 
include a number of alloys, are compared with the values calculated by various theories of 
metallic reflection. [he classical Hagen—Rubens formula is found to be valid for the 
alloys except where there are complications such as the internal photoelectric effect. The 
absorptivities of the pure metals are in fair agreement with the values calculated by the 
anomalous skin effect theory neglecting relaxation effects. 


$1) INTRODUCTION 


HE Classical electromagnetic theory of light leads to an expression for the 
absorptivity A and reflectivity R of a metal surface at normal incidence 
A=(l—R)=2(vlet oe os yl ee (1) 

where v is the frequency of the radiation and o the d.c. conductivity of the metal 
in electrostatic units. This suggests the interesting possibility that super- 
conducting metals might be particularly good reflectors of high-frequency 
radiation. In the infra-red region, Daunt, Keeley and Mendelssohn (1937) 
carried out absorption measurements on lead and tin using the radiation from a 
black body at room temperature as the source. Estimates of the relative 
absorption in the normal and superconducting states of the metal were made by 
measuring the rates of evaporation of liquid helium contained in an attached 
calorimeter. [From these observations it was concluded that the transition to 
superconductivity does not affect the infra-red absorptivity of a metal. 
Experiments in the microwave region were first made by H. London (1940) 
who showed that the absorptivity of a metal, when cooled through its transition 
point, does not fall suddenly to zero as indicated by (1) but rather rapidly at 
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first and then more gradually, apparently tending to zero at the absolute zero: 
of temperature. More recently Pippard (1947a, 1950) has used a resonance 
technique to carry out extensive investigations of this phenomenon at 1 200 and 
9400 Mc/s, and has confirmed the results of London by this independent 
method. ‘These results have, during recent years, aroused considerable interest 
in the behaviour of superconductors at still higher frequencies and also made 
it desirable to approach the problem from the infra-red end of the spectrum. 
The scope of the investigation at infra-red frequencies is at present necessarily 
limited to the shorter wavelengths owing to the lack of monochromatic sources. 
of adequate intensity for wavelengths longer than about 20y. However, 
extrapolation from the work of Pippard suggests that even at such short 
wavelengths as 10, the absorptivity of a metal might show a small decrease of 
the order of a few per cent in the superconducting state. Since the work of 
Daunt, Keeley and Mendelssohn did not claim a high degree of accuracy, it was. 
decided to make new measurements on the infra-red absorption by tin with an 
apparatus designed to give results accurate to better than 1%. 

The high-frequency study of metals, besides throwing light on the 
phenomenon of superconductivity, has also led to interesting results about the 
behaviour of normal metals. During his investigation of the behaviour of tin 
in the microwave region, H. London observed that the surface resistivity of the 
metal above its superconducting transition point was higher than that predicted 
by the classical theory. He attributed the anomaly to the fact that at low 
temperatures the mean free path of the conduction electrons is either comparable 
with or much larger than the depth of penetration of the electric field into the 
metal. ‘This observation by London has also been confirmed by Pippard’s. 
(1947 b) measurement of surface resistivities of copper, silver, gold, aluminium, 
tin and mercury and particularly by the later detailed study by Chambers (1950 a) 
of copper and silver. From simple theoretical considerations Pippard predicted 
that the surface resistivity of metals at microwave frequencies should vary as. 
the two-thirds power of the frequency. ‘The detailed theory, taking account 
of the variation of the electric field along the path of conduction electrons, has. 
been worked out by Reuter and Sondheimer (1948) using the general methods 
of the theory of metals. The results of this theory seem to offer a satisfactory 
explanation of all observed data on normal metals in the microwave region. 
The verification of the predictions of the theory regarding the behaviour of metals. 
in the infra-red region has not so far been possible due to lack of experimental 
data. Therefore it was also the aim of this investigation to measure quantitatively 
the absorptivities of some pure metals and alloys and compare the results with 
those predicted by various theories, including that of Reuter and Sondheimer. 


§2. EXPERIMENTAL TECHNIQUES 
(i) Method of Investigation 
In the present method, as in that of Daunt, Keeley and Mendelssohn, the 
absorption by the specimen of radiation from a source at room temperature is 
studied. The radiation from the source passes down the cryostat through a 
German silver tube T (fig. 1), polished on the inside to increase the effective 
intensity by multiple reflections at the wall, and falls on the flat polished surface 
of the metal specimen S under investigation. ‘The energy absorbed is measured 
by the heating of one of two identical gas thermometer bulbs to which the 
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specimen is soldered and which are enclosed in a highly evacuated metal chamber. 
The gas thermometers are connected by metal capillaries to either side of a 
miniature oil manometer at the top of the cryostat. ‘The lower portions of the 
tube T, the capillaries and the metal chamber are immersed in liquid helium ~ 
during an experiment. When the shutter Sh (which is at a low temperature and 
therefore does not itself radiate appreciably) is opened, radiation falls on the 
specimen and the increase in pressure in the thermometer bulb to which it is 
attached causes a displacement of the levels in the oil manometer, which 1s 
measured by a travelling microscope. The power absorbed by the surface is 
obtained by comparing the change in oil level caused by the radiation with that 
produced by dissipating a known electrical power in a heating coil H wound 
round the gas thermometer carrying the specimen. In this way knowledge of 
the absolute temperature rise is unnecessary, and various correction factors are 
eliminated. In order to observe as accurately as possible whether there is a 
change in absorption when the specimen changes from the superconducting to 
the normal state, the cryostat was surrounded by a pair of Helmholtz coils which 
could give a magnetic field sufficient to destroy superconductivity. 
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Fig. 1. Experimental arrangement. 
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The effective intensity of the source, which is required to calculate the 
absolute absorptivities of the various metals investigated, was estimated by 
measuring the absorption by a surface thickly coated with zinc black powder, 
known to absorb the incident radiation almost completely in the spectral region 
in question. The value thus obtained was checked independently by 
measurements on alloys whose absorptivities could be estimated from the 
classical formula (1). 
(11) Experimental Details 

The factors considered in designing the gas thermometer were its rise in 
temperature above that of the helium bath, the time taken to attain equilibrium 
when exposed to radiation and the effect of dead space on the sensitivity. The 
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rise in temperature is proportional to the thermal resistance of the capillary 
connecting it to the bath and the incident radiant flux which is independent of 
the detector itself. The equilibrium time is proportional to the thermal 
resistance of the capillary and to the total heat capacity of the thermometer 
(including the specimen and gas), which was kept as low as possible since it does 
not influence the sensitivity. The capillaries were made of an 80: 20 cupro-nickel 
alloy with dimensions chosen to give a reasonably short equilibrium time (about 
3 minutes) and a sensitivity of the order of 10~* watt per mm of oil at 4-2°k for 
a gas pressure of a few cm Hg. ‘The narrow capillaries (0-5 mm) used at first 
to minimize the dead space were frequently blocked by impurities in the helium 
and so had to be replaced by wider ones (15mm), in spite of the consequent 
reduction of sensitivity by about 50%. It was not possible to use the gas 
thermometer below about 2°k, since the rapidly falling vapour pressure restricts 
the working pressure too greatly and makes the thermometer too insensitive. 

The metal shields 5, and S,, soldered to the bottom of the vacuum chamber, 
prevent scattered radiation from falling on the sides of the thermometers. 
‘These shields decreased the apparent absorption from about 350 to 140 erg/sec, 
but introduction of additional shields carefully designed to eliminate the small 
influx of radiation through the gap below the specimen caused no further 
improvement. The consistency of the results obtained with different specimens 
of the same metal was taken as evidence that the shields S, and S, were performing 
their task efficiently. The thin (} mm) silver chloride* window W, fixed above 
the specimen, is used to prevent ‘hot’ molecules of the residual gas from heating 
up the thermometer on opening the shutter. One end of the heating coil is soldered 
to the gas thermometer while the other end is taken out through a glass—metal 
seal N, and then through a thin-walled German silver tube in order to shield it 
from extraneous high-frequency radiations such as that from a neighbouring 
cyclotron. A high vacuum was maintained in the chamber during the 
experiments by first evacuating it at room temperature with an oil diffusion 
pump, then isolating it from the pump by a tap and relying on cooling in liquid 
helium to freeze out the residual gases. 


(iii) Specimens 

The specimens listed in table 1 were made as circular discs of about 3mm 
thickness and $in. diameter. Their electrolytic polishing was carried out by 
well-known techniques (Kemsley 1948). They were all soldered with Wood’s 
metal symmetrically to the thermometer. The electrical resistivities of the 
specimens (which are needed in the theoretical discussion) were, for the alloys, 
obtained from measurements on strips cut out of the specimens after the 
absorption measurements were finished and, for the pure metals, were taken 
from published data. 


(iv) Calbration and Performance 


To test the constancy of the source of radiation, the absorption by a copper 
surface was measured in the range from 4:2° to 2°K. Since theoretically there , 
is no obvious reason why the absorption should vary, any apparent variations 

* At room temperature the transparency of a thin silver chloride window is known to be 


uniformly about 70% in the range 1 to 40 4; at low temperatures the transparencies of solids in 
general increase and extend to longer wavelengths. 
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must be attributed to variations of the effective intensity of the source. The 
result of this experiment, illustrated in fig. 2, indicates a gradual increase in 


absorption of about 6 or 7%. Detailed examination of this variation in terms of — 
the level of helium in the bath and the rate of pumping suggested that the — 


effective intensity depends in a complicated way on the distribution of 
temperature in the top portions of the radiation tube and therefore is not exactly 
reproducible from run to run. When the absorption was not expected to vary 
with temperature, the mean of all the readings in any one experiment was taken 
as the appropriate value and this should be accurate to about 3%. This’ was. 
confirmed by the consistency of results obtained on the same specimen in 
different runs. 


Table 1. List of specimens 


Specimen 


Electrolytically pure copper 


Electrolytically pure copper 


Tin J.M. 2356 (<.0:004°%, impurity) 

Lead J.M. 560 (<<0:005°% impurity) 

Aluminium from British Aluminium Com- 
pany (<0:01% impurity) 

Tin J.M. 2356 alloyed with 1% and with 
5-4°% chemically pure indium 

57% Kahlbaum bismuth+43°% chemically 


Method of preparation 


Machined from rod, annealed and electro-- 
polished. 

As above but polished with emery paper 
then with ‘ Brasso’ instead of electro-- 
polishing. 

Cast against glass and electropolished. 

Cast against glass and electropolished. 

Cut from sheet, annealed and electro-- 
polished. 

Cast against glass and electropolished. 


Cast against glass. 


\ 
| 


pure tin 
Brass 
Eureka and Staybrite 


Cut from sheets, ground with sand paper 
and electropolished. 


Note: J.M.=Johnson Matthey. 


} commercial stock 


Temperature (°K) 


Fig. 2. Variation of apparent absorption A by copper surface expressed as a ratio A/Ay where 


Ay is the lowest observed value. 


The absolute value of the effective intensity of the source, which was. 
required to estimate the absorptivities of metal surfaces, was determined by 
measuring the power absorbed by each of two surfaces, one blackened with 
zinc black and the other with Aquadag; these were found to be 11000 erg/sec 
and 8000 erg/sec respectively. Since at room temperature zinc black is known 
to absorb almost completely and Aquadag from 60 to 70% of the incident 
radiation in the infra-red region under consideration, these measurements may 
be regarded as reasonably consistent and the incident power was assumed to be 
11000erg/sec. ‘Taking into account the variation of intensity during each 
experiment, this estimate of the effective intensity is probably accurate to only 


about 3% and the final absorptivities quoted in column (3) of table 2 to only 


about 5%. 
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The mean wavelength to which the absorption measurements refer was 
deduced using Planck’s formula for the spectral distribution of intensity in black 


body radiation. If the incident intensity is given by | E,, dv, the total energy 
0 


absorbed by the metal surface can be represented by 


| cv“h, dv, 


0 
where c and a are independent of the frequency v. The experimentally measured 


absorptivity is the quotient 
| evtE, dv : | E,dv 


0 0 


which may be regarded as the absorptivity cv,“ at some mean frequency 1. 
It is found by integration that the corresponding mean wavelength A, does not 
vary appreciably with a, being 14-04 when a=4 and 13-6 when a=3. 


§3. EXPERIMENTAL RESULTS AND DISCUSSION 
(1) Superconductors 

Although a number of superconductors such as tin, lead and alloys of tin 
with indium and bismuth were investigated, attention was mainly focused on tin. 
The absorption measurements on tin in the range 4:2° to 2°K do not indicate 
any change corresponding to the transition to superconductivity at 3-7°K but 
only the gradual apparent increase in absorption which, as with the copper 
surface, could be attributed to an increase of the effective intensity of the source. 
Thus this experiment indicated only that the absorption did not vary by more 
than about 5% down to 2°x. To achieve the higher accuracy permitted by the 
detector the cryostat was cooled to a temperature (1:95°K) at which the detector 
still works satisfactorily and the absorption was measured at regular intervals 
of time alternately in the superconducting state and in the normal state produced 
with the aid of a magnetic field. As can be seen from fig. 3 the two sets of readings 


o Normal:magnetic field on 
° © Superconducting: magnetic field off 


Absorption (erg/sec) 


oe aL 
4 5 


3 
Time (hours) 


Fig. 3. Absorption by a tin surface in the normal ‘and superconducting states at different times, 


corresponding to the normal and superconducting states plotted against time lie 
on the same curve within the accuracy of the detector, indicating that the 
absorptivity of tin is the same in both states. The probable uncertainty in the 
constancy of absorption was estimated as less than 0-3°% from the deviations of 
one of the two sets of points from a smooth curve drawn through the other. 
During the experiments it was noticed that switching the magnetic field on 
and off, besides causing an eddy current heating of the thermometer (and the 
specimen in the normal state) which gave rise to a kick in the oil levels of about 
1mm, produced a kick of about lcm, caused by the latent heat accompanying 
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the transition of the specimen from the normal to the superconducting state or 
vice versa. The points shown in fig. 3 were obtained after waiting the necessary 
time for thermal equilibrium to be restored. 

In the case of the other superconductors the magnetic method was not used 
and hence the constancy of their absorption can be inferred only to within 
about 5%. For lead this is not particularly significant since the temperature 
range investigated is well below the transition point 7-2°K and no large change 
with temperature was to be expected; for the alloys, however, an appreciable 
change might have been expected because of their high residual resistance in 
the normal state. . 

Using dimensional arguments Pippard (1949) deduced that just below the 
transition temperature the gradient dr/dT in the resistivity-temperature curve 
should vary as »~?°(r=R/R, where R and R, are the surface resistivities in the 
superconducting and normal states respectively) and this was confirmed’ by 
observations on tin at the two frequencies 9400 and 1200 Mc/s. At 9400 Mc/s 
dr/dT at the transition temperature is about 20 and it should therefore be about 
0-12 for the infra-red frequency used, if extrapolation according to the v-?? law 
is permissible over such a wide range. Even though of course dr/dT becomes 
smaller as the temperature falls below the transition temperature, this would 
suggest that a change of a few per cent in 7 and therefore in absorptivity might 
be expected when the specimen becomes superconducting at a low temperature. 
The reason why such a change is not in fact observed is probably because the 
theoretical argument is based on the ‘two-fluid’ model which assumes the 
coexistence of ‘normal’ and ‘superconducting’ electrons. At these high 
frequencies relaxation effects in the normal electrons should probably also be 
considered, but this would lead to an even more pronounced fall of resistance 
in the superconducting transition. However the two-fluid model already 
disagrees with some aspects of Pippard’s observations at 9400 Mc/s and might 
be expected to break down altogether at still higher frequencies due to the 
setting in of new absorption mechanisms. It is possible that the frequency at 
which such mechanisms set in is of order of magnitude RT,/h (where T, is the 
transition temperature, Rk the Boltzmann constant and h Planck’s constant) 
which corresponds to the rather inaccessible region of wavelengths of a few 
millimetres. 

(11) Normal Metals 

In table 2 the experimental absorptivities of the different metals are compared 
with those calculated by various theories (columns (4) to (8)). The formula of 
classical theory has been used in column (4) while in column (5) the absorptivity: 
is assumed to be given by 2(v/c)!?{—wr+(1+w?r?)!2\2 where w=2nv and 
7 the ‘relaxation’ time is equal to omjne®, n being the number of conduction 
electrons per cm*, e the electronic charge and m its mass. The calculations in 
columns (6) to (8) are based on the theory of the anomalous skin effect in metals 
developed by Reuter and Sondheimer. The parameter p in these columns is 
the fraction of the conduction electrons reflected specularly on collision with the 
metal surface. Chambers (1950 b) has shown that his results on the electrical 
resistivities of thin wires are in agreement with a value of p=0. For the optical 
region Reuter and Sondheimer have made their computations only for p=1. 
The values in columns (6) and (7) are based on fig. 1 of their paper which 
represents their results if relaxation effects are assumed negligible. The values 
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in column (8) are obtained by interpolation from fig. 2 of their paper which 
represents their results for the optical region where relaxation effects play an 
important part in determining the surface resistance. In all these calculations, 
the number of free electrons per atom has been assumed to be 1; the numerical 
values are not much changed if is changed by a factor of 2. It may be otend 
that these assumptions if used to calculate the resistivities at microwave 
frequencies give fair agreement with the available experimental data .(Pippard 
1947 b, Chambers 1950.a, Maxwell, Marcus and Slater 1949), except where 
a poor surface was used. 

It can be seen that the calculated absorptivities assuming ve 0 and no 
relaxation effects (column (6)) are in fair agreement with experimental values 
while those in other columns, for most cases, lead to absorptions lower than 
observed, this discrepancy being more marked in the case of the pure metals 
and the 1% indium alloy for which the resistivities are low enough to produce 
the anomalous skin effect. These features suggest that relaxation is unimportant 


Table 2. Experimental and theoretical percentage absorptivities at A=14-0 yu. 
(1) (2) (3) (4) (5) (6) (7) (8) 


d.c. resistivity — Theoretical —_—_—_ 
Metal «10%ohmxcm _ Experi- Classical Classical+ Anomalous Skin effect Anomalous-+ 
at 4-2°K mental theory relaxation p=0 p=1 relaxation 

Copper 0-011 0-62 0-10 0-003 0-94 0-84 O-1 
Copper (M.P.) 0-011 1-47 0-10 0-003 0-94 0-84 0-1 
Aluminium 0-0161 key ti 0-13 0-004 1-03 0-91 0-08 
Lead 0-0034 1-15 0-06 0-0005 1-18 1:05 0-001 
Tin 0-0021 1:24 0:04 0-0004 1-14 1-02 0-01 
Sn+1% In 0-44 11-25 0-65 0-08 1:26 1-07 O-1 
$Sn+5-4% In 2°82 1-74 1-64 0-18 1-90 1-68 0-19 
Brass 3°57 1:78 1:84 0:96 104 ea ane 0-96 
Eureka 44-9 6:32 6°33 5-62 6:33 6°33 5-62 
Staybrite 47:8 8-23 0258 5:66 6°53 6:53 5-66 
Bi-Sn eutectic (C.S.) 7-08 7:63 PIII 1:20 2:69 PONS 1-20 


M.P.=Mechanically polished with ‘ Brasso’. C.S.=Cast surface. 


even at 14 in all these metals except perhaps for copper which approximates 
better than the others to the free electron model on which the theory is based. 
That the values in the last column are not appreciably influenced by the value of 
chosen can be inferred from the fact that, even for copper for which the 
disagreement with experiment is least, would have to be increased five or 
six times to reduce 7 sufficiently to give an absorptivity in agreement with 
experiment. Extension of the theory to the case p=0 might explain why 
relaxation effects are unimportant at 14,, but it is possible that some more 
fundamental revision of the theory may be required. 

The disagreement with the classical theory shown by the alloys Staybrite 
and the Bi-Sn eutectic is significant. Pure bismuth is known to have a high 
absorptivity (20°%) at room temperature (theoretical value ~10°%) and this 
suggests that the higher observed absorptivities of these alloys are due to 
phenomena such as the internal photoelectric effect. ‘he agreement with 
theory shown by the other alloys, brass, Eureka and the 5-4°% indium specimen, 
and the agreement of the absorptivity of Eureka with that (5-8°%) deduced by 
Hagen and Rubens from reflectivity measurements at room temperature (the 
electrical resistivity of Eureka is practically independent of temperature over 
the range under consideration) seems to confirm the estimate of the effective 
intensity of the source. 
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Finally, the high value of the absorptivity of the ‘Brasso’-polished copper, 
which is not very different from that (1:29%) obtained by Wexler (1951), 
emphasizes the importance of using electrolytically polished surfaces for such 
optical measurements. This is in agreement with the conclusion of Chambers. 
that a layer of high resistivity is produced at the surface of a mechanically polished 
surface. 
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ABSTRACT. ‘The Hamiltonian operator for an exchange antiferromagnet is expressed’ 
in terms of the deviations of the spins from the basic ordered antiparallel arrangement. 
By a Fourier transformation the Hamiltonian may be approximately represented as a 
sum of harmonic oscillator terms, corresponding to spin waves. The role of anisotropy 
energy in fixing the domain axis is discussed, and the energy levels are calculated for the 
case of an external field parallel to this axis. The ground-state energy includes the 
zero-point energy of the spin waves and is in good agreement with other theories. 


Formulae are quoted for a Hamiltonian including dipole-dipole interaction between 
non-isotropic spins. 


$1. INTRODUCTION 
Nk of the basic methods for the theoretical treatment of ferromagnetism 
is that of the ‘spin wave’, a concept due to Bloch (1930, 1932). 
Starting from the basic ordered state of the lattice, spin deviation 
operators are defined to describe the disordered states, and the Hamiltonian 
is expressed in terms of these. By a Fourier transformation of these deviation 
operators (and after making certain approximations) the Hamiltonian can be 
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made to separate into a set of independent harmonic oscillator terms, each term 
corresponding to a wave of spin deviations propagated through the crystal by 
the exchange interaction between neighbours. In this way Bloch was able to 
deduce the temperature variation, and Holstein and Primakoff (1940) the field 
dependence, of the intrinsic magnetization. 

The particular value of this method is that it is properly quantum-mechanical, 
and does not make the rather drastic assumption that the spins may be replaced 
by Ising dipoles. Moreover, it is not restricted to nearest-neighbour 
interactions—even long-range dipole-dipole interactions may be considered 
explicitly. The only limitation to its application is that the mean value of the 
spin deviation should be small; it is essentially a low-temperature theory. 

The problem to be considered here is the possibility of extending this 
powerful technique to the case of antiferromagnetism, a term which is generally 
used to describe the phenomenon whereby certain substances tend to fall into 
ordered states in which the spins of neighbouring atoms are antiparallel 
(Van Vleck 1941, Grenoble Conference 1951). Such a state would seem natural 
for a substance in which exchange forces favour antiparallel spins on neighbouring 
or next-neighbouring sites; it has also been shown that for classical dipoles the 
antiparallel state may be the lowest, depending on the crystal structure and the 
shape of the specimen (Sauer 1940, Luttinger and Tisza 1946). 


§2. THE ANTIFERROMAGNETIC LATTICE 


‘The specimen is to consist of a crystalline lattice whose N sites are located 
by the values of the vector 
R,=lu,t+hu,t+iu, .seeee (1) 


where J,, J, /3, are integers and u,, uy, us are the basic vectors of the lattice. 
On the /th site there is to be found a spin with operator S, of total spin S (in units 
of h). Between the /th and mth spins there is an exchange interaction 
—2JimS1+S, depending on the magnitude of the relative distance operator 
R,,=R,—R,,.. Thus the major term in the total Hamiltonian of the system is 


Hy, ay ey, m2I mS, : Sin teens (2) 


The problem of finding even the lowest eigenstate of #, is exceedingly 
difficult when /J,,, is negative. But the lowest states for classical or Ising dipoles 
with the given Hamiltonian are known, and it is at least plausible that the 
eigenstate is derived from this by a small perturbation. In the ferromagnetic 
case the same assumption is made (quasi-saturation), although here also the 
completely saturated state .s not an eigenstate of the complete Hamiltonian, 
owing to the magnetic interaction terms. ‘The whole question has been 
considered in detail by Anderson (1951a), using the Kramers—Heller spin-wave 
technique, and the assumption is thought to be justified. 

To describe this basic state we must introduce a special vector w such that 
exp (7w . R,) alternates in sign from one lattice point to the next. An expression 
of this sort is sufficient to describe the antiparallel ordering due to strong 
nearest-neighbour interaction in the simple cubic and body-centred cubic 
lattices, but in other cases it may be necessary to take a linear combination of 
such terms (Luttinger 1951). This does not lead to any essential difficulty, but 
it will not be considered here because it introduces unnecessary complicaticns. 
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The vector w is best defined in terms of the vectors reciprocal to u,, Uy, U3, 
Le. W=W,V1+ WyVo + W3V3, where 
Uy X Us u, x U, ie Ue 
V7,= —_—_., v.= —— _., = ——_—_ .. 
u, . Uy X Us U,. Uy X Us a Rea 
In this way, R,. w =1,w, + /,w, +/,w3, and in the typical case of the simple cubic 
lattice with next neighbour interactions w=7(i+j+k). 

Besides this alternating function, we also need to prescribe a direction 
along which the spins would be aligned in order to fix the basic state exactly. 
It is usual to assume that a crystal of an antiferromagnet consists of domains, 
in each of which the spins are parallel or antiparallel to a given direction, this. 
direction varying in a random manner from domain to domain. Yet the operator 
H » is isotropic, and there is no recognizable direction of alignment at any point 
of the crystal. When a field is applied externally, this provides a preferred 
direction, and the whole system will be affected equally, swinging round until 
it is in a position of minimum potential energy and destroying any assumed 
domain structure. In order to overcome this difficulty, it is necessary to 
assume that there is an additional anisotropy energy fixing certain directions 
of easy magnetization (Nagamiya 1951, Kittel 1951). 

It is not easy to decide the origin of this term. If it is to be correctly 
described, with, for example, cubic symmetry in the cubic lattice, it would 
contain terms of the fourth degree in the spin operators and lead to very 
complicated expressions. The best that can be done for the present is to 
introduce into the Hamiltonian a term 


H ,=—YjgBSPH, exp(iw-R), = «2... (3) 


that is, a fictitious internal field staggered so as always to be parallel to the spin 
in the basic Ising state.* 6H, need only be very small compared with /,,, in 
order to make sure that the spin system as a whole maintains its direction. This 
direction, the z axis, will evidently be one of the crystallographic axes of the 
crystal. A more general condition for domain stability will be discussed later. 

To calculate the parallel susceptibility we also need a magnetic field H, 
(due to external sources) along the z axis, making the total Hamiltonian 


He = — Bd 2S eq Sy0 Sin — Te BS{{H, + Hy exp-(iw.R,)}.  - «2... (4) 


This expression does not include all possible types of interaction, but covers. 
those which have been generally considered in the literature. The generalization 
to non-isotropic spins is discussed briefly in the Appendix. 


§3. SPIN-DEVIATION REPRESENTATION 

We shall use a representation in which the S;~ components are diagonal. 
The wave function for the whole lattice may then be considered as a linear: 
combination of all the various product wave functions, one from each site, 
which can _be labelled by the values of S,{%. Let us for the moment consider 
only the sites for which the value of S; in the assumed basic state is positive, 
i.e, those which ‘point to the right’ in the lowest-energy configuration. 
We define the following operators: 


St=S,9+iS, S-=SO-iSM, w=S-S@  ...... (5) 


* As usual, B is the Bohr magneton, and g the Landé factor for the spin. 
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n, is the spin deviation operator. ‘The states, which were formerly labelled by 
values of S{ may now be relabelled (on the right-hand sites, that is) by the 
eigenvalues », of n,. Thus 


n, |S). ee aa . . Sy) =n, |2, . Ceulhieewetny) 


okey 1/2 
But Si*[m)=(28)!# (1— "SAY" (a= 1) 


S 
Sy |2,) =(2S)#2(n, + 1421 —n,/2S)2[m,+1) = nace (6) 
and if we introduce creation and destruction operators a,, a,*, where 
a,* |n-=(m,+1)¥?|n,+1), a,|2,)=(m,)¥?|m,—-1) — ...-.. (7) 
we find Nj =a,"a, | (8) 
St =(2S)P(1 —a,*a,/2S)'?a,, Sp =(2S)?a,*(1 —a,*a,/2S)'? 


For the left-hand spins, i.e. those for which exp (‘w.R,)=—1, we define 
our spin deviation operator in the opposite sense: 


Hees tio car, Se oe ee. (9) 

The states can be relabelled in terms of the eigenvalues m,’ of n,’ and we have 
ny |m') =n |m/') 

S,+ | m/ > =(2S)27(m/' + 1)42(1 — n,'/2S)1? | n' +1) fotttees (10) 


Sp |m) = (2871 —(m'— 1)/2S}? |m'+1). J 


Creation and destruction operators are again defined by the equations 


a,*| ny) = (my +1)? [n/ +1), a,|m>=(m/)?|m/-1y ...... (11) 
to give N, =a,*a; | (12) 
Opt =(29)""a,* (1 —a,*af2S)5 0S = (25 )""1 — a 70,28) a:. i 


This analysis parallels that of Holstein and Primakoff in the ferromagnetic 
case, and with them we shall assume that factors such as (1 —a,*a,/2S)/? may 
be dropped, since their expectation value is (1—2,/2,S)"? and the expectation 
value of the spin deviation is to be assumed to be very much smaller than its 
maximum 2S. ‘This is effectively assuming a low temperature, but the effect 
of zero-point motion must not be neglected, and the whole question is discussed 
in part II of the present paper (Ziman 1952). In the same way we shall neglect 
terms in the Hamiltonian of higher order than the second in a, a,*. 

To this approximation eqns. (5), (8), (9) and (12) may be consolidated as 


follows: Si = (4S)1?(a,+ a,*) 
Op mut exp (WR a0) (ap ,2) pe eas wes (13)5 4 
SiO =exp (tw. R,)(S—a,*a,), 
and putting these into the Hamiltonian (4) we have 
KE = — Ly mJ (ES) {(Qqm + UG" gy * + YAn,* + 27" Ay, 
— Gyan LA + Ay Ag, *® — AA, — Ay* Ay + 24;* A, + 2A, *Am)} 
+ epi yi exp ((wR)iagAatO wee (14) 
where C= —2) J mom? — 1 E8S{H, +H exp(tw.R,)} ...... (15) 


and Cum = Xp (zw. R,,,,). 
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The term in C which depends on H, will vanish; it represents the magnetic 
moment of the basic Ising array, which must be zero. We need no longer 
include C in the Hamiltonian, since it is merely a constant energy. 


$4. SPIN-WAVE REPRESENTATION 
In order to reduce the double summations in (14) we perform a Fourier | 
transformation on the operators. Let 


a, = N-N2¥).a, exp (7K,. R,), . a,= N-?X,a, exp (—7zK,. R,), (16) 
6 = IN Ad, exp (21K, 6 Ry a,* = N-42X,0,* exp (2K,. R;). oe 


K, is a wave vector, and the periodicity conditions require that it shall be of 


the form DY 2. ett 
K, =27 (a + Get zvs) sat los (17) 


where G,, G,, G, are the dimensions of the crystal measured in lattice spacings 
along the basis vectors and X,, As, Az are integers such that —}G, <A, <$G,—1, 
— 4G, <A, <4G,-1, —43G, <A, <4G,-1. 
To express the Hamiltonian in terms of these new operators we must 
evaluate terms such as 
Ue f(Rim)Qq"@m=N* XL  f(Rim)oy*a, exp {2(K,.R,—K,,.R,,)}. 
l,m l,m, A, 2 
Putting R, =R,,, we obtain 
ee = N4f(R,,)ay*a, exp {a(K,—K,).R,} exp dee -R,}- 
Now the summations over / and h may be carried out independently so long as 
effects at the surface of the specimen may be neglected. In an antiferromagnet 
this is always allowable, since there is no intrinsic magnetization leading to 
shape-dependent effects. The summation over / merely gives N8,,, and the 
expression becomes 


UAUpf(R,) exp (CK, Ry) }a,*a;. 
For terms containing the factor exp (‘w.R,) the result is somewhat different. 
has, 


x f( Rin )Qy* Qn exp (7w s R,) = py NTF (Rin )&x ay 


lm lm, Ayn 


exp {u(K,.R,—K,.R,, +w.R,)} 


m 


= LY N*f(Ry)ox*a, 


Ly A, 2 
exp (7(K,—K,+w).R,} exp (7K,.R,,) 
=U i2pf(R,)a,*a, exp ((K,.R,,)} 
where A’ is obtained from A by the equation 
KK. =Kjkw i es (18) 
(the sign being chosen to bring K, into the region defined for K,). 


Terms in aa,, lead to terms in «,«%_;, as can readily be verified, and by 
carrying through all the summations and transformations we finally obtain 


KH = 2),(Aj0;* "+ $Bia,0_,+4B,* XK” *o_,* + 98H «,* ") eehetiel sk (19) 
where A,=gBH,+2%,2SJ,{o,—k(1+¢,) cos(K,.R,)} +... (20) 
B,=—%,Sd ,(1—¢,) exp (7K Re ee (21) 


B,* is complex conjugate of B,; J, =J(R,); o,=0(R,) =exp (iw. R,). 


Antiferromagnetism by the Spin Wave Method: I 545 


In order to write this Hamiltonian without any concealed repetition of 
terms, the K, space is dissected into four distinct volumes, such that when K, 
lies in one, Ky, K_,, K_, lie one in each of the others. Thus, in the standard 
example of the simple cubic lattice, we could take the regions defined by 


0<K@<h; —dn7 <K_<0; —27<K,%<—47; dp <K_,) <n, 


The summation is now to be taken over the first region only, the operators 
being written with their appropriate labels to cover the other three regions. 
That is to say, 


H =X {A ja, %0,+ A_jx_j%a_,+ Ayoy ay + A_yo_ yoy 
+ Bix ,%_2 + Bete oa + B yh yh_y + Beh hag 
+ gBH (a, *% oy +o yay tay a, tay *a_y)}. cn eees (22) 


This expression contains no double sums, and no cross-products except those 
involving four operators all belonging to a definite value of A. Each such set may 
be diagonalized by a procedure due to Polder (1949), which leaves the 
commutation properties unaltered. We find that 


He =E,fu;*(Ce*Q)VCQ)+4) +42 (C#OVC_Q) +9) — (Apt Ay} 2-23) 
where pe CA? | Balt) eRe ets... (24) 
and EPG ECHNOMMSSE et (25) 


The summation is here over the half space 0 <K,<z. 

These equations have a simple physical significance. ‘The operator 
C..*(A)C(A) has positive integral eigenvalues, N,*=0,1,2,.... The 
corresponding eigenfunctions are standing waves of spin deviations, with 
wave vector K,, propagated throughout the lattice; each spin wave quantum 
contributes an energy p,~. The two signs in (24) show that the spin waves fall 
into pairs, whose energy levels are altered in opposite directions by the external 
magnetic field, according to whether the spin deviations they cause are aided 
or opposed by this field. 

In this spin-wave representation the role of the anisotropy field is clear. 
If both H, and Hy, were zero, reference to (20) and (21) would show that for 
certain values of A (in particular for K,=0, +w) the spin waves have zero energy. 
This would mean that any number of ‘quanta’ could be put into those modes, 
turning the system out of the z direction, making the approximations invalid 
and the quantization in a particular direction meaningless. ‘The significance 
of H, is that it prevents .,~ from vanishing for any 4, and we may assert this as 
a general requirement rather than proposing a definite form for the anisotropy. 
We should also note that when H, is sufficiently large we may find that x, vanishes. 
This would correspond to the tendency for a domain to set itself perpendicular 
to the applied field, as pointed out by Hulthén (1936). But the field necessary 
for this is much larger than H,, since (A,?—|B,|”)"? is at least of the order 
(gBH,J,)!2 and J,> BH. | 


So. LHE GROUND STATE 
We have shown that the energy levels of our system may be represented by 
the formula B=¥,{(N,*+4+4)eat+(Niot Dar —(Art Av) eens (26) 
It has been pointed out that alternating Ising model basic state with which we 
started is not an eigenstate of the Hamiltonian. But if, in (26), we were to take 
PROC. PHYS. SOC. LXV, 7—A 36 
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each N,* to be zero, we should obtain the lowest energy level of the system, and 
the corresponding eigenstate would be our estimate of the lowest, or ground 
state. This differs from the basic state by the zero-point motion corresponding 
~ to a half quantum in each spin-wave mode. It would be possible to represent 
this state in terms of spin deviations, but the result would not have very great 
significance; the spin-wave representation is the more natural description. 
What is of particular interest is the actual value of the energy of the ground 
state. Using the parallel ordered state (energy — NX,,J,,S?) as a reference level, 
and remembering the term C dropped from (14), we have (in zero external field) 


Ey =N2XpJ yS*(1—o7) +32 a{(Mat ter )—2WAs+ Az} wrens (27) 


This may be evaluated by replacing the summation by an integral in any particular 
case. It has been shown by Anderson (1951 b) that for a system with nearest 
neighbour interactions J between two distinct sub-lattices 


E 2O7INSl4 7075) = Se ee (28) 


where Z is the number of nearest neighbours and y lies between 0 and 1. 
From (27) we have calculated the values of y shown in the table. 


Lattice type Zi ¥ 
Linear chain 2 0-726 
Square plane + 0-612 
Simple cubic 6 0-57 
Body-centred cubic 8 0-59 


(The anisotropy has been neglected in this calculation.) 


These results agree entirely with those obtained by Anderson (1951a) using 
the Kramers—Heller technique. It is encouraging to note that they lie near the 
centre of the range of variation allowed by (28) and that the estimated ground 
state energy for the linear chain is close to the value corresponding to y =0-7726 
found by Hulthén (1938) by an exact calculation for spin 4. The difference 
between E, and the Ising value is quite small for large co-ordination number; 
for the simple cubic lattice it is about 9-5°%. This suggests that the basic state 
is a good starting point for an approximate method, in agreement with the 
neutron diffraction experiments (Shull, Strauser and Wollan 1951) which have 
demonstrated the reality of the alternating lattice picture. 


APPENDIX 


These calculations may be generalized to include more complex forms of 
interaction. ‘Thus, in place of (4), we may take the following Hamiltonian: 


H = — Fy m2 mS s Sin ae Lg PS{H, ap Ay exp (tw $ R;)} 
os $y mp Rig °{R *S)- G’ * S., - 38, 2 G : RinSm e G V Rin} Bernese toe (29) 


where G is the tensor of the Landé g factors for non-isotropic spins. This is 
the most general form for the inclusion of dipole-dipole interactions. In order 
to make certain sums vanish, we must insist that G shall have its principal axes 
along the crystallographic axes v, y, z, but this is not a very stringent restriction. 
Putting 


lm 


G =g ali + oe os ieee (30) 
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and making the same substitutions and transformations as before, we get, 
He =X (Aja, *0, + 3B x j0_,+ $B, a, *a_y* + Cpa *ay +4D 0,0 y 
ED Saher bigs pak OM Sleep. (31) 
where 
A, =8,BH, + X,[{2SJ,— SPR, *g,"(R,? — 3Z,7)}o,, — SJ,(1+0,) cos (K, - R;) 
+3SP°R, > {ga(R,?—3X;7) +2,70,(R,?—3Y,7)} cos(K,.R,)] .(32) 
B,=%,[—SJ,(1—o,) + FSP Ry {8 2 (Ri? — 3X 77) — 8,70 (Ri —-3Y,7)}] 


SCeLO KG Ry ne eee (33) 
C,=82PH, +12, 3SPPR yg a8 —on)XnY, exp (1K,-Ry) vee (34) 
Eth peop Re eo lo A KY 1 exp (@K,.R,) So ven nes (35) 


(X,, Y;,, Z, are the components of R, along the crystallographic axes—not 
necessarily the generating vectors of the lattice). 

The same type of reduction process leads to an equation identical with (23), 
except that ,~ is now defined by 


(u,*)? =31(A?— | B,P)+(AxP?— |ByP)+2|C,?-2D? 
+ {(A?— |B, P-A,?+ | By P)P?+4E (4, + 47? —[B,— BrP)? 
+4(F,[4,+ 4,-B,+B,]-1D,[A,-A,—B,+B,)) 
x (F[4,+A,+B,+B,]+1D,,[4,-A,+B,-B,])P?] oe. . (36) 
where C,=£,+7F,. 

It is not obvious from these equations whether j.,~ will always vanish for some 
values of A when H, =0, and it may be possible to interpret the anisotropy energy 
in terms of the dipole-dipole interaction, but this is a topic which will not be 
pursued here. 
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Antiferromagnetism by the Spin Wave Method 
II: Magnetic Properties 
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ABSTRACT. From the energy levels obtained in a previous paper by the author the 
partition function and magnetic susceptibilities are calculated. x, behaves like aT? at 
low temperatures, whilst y, is constant. 'The magnetic resonance absorption frequencies 
are also calculated and shown to depend strongly on the anisotropy field. By representing 
the true Hamiltonian in terms of the operators of the approximate Hamiltonian, correction 
is made for the disordering of the system by the spin waves, and susceptibility curves 
agreeing well with experiment can be calculated. "The method is applicable to different 
types of lattice and of interaction, is not confined to zero fields, and is valid up to near the 
critical temperature. 


$1. THE PARALLEL SUSCEPTIBILITY 
N the previous paper of this series (Ziman 1952, to be referred to as I) we have 
| shown that the energy levels for an antiferromagnet near the ground state 
may be expressed approximately in the form 


B=2,(Nat+e)eatt+(N7 +d)pa hw wee (1) 
where Ha (Al = | Bp) eee. 6 Ayeyett na ee (2) 
and the N,* are integers equal to or greater than 0. (Constant terms have been 


dropped.) 
The partition function may thus be written : 


(P.F.) =e, 2nexp(—E/RT) 


Nyt, Na : 
= 2X [exp{—2%,(N,*++4)u,t/RT} exp {—2,(N7-+4))u,/RT}) 
Ns Ny 
=I1,{4 cosech (55/2 R)icosechi(a,5/2k2 \\a) on 2 eee (3) 
Hence the magnetization per unit volume is given by 
kT oO (de) 
Me, aH, In (P.F.)= — 57%, {coth (u,+/2kT)—coth (u;-/2RT)}, ...... (4) 


where V is the volume of the crystal, and y,, the differential susceptibility, is 
obtained from 
oM gp 
Xi oH, ~ AVRT 
In the limit, when the crystal is large, we may replace the sum over A by an 
integral through the volume spanned by the wave vector. Thus, since 
K , = 277{(Ay/Gy)vy + (Ag/G2)vo + (As/G3)¥3}, where A,, Ag, As, are integers, we have 


my gol 
Xi DV RT 6473 


2, {cosech? (u,+/2RT) + cosech? (u,-/2RT)}. ...... (5) 


[|| fcosech® (u,+/2kT) + cosech®(;-/2RT)}dK;, ...... (6) 


—m 
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where N=G,G,G; and we have extended the integration over the whole space 
of K, and divided by 2, since the summation was over only half this volume. 

This formula is quite general. It represents what is usually known as the 
parallel susceptibility—the susceptibility for a magnetic field applied parallel 
to the domain axis. The only difficulty is to evaluate the triple integral; in 
simple cases this can be done numerically. In the absence of anisotropy energy, 
or in very strong magnetic fields, the integrand will have singularities where 
u,~ =0 and the integral may not converge. This is a convenient criterion for the 
stability of the system. Thus, if x, does not converge even in a vanishingly 
small field, it suggests that the system is not stable against thermal agitation, and 
antiferromagnetism will not appear. For example, it is easy to see that the 
linear chain, and any plane lattice with only short-range interaction, should not 
be antiferromagnetic (cf. Yin-Yuan Li 1951). 


2-0 35. 30 


I'5 
RkT/28\J| 


Fig. 1. Parallel susceptibility (uncorrected) for simple cubic lattice. 


- For the simple cubic lattice, with nearest neighbour interaction J and 
negligible anisotropy field, we have 


pat =4S|J|{9—(cos K,% + cos K,@ + cos K,®)?}2 + 98H, ..... a(2) 


and x, convergesif H,=0. Sucha lattice is thus spontaneously antiferromagnetic ; 
for fig. 1 we have computed y,/xo as a function of RT/2S|J|, where 


Xo = Ng?B?/24VS|J|. eee. (8) 


At low temperatures x, behaves like (RT/2S|J|)?xo/61/3, but becomes 
approximately linear in 7 as the temperature increases. No singularity 
corresponding to a Curie point appears; this is not surprising since the analysis 
is only valid at low temperatures. If we had included a small anisotropy term in 
this calculation the susceptibility would be only slightly reduced; its effect is 
only appreciable near poles of the integrand, and these contribute important 
terms only at very low temperatures. 
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§2, THE PERPENDICULAR SUSCEPTIBILITY | 
We now consider the effect of an external field H, in the x direction, — 
i.e. perpendicular to the domain axis. This contributes a term 
GO == ep USOT Th YE A eae (9) 
to the Hamiltonian, and we shall treat this as a perturbation on our main 
Hamiltonian (for simplicity, in the absence of H,). From eqns. (13) and (16) 
in I we have 
KH” = — HygB(ZS)"? X, (a, + a,*) 
= —H,,gB(4S/N)1? x, , {exp (—iK,.R,)a,+exp (4K,.R,)x,*} 
HT OG RING Aer ogy) a SOURS ees (10) 
where a, is the operator «, for K,=0. 
This commutes with the whole of # except . 
He 9 SA gto to he Bol holo t Corton ea 0 eee (11) 
As shown by Holstein and Primakoff (1940), this is reduced to diagonal form by | 
the transformation 


p= 1 CotiCe*, Oy = LCo* Cs; 
where 
nog eG eer cf VA (Ag = Sot | 

b= 15 (4,2—B ey? and J,=— 2 (A?—B AP soba (12) 
(This choice of signs, not clearly indicated by Holstein and Primakoff, is 
essential to the reduction—which is only a special case of the method used 
previously in I.) 

We may write (A)?— B,”)" =, which is small, but not zero if our domain 
is stable. Thus 


HK’ = — H,gB(ZNS)"(Co + Co*){([Ao + Hol) 240)? — (Ao — Ho] /2H40)"?} 

= —H,gB3(NSpo/Ap)"(Co+Co*) nae ees (13) 

where we have neglected (u/A,)? compared with unity. 
The part #, of the Hamiltonian now becomes 
KH = bo(Co*Cy + 4). ceeeee (14) 

If the eigenfunctions of C)*C, are | No), then 

Co|No)=(No)"? | No—1), Co*| Nod =(No+1)!?| No+1) 
and clearly «N)| #’|N,)=0 for all No. 


But the second-order perturbation does not vanish, and the state |.Nj) is 
shifted in energy by the amount 


(Nol 4’ |.No') (No | #’| No» 
AE= ~{ — Ss 
No +N, bo( No — No’) 
- ET ,7g?B?N Sto y (No| Cot Co*| No’) <No’ | Co+Co*| No) 
4A oto No #N, Ny — No. 
= — NSg?8?H,?/4Ap. 
Thus the application of a perpendicular field alters the energy of the system 
by AE, regardless of the state of occupation of the levels, and we have 
ees af 02 2 N 2628 
Ki Van ae 2A, Oe, CO ita Stee ee 


..(15) 
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In the case of the simple cubic lattice 


vain. een Meee teers (17) 
Hulthén (1936), using the Kramers—Heller technique, found 
X, =Sxo{l —(RT/2SJ)2/3 1/3}, spuleueren (18) 


which differs from (17) by the temperature-dependent term. The discrepancy 
presumably lies in the fact that the method used by Hulthén is a quantization 
of a semi-classical model, and his operators é and 7 are harmonic conjugates only 
when the system is unperturbed; when the state is partially disordered by 
temperature agitation their commutation properties do not remain constant. 
It will be shown later that corrections for the approximations in our Hamiltonian 
do not alter (16). It is interesting to observe that the anisotropy energy does 
not appear in (16)—or only as a small correction in (13). Yet it must always 
be included as a convergence factor, even if negligibly small. | 


§3. MAGNETIC RESONANCE 


Suppose that we have a crystal in which there is a constant magnetic field H, 
along the domain axis, and an oscillating field H, exp (ivt) perpendicular to this. 
We have shown in (10) that the perturbation due to this is of the form 


HM = ~—H, exp (tvt)gB(gNS)(aotiay*); cee wee (19) 
and this can only induce transitions between the levels due to #)+ gp Hu *a9, 
where «=a, for K,=w. But the spacing of these levels is uniform, and given 
by 4»+g8H,. ‘Thus resonance absorption can only take place at frequencies 


hya=(Ay Be) oe ee (20) 
This result is interesting in that it depends strongly on the anisotropy field. 
Reference to I (20) and (21) shows that 


hv=ep1/{(H, +2Hy)Hy} + ep ee (21) 
where g8H, = —,SJ,(1—¢,), and plays the role of an effective internal field 
due to exchange. An identical result has been obtained by Kittel (1951), using 
a sub-lattice internal field model, and it should be capable of experimental 
verification. 

If the constant field does not happen to be parallel to the domain axis, only 
its component in that direction is effective in splitting the levels, whilst only 
the component of the alternating field in the xy plane is capable of inducing 
transitions. Experiments on powders, or on crystals in which the domain axes 
are randomly oriented, will thus show very broad lines; the two components 
of (21) may not be distinct (Nagamiya 1951). 


§4. CORRECTION FOR SPIN-WAVE SATURATION 

Up to now certain terms in the Hamiltonian have been ignored. ‘These 
_are terms introduced by factors such as (1—a,*a,/2S)"", which were dropped 
from I (13) and (14). So long as the average spin deviation is small compared 
with 2S this is perhaps permissible, but we must now try to correct for the 
approximation. 

The neglected factor has a simple physical significance; it represents the 
fact that a given spin cannot acquire more than 2S deviations; it corresponds 
to interaction between the spin waves, which have hitherto been considered 


552 F. M. Ziman 


independent. We have assumed that any number of spin waves may be put — 
into the system, whereas, in fact, beyond a certain point the lattice becomes _ 
saturated with them, and an additional quantum does not so much increase the _ 
energy but redistribute it from point to point. In the summations over V,* in _ 
the partition function we have allowed an infinity of values, where in fact we 
know that there cannot be more than N*5*? states for the system as a whole. 

Referring to eqns. (4), (5), (8) and (12) of I, we can see that the complete 
Hamiltonian should be written 


12 ek a TN ie 
Mies Pali es k ~1m){(1= § a ) (1 = +) dd 
+ AyAy,* (1 - oe) (1 — ‘ost) + (1+ oy) {(1 ~ ve) ape 
SB $8)" 
— Tym {aa (1 _ oe] + Qn * en (1 _ ost 


+328 (As: Hyoy}a"a, Cs 0 Oe 00 ie ee ee eee (22) 
We now propose to replace 1—a,*a,/2S by $6=1-— (a,*a,),y/2S where 
(a,*a,)y iS the expectation value of a,*a,, averaged over all the sites, in a given 
spin-wave state. We ignore the non-commutation of a, and (1—a,*a,/2S)", etc., 
and also the correlation factor which makes (a,*a,a,,*dm),y differ from 


(4;*Q,)ay<4m* Am Yay When I and m are close together. (This is justified in the 
- Appendix.) Clearly 


HA * ~$Hy + Ygb{H, +H, exp (tw.R,)}a,*a, 2200 (23) 
where ¢ is a number which depends on the state of perapetog of the 


eigenstates |V,*) of the approximate Hamiltonian #%. 
We calculate ¢ as follows. The average of a,*a, over all sites is defined by 


(4,*@))ay = N- d,a,*a, 


=N-? & exp {i(K,—K,). R,}a,*a,. 
laeAsrA 


= N-1 Digh,* a. ececee (24) 
To find the expectation value of this operator we must transform to a 
representation in terms of C(A). In the general case, when H, does not vanish, 
this may be rather complicated, but it is sufficient to ignore this term and 
consider the system in the absence of a field. The transformation is then of 
the form used by Holstein and Primakoff (e.g. eqn. (12)), and we have 
(4;*@y)ay = N~ Di {h(C,*C, + C¥ ,C_,) +(C,C,* + C_,C*,) 
+h1(C"C)"+ CC —CCe 02 Co) eee (25) 
where the summation is over the half space of K, and J, and J, are defined as 
in (12), but with A, in place of Ab, etc. 
The terms in C,*C,*, etc., do not contribute to the expectation value 
because (N,|C,*C,*|N,)=0, and we have 
(Na | (ar*4)ay | Na) =N-* 2,4 Ny|L2C,*C, + b?C,C,* | N;) 
=N*2,((N,+ 3) Aj/Ha-3} weeeee (26) 


where now the sum is over all A, and u,=(A?— Peay ocee 
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We have seen that the anisotropy field is important chiefly as a convergence 
factor; its magnitude is small, and does not enter explicitly in most cases. 
We shall therefore, for simplicity, assume that it may also be multiplied by ¢, 
and the corrected formula for the energy levels becomes 


E* =[2(N,t+ Ni +1)u,)[1 -(2NS)72,{(N,++ Ne + 1)Ag/ea— 13] 
+2,(N,+-—N,)gBH,, ....-. (27) 


where the sum is again over the half space of K,. 

Use of this result in place of our previous equation (1) should lead to better 
estimates of the properties of the system, and there seems no reason why it should 
not apply over a wide range of temperature, up to near the Curie point. The total 
number of levels is properly restricted by the limit on ¢. We may note that 
the domain stability condition, 1,>0, is also required in order that the sum 
in (26) should be finite, although when the sum is replaced by an integral this 
may not be significant numerically. 


§5. THE GROUND STATE 

As before, this is the state in which each N,* is zero. We shall be left only 
with the zero-point motion, but this will be effective in making ¢ differ from 
unity. Indeed, for the linear chain, the integral corresponding to (26) does not 
converge, and the correction is indeterminate; the agreement of the uncorrected 
result with that of Hulthén (1938), previously noted, must be fortuitous. But 
in the other cases we need only multiply y by the zero-point value 45 of ¢ to 
correct the table given in I. The correction depends on the actual value of S, 


and may be written 
y =doy=(1—m/2S)y scene (28) 


where 2=0-195, 0-075, 0-078 for the plane quadratic, simple cubic and 
body-centred cubic lattices respectively. 

We may also obtain an expression for the average value of [S{°|. The 
results are identical with the formulae deduced by Anderson (1951). 


§6. CORRECTION TO x, 


Reference to (9) and (10), and comparison with I (8), etc., will show that we 
have dropped a factor \/¢ in #’. ‘There should also be a factor ¢ in the 
perturbed Hamiltonian #, and since #’ is squared in calculating the 
second-order perturbation, these factors exactly cancel out. ‘This is the case 
whatever the state of occupation of the levels, i.e. whatever the temperature. 
Thus (16) should be true at all temperatures—or at least until the parallel order 
has entirely broken down. Even taking into account the correlation effects, 


this is still true. 
SiCORRE CARLO NE LO Xu 


Although (27) is probably a good description of the distribution of the 
energy levels, it is impossible to put this into an expression for the partition 
function and perform the summations exactly. We must approximate, as 
follows. In ¢, instead of N, we write the thermodynamic average value of N, 
at the given temperature and in zero field. ‘Thus, in the partition function, 
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¢ is a constant for the summations over N,, and the average value of N, at a given 
temperature is obtained from 


(Me hese > S Sy, exp — {2,(Ny-+F)pug/kT) =F coth (fy,/2kT). .(29) 


This is the value of N,+4 to be put into (26), so that we obtain the implicit 
equation 


$=1—(2SN)-1X,h{coth (¢u,/2RT)A,/u,-1} —.... (30) 


from which ¢ is obtained as a function of JT. (This device for approximating 
to the partition function is not altogether satisfactory when it comes to 
calculating thermal properties, because the free energy and internal energy 
obtained by the usual thermodynamic formulae do not quite obey the 
Gibbs-Helmholtz relation, but the difficulty only arises when we wish to 
differentiate. with respect to 7.) 

The determination of the parallel susceptibility can now be done without 
difficulty, and we obtain, in place of (6), 


N 1 gee bus+gBH. 
Se Ns 2 Z 
Xu V 6473 RT | | J {cosech ( 2kT ) 
bta—&PH, 
2 
+-cosech? ( ORT aK yea rioes (31) 


It is usual to assume that the total susceptibility y of a substance containing 
randomly oriented domains is given by y=4y,*+3x, (Van Vleck 1941), and 
in fig. 2 we have plotted x/Sy. as a function of RT/2S|J| for the simple cubic 


Calculated by Spin Wave Method 
Experimental Curve (arbitrary temperature scale) 
High Temperature Theory (Yin-Yuan Li) 


“0 05 


10 1:5 2-0 2:5 
kT7/28\.7 | 


Fig. 2, Temperature variation of susceptibility for simple cubic lattice. 


lattice with various values of S. For comparison a typical experimental curve 
(Van Vleck 1941) is shown. ‘The comparison cannot be made precisely since 
we have no means of deciding the exact value of the critical temperature by 
our method. The following points may be noted: 


(a) The shape of the theoretical curves agrees with experiment over a wide 
range of temperature. 


‘ieee RT \2 
b) A A oe | 
(6) As T->0, Xi > 7g losin) Xo: 
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(c) The susceptibility rises rapidly as the temperature approaches a definite 
limit, which may be taken as the upper limit for an estimate of the Curie 
temperature. Thus, when the spin is 4, this limiting temperature is given by 
RT/2S|J|=1-7. This lies rather below the estimates of other theories, the 
lowest of which is that of Yin-Yuan Li (1951), who found RT,/2S|J|=2:-0. 
He also found that the corresponding susceptibility y, is given by y,=0-871 Sx, 
and we have plotted in fig. 2 the variation he found for the susceptibility 
above that point. It is possible that his method, which is an accurate 
quantum-mechanical treatment of the cluster of atoms surrounding a given 
site, gives very good values for the susceptibility, which depends chiefly on the 
short-range order above the Curie point, but that it is not accurate in fixing 
the critical temperature itself, since this depends on cooperative phenomena, 
involving the whole lattice. 

(d) It is suggested that if his calculations could be extended below this point, 
ignoring the possibility of the existence of an ordered state, they would define 
a curve which would intersect our curve at a point very close to the true critical 
temperature. Indeed, a rough graphical extrapolation of his published curve 
suggests that it should cut at a value of y/Sy somewhat below unity; this would 
help to explain the observed ratio of y(T =0) to y,, which is always larger than the 
value 2/3 predicted by internal field theories. ‘The reduction in the value of T, 
(which appears to lie near 1-2 x 2.S|J|/R) may also explain the large experimental 
ratios of 6/T., where @ is the constant in the approximate high-temperature 
Curie-Weiss law, x = y’/(T' +8) (see, for example, Anderson 1950). The precise 
values of these ratios probably depend, however, on the spins, range of 
interaction, and lattice arrangement of the atoms. 


§8. CONCLUSION 

Throughout this work we have not attempted to discuss any particular type 
of system in detail, except the simple cubic lattice, which has served as a 
convenient model and example. We have tried, rather, to construct a theoretical 
framework of sufficient generality to cover a wide variety of possibilities. 
Application of the method to more complex lattices is at present under 
investigation, but the following summarizes the range of application of the 
theory: (a) Any type of lattice with one magnetic ion per unit cell may be studied. 
(6) Any type of interaction Hamiltonian may be used, and the range of the 
interaction forces need not be restricted to a single lattice spacing (although for 
long-range forces, involving many atoms, the computation may become too 
laborious). (c) The method is not restricted to vanishingly small external fields, 
as long as due account is taken of anisotropy energy in maintaining stable 
domains. (d) The calculations are reliable up to temperatures approaching the 
critical temperature, and because we have used a representation in terms of 
operators with exact commutation relations it is possible to make further 
corrections to the energy levels by simple perturbation methods. 
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APPENDIX 


We wish to calculate the expectation value of (4,*a;4m*m)ay) Where J and m 
are neighbours and the average is over all sites in the lattice. Thus 


(2;* {An * An ay = N+ Day* 07, 1415 
where the (/+1)th site is the neighbour of the /th, and distant vector a from it. 
In terms of the Fourier transform operators 


(@,*a;Am " Ora)ay, 
= N73 Dayo, *%,0,*%, exp 7{R,. K,—R,. K, +(R, +a). K, — (R,+a).K,} 
= NN Daypota*e,t,*%, exp {ia.(K,—K,)}d(A—p+v—o). 


This operator is to be evaluated in a given state |N,). From the properties 
of the operators the expectation value vanishes unless A=, v=o, or A= —», 
fp=—a, or A=o, w=v. But consider, for example, this last case: 


(Ny N,,| N-? Dye ;0,,*0, exp {ia.(K,—K,)}|N,, N,) 


2 
= | y3{ (N+) 52 — jh exp a. Ks 
In the state in which N,+ 4 is given by the average (28) this is zero, since there ~ 
are always pairs of terms for which A, and mw, are of the same but exp(za.K,) 
of opposite sign. Since the correlation factor is only a small correction to the 
Hamiltonian in any case, the most probable distribution of occupation numbers 
will be that of (28), and there will be no contribution from these terms. The 
only term which does not thus vanish is that in which A=, v=o, and we have, 
immediately, 


(Gy AA Ay, Day = N~? Ey (0303) (Oey * Oy) = (KG,* Gy ))? 
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An Anticoincidence Experiment on Cosmic Rays at a 
Depth of 50 Metres Water Equivalent 


The existence of a neutral component of cosmic rays underground has been suggested 
several times in the past but always rejected. Quite recently it has been reviewed again, ina 
new form, by Hayakawa (1952) with the aim of explaining the pairs of penetrating particles, 
which, according to various authors (Braddick et al. 1951, George and Trent 1949) are 


produced underground by the incident charged particles with a cross section Cexp Of about 
5 x 10-?° em?/nucleon. 
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We have recently investigated this problem and obtained results which seem to indicate 
that the cross section of the phenomena described above is appreciably smaller 
(Sexp~5 X 10-8° cm?/nucleon) than was thought before (Amaldi et al. 1952). If our 
conclusion is definitely proved it will no longer be necessary to invoke the existence 
of particles other than 4 mesons in order to explain the observed pair production. 

However, at the end of his paper Hayakawa suggests, as a test of his theory, an anti- 
coincidence experiment which happens to have been made by us during our work 
underground. We report here only on this point while a more extensive discussion of our 
results will appear in J] Nuovo Cimento. 'The figure shows our experimental set-up. We 


oe S 
eZ HK" oaecn WK" He oc" Sse om MHH& 
WIMMMEE@E eE=— |1=C 000 

B 


A 


Experimental set up of counters (not to scale). 


consider in this letter only the fourfold coincidences C, between the four sets of counters 
BCDE (all counters of each set are in parallel; between the various counter sets there is a 
total of 15 cm Pb) and the fivefold coincidences C;, counted simultaneously, between the 
four preceding sets and the counters A all in parallel. Measurements have been made with 
and without the layer of Pb & whose thickness As was either 5 or 7:5 cm. 

Assuming with Hayakawa that underground there are, besides charged particles, also 
neutral particles which can undergo a charge exchange process with cross section o9, we have 
Bae our the deed Ole NeG 4 CA 0) Nba Omen) (1) 
with the lead 2 CAs) =[Neet+NywooNAs] Gy, C3(Az)=enaNcG, —... .-. (2) 
where G, is a geometrical factor of the telescope, N, and N,o the numbers of charged and 
neutral particles incident in the vertical direction per cm? and steradian, 7,4 the efficiency of 
the set of counter A, € the absorption of & and N Avogadro’s number. 

Combining (1) with (2) we get 
Nieto 1 eek C;(As) 

C,(As)—C;(As)]— 
Ne REN Cy C,(0) [CyC 5) 5 =)] C,(0) 
Equation (3) is written in terms of the differences between the fourfold and the fivefold 


coincidence rates in order to take advantage of the fact that with our coincidence circuits the 
fivefold coincidences were not statistically independent of the fourfold coincidences. 


[C,(0)— cxon. Ri (3) 


, fe c — 
As(g/cm?)|| ¢ (min) C, | (optante lenin) Cs | (counts froin) Saks 


12°3742- 0-024 0-431 -+-0-0045 
12-:215+0-024 0-415 +0-0044 
12:019+0-026 0:420-40-0049 


12-805-++0:024 || 269-798 
12:630+0-024 || 256-728 
12-439+0-026 || 218-786 


0 21-803 || 279-188 
Sy) 21-018 |] 265-458 
85 18-203 || 226-433 


Our results are given in the table. Introducing the appropriate values in eqn. (3) we get 
for As=57 gjcem®, Nyeo/Ne=(—2°5+1°5) x 10-** cm?/nucleon, 


for As=85 gjcm®, Nyoo/No=(1°5+ 1) x 10-?° cm?/nucleon. 
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These two values are equal inside the experimental errors and therefore we can take their 
mean value which is 
aA) (Sila eS) <a cm2?/nucleon. Sette hn (4) 


Cc 
Now, according to Hayakawa (1952) 
Ne=NutNet and oo=fo,, wets ees (5) 


where N,, and Nx+ are the numbers of uw mesons and of k= mesons, cy is the cross section for 
production of pairs of penetrating particles by «+ mesons, and f a numerical factor of the 
order of a hundred (charge exchange is a first order, pair production is a third order process). 

Noting that the experimental cross section for pair production is connected with oy, 
by the obvious relation Net 


a 6 
exp N, ce Wee on (6) 

and combining (4) (5) and (6), we get 
Nxt F Fexp (7) 


Ne  (=121-5)x10-= 


Introducing the numerical values f100 and ogxyp=5 x 10~-*® or 5 107° we have 
Nxt/Neo = 150 or Nyt/Neo 2 15. 

Both values are much higher than those expected according to Hayakawa’s theory, 
which suggests a value of Ny+/N of the order of 1. 

It seems to us that we can conclude that the theory of Hayakawa does not agree at all 
with our anticoincidence experiment, whatever is the value of the cross section for pair 
production by cosmic rays underground. 


Istituto di Fisica dell’ Universita, E. AMALDI. 
Centro di Studio per la Fisica Nucleare del C.N.R., C. CasTAGNOLI. 
Roma. S. SCIurTh 
Istituto di Fisica dell Universita A. GIGrr: 
Pavia. 


15th April 1952. 
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Note on the Differential Momentum Spectrum of 
Cosmic-Ray Protons at Sea Level 


Rossi (1948) has used the data of Rochester and Bound (1940) to determine one point 
of the differential momentum spectrum of cosmic-ray protons at sea level. Due to a 
misprint in the text of Rossi’s article, several authors (e.g. Mylroi and Wilson 1951, Rochester 
and Rosser 1951, and Puppi and Dallaporta 1951) have given an incorrect value for this 
spectral point. ‘The misprint in the text of Rossi’s article states that the number of protons 
within a range interval of 2 cm of lead is 5 x 10-7 measured in units of cm-2 sec-! sterad—. 
As seen in the table accompanying the text, this figure actually is the value of the differential 
range spectrum of protons in g~'sec™!sterad1. Thus the flux of protons through 
Rochester and Bound’s apparatus is 8:6 x 10-§ cm=? sec} sterad-!. 

T'aking into account the wall thicknesses of the counters, cloud chamber and housing 
around the apparatus, all protons with momentum in the range 0:34—0-53 x 109 ev/c 
would be stopped in the 2 cm lead block. The flux of protons in this momentum range 
corresponds to a vertical intensity of J,=4°5+ 1-6 x 10-8 em-? sec™! sterad-! (mev/c)-1. 
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Because the purpose of the experiment of Rochester and Bound was to measure the 
attenuation of protons in lead rather than to determine the intensity of protons, a very 
stringent criterion was placed upon those tracks which were accepted as protons. The 
eight particles considered were probably all greater than about four times minimum 
ionization. The momentum of a proton with this ionization is only 0:40 x 10%ev/c, and 
not the upper limit of 0-53 x 10° ev/c used in the above calculation, which corresponds 
to 2:5 times minimum ionization. If now we take the momentum interval for counted 
protons to be 0-34 to 0:-40X10* ev/c instead of 0-34 to 0:53 x 10° ev/c, we obtain a 
corrected value of the spectrum J, =14+ 5 x 1078 cm~* sec sterad-! (Mev/c)-. 

The table gives a summary of the present data available in the low momentum region 
of the differential spectrum for protons. Two sets of values of the spectrum are given for 
the data of Mylroi and Wilson. The first is obtained from the differential meson spectrum 
given by Rossi (1948), and the second by assuming that Rossi’s meson spectrum is too low 
by 20°, as indicated by York (1952). It is to be noted that this improves the agreement 
between Mylroi and Wilson’s values for the proton spectrum and those of Merkle, 
Goldwasser and Brode (1950), Goldwasser and Merkle (1951), Ballam (private 
communication) and York (1952). 


Momentum interval Corrected Iy 
Observer (x 10-* ev/c) Ty (x 108) (x 108) 
Rochester and Bound (1940) 0:34-0:53 4:51+1-6 — 
0-34-0-40 — 14+45 
Ballam 0:42-0:91 30+8 —_ 
0:51-0:96 31+14 — 
Merkle, Goldwasser and Brode (1950) 0:48-0:88 2343 _ 
Goldwasser and Merkle (1951) 0-78-1-10 20+5 — 
Mylroi and Wilson (1951) 0:53-0:71 16:8-+4:4 20:1+5:°3 
0:71-0:93 8941-6 10:7+1:9 
0:93-1:35 5°5--0-9 6:6+1:1 
York (1952) 0:55-0:75 28+3 — 
0:75-0°88 10:7+0°8 


The author is indebted to Professor Rossi for pointing out the inconsistency between 
the text and table in his article and to Dr. G. D. Rochester for supplying the original data 
from the experiment of Rochester and Bound. The author would like to express his 
thanks to Dr. Rochester and Professor R. B. Brode for several stimulating conversations 
and for checking the calculations used in the table. He would also like to thank Professor 
P. M. S. Blackett for extending the facilities of his laboratory during the author’s stay in 
Manchester and the United States Educational Commission in the United Kingdom for 


financial support. 
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Paramagnetic Resonance in Ferric Rubidium Sulphate Alum 


Ferric ammonium sulphate alum, in which the ferric ion is in a °S state, has been 
extensively used in adiabatic demagnetization experiments. The sextet level is split by the 
crystalline field; if the field has cubic symmetry a doublet and a quadruplet remain, whose _ 
separation, according to specific heat measurements (Cooke 1949) is 0-14cm™!. Bagguley 
et al. (1948) found, however, that paramagnetic resonance experiments indicated a splitting, 
assuming a cubic field, of less than 0-1 cm“; they suggested that a trigonal field was present 
which further splits the quadruplet into two doublets. In two salts of manganese, where the 
Mn?** ion is isoelectronic with Fe**+, the trigonal splitting is much larger than the cubic 
splitting (Bleaney 1950, Bleaney and Ingram 1951). 

Various attempts (Bijl 1950, Wiedner et al. 1949, Ubbink et al. 1951) to clarify the 
position in the iron alums by the use of diluted salts have not been completely successful _ 
because the side lines are unduly broad, presumably because the splitting of the levels by the ~ 
crystalline field is not uniform throughout the crystal. Meyer (1951) tried to explain the 
observed spectra using a combined cubic and trigonal field similar to that found to be 
successful in the manganese salts, but some discrepancies remained, though no exact test 
of the theory was possible owing to lack of experimental resolution. 

We have found that high resolution can be obtained with diluted ferric rubidium sulphate 
alum. ‘The results show that the cubic axes of the ferric ion are not identical with those of 
the crystal, but are rotated by an angle of 7$+ 4° about the trigonal axis. ‘The splitting can 
be interpreted using the Hamiltonian for manganese (Bleaney and Ingram 1951) with an 
additional trigonal term corresponding to the spherical harmonic Y,°, whose magnitude is 
denoted by the parameter F. In a strong field making an angle 6 with the trigonal axis, 
the spectrum consists of five lines, whose positions are given by the formulae (correct only 
to the first order) 


hy=gBH + {2pa+2D(3 cos?@—1)+4F(35 cos*@—30 cos?6+3)} 
hy=gbH+ {—3pa+D(3 cos?6—1)— 3F(35 cos*#—30 cos?0+ 3)} 


hy=gpH 
The values of the constants are as follows (a, D, F in cm7) : 
Temp. (°K) g a D F 
90 2:004 + 0-003 —0-0134 +0-0022 —0-0003 
20 2-004 + 0-001 » —0-0134 +0-0031 —0-0003 


The relative signs of the constants are found from fitting the spectrum, the actual signs 
from the varying relative intensity of the lines owing to their unequal populations at low 
temperatures. 

The rotation of the cubic axes is deduced from a detailed examination of the spectrum at 
intervals of 5° in the (100) and (110) planes, and is slightly smaller than that (9$°) found by 
Lipson (1935) from x-ray analysis of aluminium rubidium alum. The excellent agreement 
with the measurements that can be obtained with the values of a, D and F given above is 
illustrated by the following table. This gives the splitting between the central line and the 
side lines at 20 °K, second-order displacements being eliminated by taking the mean of 
corresponding lines on either side. 


Direction Multiplicity Experimental Calculated 

splitting (gauss) splitting (gauss) 
(100) Quadruplet 264 262 
on 329 B27 
(111) Singlet 167 169 
; “ 320 319 
Triplet 233 235 
Pa 124 125 
(110) Doublet 123 128 
ff 123 Lunresolved 125 
46 45 43 
ae unresolved from 4+ 


central line 
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In zero field the levels form three doublets, calculated to lie at 0, 0-0026° and 0:0633°x 
respectively. The corresponding contribution to the magnetic specific heat tail is 
CT?/R=2a?+56D?/9+14F?/27—28aF/27=8:°58 x 10-4 (°k)?. 

Preliminary measurements on ferric potassium selenate alum and ferric methylamine 
sulphate alum indicate that a is about the same as in the rubidium alum, but the values of D 
are roughly 0:01 cm™ and 0-19 cm™! respectively, the sign in each case being the same as 
that of a. 
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A Note on the Theory of Conduction in Metals 


The object of the present note is to draw attention to a feature of the theory of metallic 
conduction which appears to have been overlooked. In the usual version of the theory 
(see, for example, Wilson 1936) it is assumed that electrons can be scattered by lattice 
vibrations with all wave numbers |q| <q, where gj corresponds to the cut-off frequency 
of the Debye spectrum. However, in any collision involving a (quasi-free) electron with 
wave vector k and a lattice wave q, the laws of conservation of energy and momentum 
require that |k|=|k+q]|, and thus an electron in the state k can only interact with lattice 
waves for which |q|<2{k|. This fact is well known in the theory of semiconductors, 
and has also been taken into account in Froéhlich’s theory of superconductivity (Frohlich 
1950); but its consequences for the theory of normal metals do not seem to have been 
noticed. In a metal it is only necessary to consider electrons with wave vectors |k|=K,y 
corresponding to the surface of the Fermi distribution. It follows that, if 2Ky)<q, the 
upper limit of integration in the collision operator is 2K, and not go. A number of 
interesting consequences follow. 

(a) If 2Ky<qo the usual expression for the time of relaxation at high temperatures 
(Wilson 1936, p. 208) must be increased by a factor (q9/2Ko)*. Since Ko/qo=(4Na)", 
where N, is the number of conduction electrons per atom, the condition 2Ky<qp is 
equivalent to N,<4, and the factor (qo/2Ko)* equals (4N,)-**. The high-temperature 
electrical and thermal conductivities are to be multiplied by the same factor. 

(b) The usual interpolation formulae for the temperature variation of the electrical 
and electronic thermal conductivities of metals (discussed, for example, by Makinson 1938) 
involve the parameter D/l=qo'/2K fe =244 Ng *. ‘These formulae are clearly only valid if 
D/€<2, and so do not apply to metals in which the number of conduction electrons is 
small; in particular, Makinson’s discussion of the thermal conductivity of bismuth is not 
valid. 

(c) The modified interpolation formulae which hold when D/¢>2 are easily obtained 
by the usual methods (Wilson 1937, Sondheimer 1950). In the notation used by these 
authors it is found that, for an ideally pure metal, 


1 3ahD (T\? ©’ 
+= Fela) x(F), 5 Oro toA0uc (1) 


1 - seurtna( 5) [2 ) + 2(2) ee(z) -(G)}} 


PROC. PHYS. SOC. LXV, 7—A 37 
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where @’=4/(2¢/D)®. It follows that the temperature variation of o and « is qualitatively 
unchanged; the values of the conductivities at sufficiently low temperatures are also not 
affected, as is obvious since only long lattice waves are then excited in any case. The ratio 
of the high- to the low-temperature conductivities is thus increased by the factor | 
(D/2¢)?=(4N,)—8. Further, the transition from ‘high’ to ‘low’ temperature behaviour __ 
takes place at a temperature ©’ which is less than the Debye temperature 9. 

The present note has been written mainly in order to warn against the use of the usual 
theory in cases where it does not apply. The new formulae do not, unfortunately, have 
more than qualitative significance, since the metals in which the number of conduction 
electrons is small (the most popular example of which is bismuth) are just those in which 
the simplifying assumptions of the theory will be furthest from the truth. 

It is perhaps, however, worth while considering briefly the implications of the new 


theory for the properties of bismuth. In this metal N, is of the order of 107%, so that — | 


D/€~80 and q/2Ky~6. (Makinson (1938) adopted a much higher value of Na 
apparently confusing the actual number of free electrons in a metal with the ‘ effective ’ 
number in the sense defined by Mott and Jones (1936).) The most striking result of the new 
theory is that the time of relaxation, and with it the free path, of the electrons in bismuth 
are likely to be much larger than the corresponding quantities in ‘ good’ conductors. 
A quantitative estimate of the absolute values would, of course, require a much more refined 
theory, in view of the many uncertain parameters involved and the possible effect of 
additional scattering processes (inter-band transitions, ‘ Umklapprozesse’) which have 
been neglected. It is remarkable, however, that experiments on thin films and wires of 
bismuth (Eucken and Forster 1934, Justi, Kohler and Lautz 1951) indicate that the free 
path at room temperatures is of the order of 10-4 cm, which is about 20 times larger than 
the free path in silver. If the present suggestion is correct, this is simply due to the fact 
that slow electrons are prevented by the conservation laws from interacting with any but 
the longest lattice vibrations. The fact that bismuth has, in spite of this, a small electrical 
conductivity is, of course, due to the fact that the small number of conduction electrons 
outweighs the effect of the long free path. 

The present theory further predicts that the electrical resistance of bismuth should 
vary linearly with temperature down to temperatures of the order of 10° K, which does not 
agree with observation. This indicates the grossly oversimplified nature of the present 
model. The thermal conductivity should be independent of temperature, and the Lorenz 
number should have its normal value, down to temperatures of the same order. It is, 
unfortunately, difficult to test these predictions owing to the presence in bismuth of a large 
lattice component of thermal conductivity. 

The considerations outlined above also have a bearing on the interpretation of recent 
experiments on the thermal conductivity of monovalent metals. This question is taken up 
in the following note. 
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The Thermal Conductivity of Metals at Low Temperatures 


Recent experiments on the thermal conductivity of various metals at low temperatures 
(Hulm 1950, Berman and MacDonald 1951, 1952) have made it possible to check the 
predictions of the theory (Makinson 1938, Sondheimer 1950) in much greater detail than 


hitherto. In the present note we wish to make some comments on the comparison between 
theory and experiment. : 
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The electronic thermal conductivity « at sufficiently low temperatures is given by the 
formula RP a 4 cl ESP a lepton RNR Va. be (1) 


where the first term represents the effect of the lattice vibrations and the second term that 
of impurities and other static lattice defects. In considering the theoretical value of « 
(following Hulm), it must, however, be remembered that the theoretical formulae are 
different according as D/{22. (For the significance of this parameter see Makinson 1938, 
and the preceding note, Sondheimer 1952.) In fact we have 


a=6CF5(00)/7Dk 0” DICER AR, OR Ese eats - (2 a) 

; a=3DI,(00)/27?Lk 0” La OT aay Ae OE ey Peat (2 b) 
where «,, is the thermal conductivity at high temperatures and the notation is otherwise 
the same as Makinson’s. If « is regarded as a function of D/{ it has a minimum value 
O&min= 375 (0 )/7?K —O? =37-8/K .O?, which occurs when D/f=2. The observed values of « 
are in all cases rather smaller than this minimum value; some typical examples are shown 


in the following table. (The values of © have been taken equal to those which give the best 
fit for the electrical resistance.) 


Na(@=202) Cu(® = 330) Pb(®=88) 
DOr Orrin 7-1 0-90 140 
10%a (observed) 3°8 0-23 22 


(« is measured in watts/cm/°K) 


Thus no value at all of D/¢ will fit the observations. Hulm, who assumed (with Makinson) 
that eqn. (2 a) was valid for all values of D/£, concluded that the experiments appear to 
indicate an abnormally large value of D/ ¢ (and hence a small number of conduction 
electrons) in the metals studied; we now see, however, that, whatever the explanation of 
the discrepancy, this interpretation is not admissible. 

The theoretical value of D/{ in monovalent metals is 2-1/?, and it is clear that (as has 
been emphasized by Hulm) the observed thermal conductivity at low temperatures is 
considerably greater than the theoretical value even in metals like sodium and copper, to 
which the theory is most likely to apply. ‘The non-existence in monovalent metals of the 
theoretically predicted minimum in « at higher temperatures (Berman and MacDonald 1951, 
1952) is evidently a related phenomenon. It should be noted, in this connection, that the 
theoretical formula for «(T)/« involves the two parameters T/® and D/¢, and the exact 
form of the temperature variation depends sensitively on the value which is chosen for D/€. 
The temperature variation of the electrical conductivity, on the other hand, is entirely 
independent of the value of D/ ¢, and the behaviour of the thermal conductivity thus 
provides a much more searching test of the theory than does that of the electrical 
conductivity. Further, the value ©, of © which has to be used in interpreting the electrical 
conductivity data (and which should presumably be used also for the thermal conductivity) 
is usually widely different from the expected value ©,, that which corresponds to the 
longitudinal lattice vibrations (Blackman 1951); thus for sodium ©,=202, ©,~270, 
and for copper ©,=330, ®,—520. In the present state of the theory, therefore, @, has 
to be treated in effect as an arbitrary parameter, whose precise value depends on factors 
which have been neglected in the derivation of the theoretical formulae. Now the usual 
relation giving D/¢ in terms of the number of free electrons per atom (Sondheimer 1952) 
can be expected to hold quantitatively only if the simplifying assumptions of the theory 
are valid; since they are not valid, it is reasonable to postulate that D/¢ should also be 
treated as an arbitrary parameter, whose value is of the same order as that suggested by 
the number of free electrons, but not necessarily equal to it. If, for example, we suppose 
that the effective value of D/ ¢ in monovalent metals is about twice the normal value, the 
agreement between theory and experiment is very much improved, especially if the higher 
approximations to the theoretical thermal conductivity are taken into account (Sondheimer 
1950). (Note here that the minimum in « is only present if D/ €<1-7.) As we have seen, 
however, it is not possible, with any reasonable choice of ©, to fit the experimental data at 
low temperatures completely with any one value of the parameter D/C. 

A more drastic procedure, analogous to that used in the analysis of specific heat data, 
is to regard both © and D/¢ as functions of the temperature rather than as constants. 
It is then possible to adjust the two parameters at any temperature to give the right values 
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of both o(T)/o and «(T)/K.. The values of D/¢ required at low temperatures to obtain 
agreement with experiment are, however, so widely different from the normal value that 
the procedure has only formal significance. It seems clear, therefore, that the discrepancies 
encountered in the case of the thermal conductivity, while not quite so unexpected as has 


sometimes been suggested, cannot be satisfactorily explained until a more exact theory has | 


become available in which the formula for x is modified substantially. 
I wish to thank Mr. A. H. Wilson, Mr. R. S. Leigh, Dr. R. G. Chambers and 
Dr. M. Blackman for helpful discussions of these topics. 


Department of Mathematics, E. H. SONDHEIMER. 
Imperial College, London. 
2nd May 1952. 
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The 7 Mev Level of '7C 


During a cloud-chamber study of internal pair emission from #*C we have found evidence 
for the decay by pair emission of an excited state at approximately 7 Mev. 

The source consisted of 200 mc of polonium deposited in a thin layer on a 1 cm? piece 
of platinum foil. It was covered by a sheet of beryllium (25 mg/cm?) and enclosed in an 
evacuated cell having a thin mica window in front of the beryllium. This window faced a 
similar one in the wall of a cloud chamber containing air at 60 cm Hg pressure. A magnetic 
field of 550 gauss was used for momentum determination. 

Out of 79 measured pairs 72 had a total energy of 4-4+0-5 Mev, while seven had a total 
energy of 7-0+ 0-6 Mev. 

We have examined and rejected the possibilities that the tracks of the high energy pairs 
are those of normal pairs which have undergone considerable multiple scattering, or that 
they are random associations of a positron with a Compton electron. 

It is unlikely that the high energy pairs arise from a 7 Mev y-ray of !?C or *C converted 
in the platinum backing. A number of authors have searched for such a y-ray from 
Po-Be sources and a limit to its intensity relative to that of the 4-45 mev line has been set at 
0-5% by Terrell (1950) and at 0:2% by Artemov and Vlasov (1951). An intensity ratio of 
about 1°% was reported by Pringle, Roulston and Standil (1950), but Bell and Jordan (1950) 
produced evidence for regarding as spurious their own similar observations made by the 
same method. 

Phillips, Cowie and Heydenberg (1951) have observed the emission of electron pairs of 
total energy 7 Mev from "B bombarded by protons. ‘They suggest that these might originate 
in a pair-emitting state of "Be. The pairs we have observed cannot be produced in this way 
since the energy available for excitation of §Be is less than 2:1 Mev. 

We therefore conclude that they most probably arise from the decay by pair emission of a 
7 Mev excited state of 1?C. 

If we explain the failure to observe 7 Mev quanta from !°C by assigning to the 7 Mev 
state the same spin and parity as the ground state (0, even) we must expect a partial decay by 
radiation to the 4-45 Mev state for which we find, by measurements to be published, spin 2 
and even parity. 

Observations of the neutron groups from Po-Be by Guier, Bertini and Roberts (1952) 
give 1 :8 as the ratio of the population of the 7 Mev state to that of the 4-45 Mev state. Since 
the ratio of 7 Mev to 4:45 Mev pairs is of the order 1:10 and since the pair conversion 
coefficient for the 4-45 Mev state is of order 10~-%, it may be deduced that the probability of 


pair emission from the 7 Mev state is small compared with that of radiative transition to the 
4-45 Mev state. 
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We should therefore expect the emission from Po—Be of a y-ray of energy about 2:5 Mev 
and of approximately one-eighth the intensity of the main 4:45 Mev line. This intensity 
would probably not have been detected in the experiments already mentioned or in those of 
Bame and Baggett (1951). 


The Clarendon Laboratory, G. Harris. 
Oxford. W. T. Daviss. 
22nd May 1952. 
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An Investigation of the Decay Scheme of **Sm 


Previous work on 47-hour 158Sm is summarized in Nuclear Data (1950). Although there 
is a considerable measure of agreement on the radiations emitted, two different decay 
schemes have been proposed by Muehlhause (Nuclear Data 1950) and by Hill and Shepherd 
(1950). These decay schemes are inconsistent, the total disintegration energies being of the 
order of 950 kev and 800 kev respectively. The integrated spectrum technique (Bannerman, 
Lewis and Curran 1951) was therefore employed as a direct means of testing the validity of 
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Fig. 1. The 1°*Sm integrated spectrum. Fig. 2. The y-ray spectrum from Sm. 


Inset. ‘The proposed decay scheme of °8Sm. 


either scheme. A thin source (~5 g/cm?) was clamped between the two halves of a 
1 cm® crystal of Nal(T1), thus obtaining practically 100°% absorption for all radiations. 
The crystal system was mounted on an EMI 5311 photomultiplier tube and the whole 
surrounded by 2 in. of lead shielding. The end-point of the integrated spectrum occurred 
at 800+ 20 kev, and therefore agrees with the decay scheme of Hill and Shepherd. All 
features of the low energy region of the spectrum (fig. 1) decayed with the 47-hour half-life 
of %3Sm. The main peak at 100 kev further confirms the scheme of Hill and Shepherd, but 
the other features indicate a more complex final solution. Further work on the source can 
be summarized as follows. 
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The y-ray spectrum. The spectrum observed in a 0:5 cm® crystal of NaI(T1) (fig. 2) | 
showed two lines at 101+2 kev and 42+0-5 kev. This second line is the EuK x-ray 
(energy 41-5 kev) arising from internal conversion. No y-rays were found above an energy — 
of 101 kev, and an upper limit of 0-005°% can be given for the intensity of the 610 kev y-ray. | 
An asymmetry on the low-energy side of the 101 kev peak was shown to be due to the | 
emission of a few unconverted 69:5 kev y-rays (see fig. 2). The relative intensities of the 
y-ray and x-ray lines at 41:5, 69-5 and 101 kev were found to be 70%, ~4% and 26% 
respectively. 

y-y coincidence measurements. The spectrum of quanta in coincidence with 101 kev 
+-rays consisted of a single x-ray peak, indicating the presence of a highly converted y-ray in 
cascade with the 101 kev transition. It was calculated from the geometry that one in every 
three 101 kev y-quanta was accompanied by an x-ray in coincidence. Thus about 70% of 
the total K x-ray intensity is due to internal conversion of the 101 kev y-ray, the remaining 
30% arising from the 69-5 kev transition. Allowing for the fluorescence yield in the 
K shell of Eu and for L conversion (K/L ~5/1, Hill and Shepherd 1950), a value of the order 
of 3 is then obtained for the internal conversion coefficient of the 101 kev y-ray. No separate 
peak due to 69-5 kev y-rays was observed in the coincidence spectrum, but an upper limit 
of one-tenth the K x-ray peak can be given for its intensity. The internal conversion 
coefficient is therefore greater than 10. 

The B-ray spectrum. A thin source (~5 g/cm?) and a 0-5 cm$ crystal of NaI(Tl) were 
used to examine the f-ray spectrum. ‘The solid angle was kept small to reduce coincidences 
from the source, and the end-point was confirmed at 800+ 20 kev. The highest energy 
8-ray therefore goes directly to the ground state of #°°Eu. ) 

b-y coincidence measurements. Under suitable conditions of geometry, coincidences 
were observed between f-particles and 101 kev y-rays. The f-spectrum obtained when 
isolated by these coincidences had an end-point at 700+ 20 kev. This indicates clearly 
that the 101 kev y-ray transition takes place to the ground state of Eu. A third f-ray of 
energy 630 kev must also exist leading to the 69:5 kev transition, but no attempt was made to 
isolate this. 

The resultant decay scheme is illustrated in fig. 2, inset, with approximate values for the 
relative intensities. "This scheme accounts directly for the peaks at 100 kev and 170 kev in 
fig. 1, since these will be the ‘ zeros’ of the two displaced 6-spectra. The other features of 
the curve can be explained by the escape of radiations due to an imperfect 47 geometry. 
The peak at approximately 210 kev, however, cannot be accounted for on this scheme, and 
may indicate the presence of a further excited level in ®*Eu. 

The writer is indebted to Dr. S. C. Curran for many helpful suggestions and discussions 
throughout the course of the work. 


Department of Natural Philsophy, R. C. BANNERMAN. 
University of Glasgow. 
8th May 1952. 
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REVIEWS OF BOOKS 


Integral Transforms in Mathematical Physics, by C. J. TRANTER. (Monographs 
on Physical Subjects). Pp. ix+118. (London: Methuen, 1951.) 6s. 


The principal aim of this book is to outline the use of integral transforms in the solution 
of partial differential equations. The solution of an ordinary differential equation can 
often be reduced to that of an algebraic equation by the use of an integral transform. In a 
similar manner a partial differential equation in m independent variables can be reduced 
to one containing (n—1) variables. The method facilitates the solution of such equations 
with given boundary and initial conditions. 
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The following integral transforms are treated: Laplace, Fourier sine and cosine, 
complex Fourier, Hankel and Mellin. The corresponding inversion formulae are obtained 
formally. Separate chapters are then devoted to applications of these transforms, each 
being illustrated by the solution of problems—many of which are taken from recent papers 
by the author. The student of electric circuits will perhaps be disappointed to find that 
practically every example relates either to heat conduction or to elasticity. 

Of special interest are the two chapters dealing respectively with the numerical 
evaluation of ‘inversion’ integrals, and finite Fourier and Legendre transforms. The last 
chapter shows how a partial differential equation in three variables may be reduced to one 
in two, which latter can then be solved by relaxation methods. 

Those with a limited. background of mathematics will need to acquaint themselves 
with the elements of analysis and with the evaluation of integrals using contour integration. 
Given this equipment the book will be welcomed as a compact and workmanlike exposition 
of a subject which has come into great prominence in many fields. H. H. HOPKINS. 


The Classical Theory of Fields, by L. Lanpau and E. Lirsuitz. 2nd Edition. 
Pp. ix+354. (Cambridge, Mass.: Addison-Wesley, 1951.) $7.50. 


The well-known vigour and originality of thought of the senior author of this book, 
and the mathematical virtuosity of his eminent collaborator, confer a special interest on this 
product of their didactic activity. The present volume gives us only a fragmentary view 
of their conception of a course of theoretical physics. Judging from the excellence of this 
part, one would wish to have the whole work as carefully translated and neatly produced. 
‘The plan of the course must be quite unconventional, for the part under review presents 
a rather unfamiliar assemblage of topics, which are only approximately described by the 
title. Essentially, what we are given is a discussion of the production of electromagnetic 
fields in vacuo by charged particles, and a treatment of the gravitational field from the 
point of view of general relativity. However, the electromagnetic section, which is built 
up from the start on the theory of relativity, also includes a brief, but masterly, exposition 
of geometrical optics and diffraction. It is gratifying to find a simple and accurate 
discussion of the nature and production of gravitational waves and an exposition of the 
cosmological aspects of gravitation theory distinguished by an exceptional lucidity and 
sobriety. As a matter of fact, one may pick up at random any section and one is sure to 
find a concise and elegant treatment, in which the essential points are admirably brought 
out. The book is almost entirely self-contained; the reader is only assumed to have 
some familiarity with the methods of analytical mechanics. The mathematical calculations, 
though condensed, are always sufficiently explicit to be followed without difficulty by the 
advanced student. A most valuable feature is the refreshing originality of the special 
problems selected for discussion. Although the scope of the book is kept rigidly within 
the bounds of classical physics, most of these problems have been chosen with a view to 
future use in quantum theory; others have by themselves considerable interest, such as 
the radiative effects of accelerated particles in the extreme relativistic domain. Even the 
specialist will here find little known results of recent investigations by Soviet physicists. 

Does then Landau and Lifshitz’s course at last provide us with the ideal introduction 
to theoretical physics? Here I must, alas! express my disappointment. ‘The stern beauty 
of a smooth flow of rigorous mathematical deductions is bought at a heavy price. It actually 
amounts to putting the natural development of the subject upside down. I know from 
our discussions in the good old days that Landau’s ‘violon d’Ingres’ has always been his 
belief that one could encompass the whole of physics*in a monumental ‘principle of least 
action’. So far as I can judge, this is the leading idea which he has unflinchingly carried 
out in this course, and we can see the result. A few specific examples will show better 
than much philosophical discourse the kind of incongruities which this attitude entails. 
To write down the principle of least action governing the relativistic motion of a particle 
we have just to look for a suitable Lorentz invariant, and there is no other than the interval. 
This sounds fine, but if we ask for the physical basis of the requirement of Lorentz 
invariance we must work our way back to page 1, only to strand on the statement: 
“« Experiment shows that . . . all the laws of nature are identical in all inertial systems of 
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reference.” The introduction of the electromagnetic field is much worse still. The 
additional term in the Lagrangian, containing the potential, is simply thrust upon the 
reader with the comforting words: ‘It is found to have the form...” And on this 


basis the first pair of Maxwell equations is ‘derived’! It need not be emphasized how 
misleading and didactically unsound such a procedure is. 

The debate is not new, of course. The idea of treating mechanics as a ‘deductive’ 
science goes as far back as Lagrange. Commenting on the Mécanique analytique, that 
profound mathematician Poinsot shrewdly observed that if one tried to introduce into 
that great book a ‘proof’ of the ‘principle of virtual velocities’, “l?ouvrage se trouverait 
fait deux fois : je veux dire que cette démonstration comprendrait déja toute la mécanique”’. 
It seems to me that there is no escape from this situation. One might say, with a grain of 
salt, that theoretical physics should be learned twice over. The first stage should proceed 
inductively and historically, so as to retrace the dialectical unfolding of the subject. This 
stage would indeed culminate in the establishment of great universal principles; it would 
show these principles in their true light, as a condensation of experience in abstract, 
idealized form, with all the inherent limitations of their range of validity. Only then would 
the ground be safe for a second survey in the reverse direction. The principles could then 
be used as a basis for a mainly deductive analysis of more delicate problems, in the 
treatment of which they would demonstrate their real fruitfulness. For this return journey 
the advanced student could find no better guide than our authors. But for his first 
introduction to the true spirit of theoretical physics he will still have to rely on the old 
masters, Abraham and Sommerfeld, who (albeit unwittingly) display that fine sense for 
dialectics which is so strangely lacking even in the best productions of our Russian friends. 

L. ROSENFELD. 


Dictionary of Mathematical Sciences: Vol. 1, by L. Hertanp. Pp. 235. 
(London: Hoffner, 1951.) 24s. 


This work claims to be the first modern bilingual dictionary of mathematics. The 
vocabulary is mainly that associated with pure mathematics and mathematical physics; 
but, in addition, a very large number of entries seem to be commercial terms and phrases. 
Moreover, many terms from such fields as physics and astronomy are included. 

There are some surprising entries to be found within the pages of this volume. It 
seems unnecessary to explain that ‘chi (griechischer Buchstabe)’ has the meaning ‘ chi” 
—nor, indeed, to give corresponding ‘meanings’ of the other letters in the Greek alphabet. 
Again, the now rarely used, ‘q.e.d.’ is given under ‘quod erat demonstrandum (cf. 
Q.E.D.)’ and again under ‘q.e.d. (abbr.) quod . . . etc. (Lat.) (was zu bewiesen war, which. 
was to be demonstrated) (abbr.) Q.E.D.’. This seems to be a sheer waste of space. 

On the other hand a careful check of about fifty technical words used in mathematics 
revealed only one omission—namely ‘ Aufpunkt’. A general browsing revealed occasionally 
odd English (?) words; e.g. catacaustic. 

The typography is extremely clear, and the general lay-out of the book is commendable. 
Against each word is an indication of what part of speech it is, and the genders of nouns 
are given—but, unfortunately, not their plural forms. 

There is no doubt that most libraries and many individuals will find it necessary to 
possess both this and the subsequent volume. H. H. HOPKINS. 


Atomic Scientists’ News, Vol. 1, No.1. (London: Taylor & Francis, September 
1951.) Published bi-monthly, 6s. per issue. 


This is the first issue of a new series of the journal of the Atomic Scientists’ Association, 
intended for the general public as well as for the members of the Association. 

Even if the articles were expressions only of personal opinion their general tenor 
would inevitably be taken by the public as the ‘ policy’ of the Association, and it is clear 
from this issue, and in particular from the account of the annual general meeting, that it is 
aware of the dangers this would involve. It can offend nobody by having no policy at all 
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on the development of nuclear energy, or it can arouse enthusiasm but only among the 
few by presenting a forthright policy with probable political implications. The journal 
is at pains to avoid either extreme by adopting a policy primarily of scientific accuracy. 
It is an open question whether this will always keep the journal free from the charge of 
being political, but it is a policy to which no honest scientist can object. That it can be 
used with effect is shown by the devastating analysis of an article on Russian Atomic Power 
in World Review. 

The longest article is a review of recent progress in nuclear physics and engineering 
and related political events, and there is a most useful comprehensive account by Sir John 
Cockcroft of the long-term possibilities of nuclear power. The reader may care to guess 
what is concealed beneath Professor Peierls’ title Bathwater and the Baby. 

It is a well-balanced issue although the scientific terminology must prove rather 
difficult for the ‘ general’ public. If it continues in this manner the journal should be very 
useful in keeping the average scientist well informed about the physical facts and the social 
implications of the most dangerous branch of his activities. NaC.B. 


Advances in Physics : A Quarterly Supplement of the ‘ Philosophical Magazine’, 
Vol. 1, No. 1, Fanuary 1952. Editor: Professor N. F. Mort, F.R.S. 
Pp. 109. (London: Taylor and Francis.) 15s. 


Advances in Physics is a supplement to the Philosophical Magazine which will appear 
quarterly, from January of this year. The Editor, in his introductory remarks about the 


new journal, describes its purpose as follows : “‘ The aim of the journal is to publish specialist 
articles, of admittedly ephemeral value, in rather narrow fields, and to publish them 
quickly ”’. 


It is natural to compare the journal with the Physical Society’s annual Reports on Progress 
in Physics. Some differences are already clear. The intended practice of publishing 
small quarterly numbers, each containing about three articles linked to some extent in their 
subject matter, rather than large annual volumes covering diverse branches of physics, will 
open the way towards giving each number a specialist appeal rather than a wide interest. 
The effect of this, together with the policy of inviting ‘stop press’ accounts of subjects at 
present in a state of flux, should be that the new journal will complement the Reports. One 
imagines that people will turn to the Reports to seé how the frontiers of physics are advancing 
as a whole, and to the Advances to join in forays and thrusts that are now being made in 
particular fields. 

The three articles in this opening number deal with problems of crystals. Dr. 
Sondheimer writes on the mean free paths of electrons in metals, Professor Seitz on the 
generation of vacancies by moving dislocations, and Dr. Frank on crystal growth and 
dislocations. 

Dr. Sondheimer gives a most interesting account of recent work on methods of finding 
the mean free path of electrons in metals. The idea in all cases is to measure electrical 
resistance under such circumstances that this path is about equal to some other characteristic 
length in a metal. The simplest case occurs with thin films at low temperatures, where the 
mean free path may equal the thickness of the metal. Another case is a magneto-resistance 
effect, where the important length is the radius of a free electron orbit in a magnetic field. 
Lastly, there is the anomalous skin effect, which involves the depth to which a high frequency 
electric field penetrates a metal. 

Professor Seitz accepts the Editor’s hint by concentrating intentionally on ideas that are 
still mainly speculative. His main point is that vacant lattice sites can be generated from 
moving dislocations. The consideration of special processes, such as intersecting screw 
dislocations cutting through one another, shows that this idea is very reasonable, and several 
experiments exist which go some way towards confirming it. In his final pages Professor 
Seitz speculates in a most interesting and provocative way about the origins of 
work hardening. 

Dr. Frank in his article describes a delightful chapter of post-war physics, his own 
dislocation theory of crystal growth and its dramatic confirmation by the work of Griffin, 
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Dawson and others. It seems hardly necessary even to read Dr. Frank’s closely reasoned 


argument, for the beautiful photographs speak so eloquently in support of the main idea of 


his theory, which is that growth is catalysed at the points where dislocation lines emerge on 
the surface of the crystal. 
The Editor is to be congratulated on having made such a good start to his new journal; 
it will whet the reader’s appetite for future numbers. 
A. H. COTTRELL, 


Fourier Transforms, by Ian N. SNEDDON. Pp. xii+542. (London: McGraw- 
Hill, 1951.) 85s. 


The subject matter of this book is the practical solution of the boundary value problems 
of physics by means of integral transforms. Despite its title, it deals not only with those 
transforms in which the kernel is a trigonometrical function but also gives a detailed 
treatment of the Hankel transform and brief accounts of the Mellin and Laplace transforms. 
The first ninety pages are used to develop the necessary formulae for the transforms and 
their inversions; the rigour of the treatment is sufficient to meet most of the requirements 


of mathematical physics. The remaining four-fifths of the book show how the transform 
method has been applied to the solution of problems in many different parts of physics. 
Successive chapters are devoted to vibrations, heat conduction, neutron diffusion, 
hydrodynamics, atomic and nuclear physics, and stress distributions. Each of these 
chapters consists of a collection of examples, mostly taken from recent research applications 
of the transform technique. The examples have been chosen for their practical importance, 
rather than as illustrations of the method, with the result that the calculations are often’ 
extremely complicated. Many of them, particularly in the field of elasticity, are the original 
work of the author. 

Several of the examples show the value of using the less well-known methods of finite 
transforms and dual integral transforms. Relatively little space is given to the Laplace 
transform, but the author shows that it is sometimes advantageous to use other transforms 
in place of the Laplace one. It would have been useful if the author had indicated how 
it is possible to decide on the most suitable transform for any particular problem. One 
of the appendices consists of a short list of some of the commoner transform pairs. A more 
complete list would have been welcomed since, although other lists have been published, 
they are not always easily accessible. 

The greatest merit of this book is the vitality which it gains from dealing with problems 
of real practical interest. Inevitably only a small proportion of these problems will be of 
interest to any one physicist, but they may well suggest other applications of the transform 
methods. For this purpose the book might have been improved if there had been rather 
fewer examples in the text but more discussion of the conditions under which the transform 
methods are likely to succeed. Jo CosBs 


Photoelectric Tubes, by A. Sommer. 2nd Edition, revised. Pp. viii+118. 
(London: Methuen, 1951.) 8s. 6d. 


This volume represents the second edition of Photoelectric Cells previously published in 
this series by the same author.. The new title gives amore accurate description of the subjects 
under discussion which include all types of electronic tubes based on photoelectric emission 
and exclude photovoltaic and photoconductive cells which depend on barrier layer 
mechanisms. A new chapter has been added (contributed by Mr. S. Rodda) which deals 
with some of the fundamentals of photoelectric emission processes. This addition should 
certainly prove helpful to the type of reader for whom the book is intended, though there are 
some parts which may prove a little difficult, e.g. those dealing with the behaviour of electrons 
in surface layers and with photo-ionization. Other topics discussed in the volume include a 
description of various types of photo-cathodes, problems concerned with spectral sensitivity 
and the matching of phototubes to particular light sources, the operating mechanisms of gas- 
filled photocells and the properties of multiplier photocells. There is also a brief survey of 
applications in photometry, sound film reproduction and television, and a useful biblio- 
graphy gives the sources from which more detailed information can be obtained. The 
book is generally well written and will commend itself particularly to research students and to 
others who need an introduction to this important subject. H. K. H. 
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Variations et origine du rayonnement cosmique, by A. DavviLurrr. Pp. 83. 
(Paris : Editions de la Revue d’Optique Théorique et Instrumentale, 1951.) 
No price given. 


The study of the variations of the cosmic-ray intensity with time has already yielded 
much valuable information about the origin of the radiation, and it is appropriate, therefore, 
that these two aspects of cosmic-ray physics should be treated together in this book. In 
the first chapter there are sections dealing with the periodic and non-periodic intensity 
changes. ‘These include the solar and sidereal daily variations, the lunar variation, the 
seasonal variation and an account of the intensity changes associated with magnetic storms 
and solar activity. The second chapter reviews some of the theories and speculations 
about the origin of the radiation which have been advanced during the past thirty years 
or so. These include the annihilation hypothesis of Millikan, the cosmological theories of 
Lemaitre and Milne, and the more recent theories which ascribe the origin of the radiation 
either to stellar magnetic fields or to weak magnetic fields in interstellar space. These 
topics have received very little attention from previous authors of books on cosmic-ray 
physics, and this book will be of value to all who are interested in the geophysical and 
astrophysical aspects of the subject. H. ELLIOT. 


Combustion, Flames and Explosions of Gases, by B. Lewis and G. von ELBE. 
Pp. xix+ 795. (New York: Academic Press, 1951.) $13.50. 


In the past, progress in knowledge of combustion phenomena and the development of 
chemistry have gone closely together, from the days of the phlogiston theory to the important 
developments in chemical kinetics and the chain propagation theory around 1930. The 
more recent advances in combustion seem, however, to be more closely linked with various 
branches of physics. ‘This trend is reflected in the present book, in which a bare quarter is 
devoted to the purely chemical processes, and the remainder largely to such subjects as wave 
propagation, turbulence, heat transfer, thermodynamics, diffusion and minimum ignition 
energy. 

While the book retains the same title as the first edition in 1938, it is in fact almost entirely 
new, nearly all of it being rewritten, incidentally from a much more mature viewpoint, and it 
is over twice the original size. ‘The object is still to provide the chemist, physicist and 
engineer with the scientific basis for understanding combustion phenomena. It will 
undoubtedly become and remain for a long time the standard reference book on the subject. 

The book is divided into four main parts, headed ‘ Chemistry and kinetics of the re- 
actions between fuel gases and oxygen ’, ‘ Flame propagation ’, ‘ State of the burned gas ’ and 
‘ Problems in technical combustion processes’. There are also some valuable appendices 
giving numerical and graphical data for thermochemical calculations, limits of inflammability 
and flame temperature. 

The first part discusses both the general theory of chemical kinetics and the many 
individual reactions occurring in the combustion of hydrogen, carbon monoxide and hydro- 
carbons. This is an actively developing field, but there now exists a considerable area of 
agreement among workers and the more controversial aspects are not unduly stressed. 
Indeed the new chemical mechanisms put forward here may help to stabilize knowledge of 
the subject. The second part, on flame propagation, occupies about half the book; this 
deals with the propagation of the combustion wave in non-turbulent and turbulent mixtures 
and detonations, and is especially valuable for its discussion of the many aspects of flame 
stability and for the integration of the concepts of minimum ignition energy with results for 
quenching distances, spark ignition and burning velocity. ‘The whole treatment here is quanti- 
tative and precise. While the emphasis is perhaps on the basic scientific theory, the results 
are illustrated all through by excellent flame photographs of various types, many numerical 
data and many graphs and figures. On the experimental side the principles of flame photo- 
graphy by direct, shadow, schlieren and interferometric methods are described, and also the 
determination of pressure, and of flow pattern by particle track photography. The sections 
on measurement of burning velocity and of flame temperature are perhaps a little brief and 
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rather lost. There are some excellent photographs and an interesting account of the cellular 
structure of some flames, and of the so-called polyhedral flames ; the importance of the 
selective rates of diffusion of the gas constituents in causing these effects is stressed. 

On the whole this is a big task very well accomplished. The authors, in addition to their 
own great experience, have been able to call on the co-operation of several of the staff at the 
U.S. Bureau of Mines, and the whole book is a thoroughly authoritative treatment of the 
subject. 'There is obviously some difficulty in arranging a complex subject like combustion ; 
some of the sections are very long, that on cool flames and two-stage ignition running for 33 
pages and that on ignition by electric sparks for 35 pages without any sub-headings at all. 
Particularly valuable features are the numerous schlieren photographs of flame structure, and 
diagrams of flow patterns. The photographs and illustrations are a great asset in a book 
of this type. A. G. G. 


Progress in Cosmic Ray Physics, Ed. J. G. WiLson. Pp. xvi+557. (Amsterdam : 
North Holland Publishing Co., 1952.) 86s. 


The output of papers on the cosmic radiation is now so large that any attempt to assist the 
assimilation of the available data is certain to be welcomed. ‘The present volume consists. 
of eight separate articles, which are written by different authors on branches of the 
subject in which they are authorities. All the articles are primarily intended for those 
engaged on cosmic-ray research, but some of them may also be useful to those working 
on other subjects. The stress is always on the physical results rather than on the experimental 
methods. 

The four most successful reviews are those dealing with high energy nuclear dis- 
integrations, the nature of the primary cosmic radiation, observations underground and 
time variations. In each case the subject matter is sufficiently limited to enable that a fairly 
complete survey can be given. They are all most valuable in showing clearly both the 
extent of the knowledge which has been acquired and the problems which still require 
further investigation. The only theoretical review gives an account of the various types 
of fundamental particle and the possible couplings between them. It is a remarkable 
achievement that an article on theoretical physics should contain hardly any mathematics; 
even so, and despite the obvious care of the author, it is still by no means easy for an 
experimentalist to follow all the arguments. An article on the decay processes of heavy 
mesons suffers from the paucity of results available—it reads more like an original paper 
than a review article. A chapter which deals with recent data on geomagnetic effects 
contains a useful account of some experiments but does not review all the results which 
give information on these problems. An article on the equilibrium of the cosmic-ray beam 
in the atmosphere is the only one which overlaps with any of the others. The scope of 
this article is very wide; possibly in consequence of this, the data have not always been 
considered sufficiently critically. 

It will be seen that although this volume is not intended to be a complete survey, it does 
give an account of most of the more rapidly developing parts of the subject. The general 
production of the book is very good. Each section has its own list of references, but there 
is also a general author index—almost a bibliography of cosmic rays—and a subject index. 
The reproduction of the photographs is so good that the varying quality of the originals 
is clearly revealed. There are very few errors or inconsistencies between the various 
authors. ‘T’he editor and publishers are especially to be congratulated on producing such 
a complicated work in less than six months. The high cost is the only consideration which 
may deter some cosmic-ray physicists from acquiring this excellent volume. ifs pagk 
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Control and Removal of Radioactive Contamination in Laboratories, 
U.S. Department of Commerce, National Bureau of Standards. Pp. iv+24. 
(Washington : Government Printing Office, 1951.) 15c. 


X-ray Calibration of Radiation Survey Meters, Pocket Chambers and Dosimeters, by 
FRANK H. Day, National Bureau of Standards Circular 507. Pp. 11. 
(Washington D. C.: U.S. Department of Commerce, 1951.) 15c. 


Energy and Angle Distribution of the Photoprotons from Deuterium, by MartTIN 
WIENER, National Bureau of Standards Circular 515. Pp. 1413. 
(Washington D.C.: U.S. Department of Commerce, 1951.) 15e. 


Nuclear Data, Supplement 2 (July 1950 to January 1951) to National Bureau of 
Standards Circular 499. Pp. 63. (Washington, D.C.: U.S. Department 
of Commerce, 1952.) Price, $4.25, includes Circular 499 and three 
supplements, the third of which is scheduled for publication in 1952. 


CONTENTS FOR SECTION B 


PAGE 
Dr. F. C. CHampion. Electrical Counting Properties of Diamonds. : 2) 465 
Mrs. K. Stratton and Dr. F. C. Cuampion. Electrical Counting eg of 
Two Large Diamonds under Beta-Irradiation : : : 473 
Prof. A. L. REIMANN and Mr. J. V. SULLIVAN. Rectification Phenomena Exhibited 
by Natural and Sulphurized Galena . : 480 
Mr. A. N. Gent and Mr. R. S. RIVLIN. Sonate, on a Mechanics of Adee: 
II—The Torsion, Inflation and Extension of a Tube . : 487 


Mr. P. GrooTennuis, Dr. R. W. PowELt and Mr. R. P. TYE. Thea “ae 
Electrical Conductivity of Porous Metals made by Powder Metallurgy Methods 502 


Mr. R. S. Barnes. Effects associated with the Flow of Vacancies in Intermetallic 
Diffusion P d ; ; : 4 : : : i : x SW 


Mr. A. R. Verma. New Observations of Crystal Overgrowth on Silicon Carbide . 525 


Mr. F. R. L. ScHoenrtnc, Dr. J. N. vAN NIEKERK and Dr. R. A. W. Haut. 
Influence of the Apparatus Function on epee Size Determinations with 


Geiger Counter Spectrometers ; 528 
The late Mr. H. Lunp and Mr. L. Warp. The Spectel Eee nf ae 

Nickel and Cobalt . : 535 
Dr..K. E. Spetts. A Study of Gipeutaton Petey within Liquid Drops a 

through a Liquid . ; Beal! 


Mr. E. W. B. Git and Mr. G. F. ALFREy. The Roe heron of Liquid Dene 7 a46 
Letters to the Editor : 
Prof. P. SeLeNy1. The Influence of Hess Vapour on Selenium Rectifiers 


and Selenium Photoelements . : ‘ : : x 4 DE 

Mr. H. A. Geppre and Mr. D. G. Kizty. The Dielectric Constant and Loss 
of Amorphous Selenium at a Wavelength of 3cm_. : : ; 5 See 
Dr. R. J. Weiss. Extinction Effects in Powders : : ; : 5 oS 
Corrigenda . : ; , : : ‘ : : ; 3 ‘ SEY 
Reviews of Books ; : : : : ; ; ; 5 3 . 2555 
Contents for Section A : : : ; ‘ ‘ : : f 5 gelete: 


Abstracts for Section A ; : : é : P : ; ‘ 5. SS 


574 


ABSTRACTS FOR SECTION B 


Electrical Counting Properties of Diamonds, by F. C. CHAMPION. 


ABSTRACT. Examination was made of over 200 diamonds, most of which were of gem 
quality and many of which weighed only a few milligrams. About 25 of them showed good 
conduction properties when submitted to ionizing radiations. For depths of crystal 
penetration which exceeded a certain amount, all the counting specimens responded to 
both «-particles and f-particles with a response which was at least qualitatively dependent 
on the energies of the particles. In some cases, however, a marked drop in the counting 
properties was easily demonstrated with a-particles under conditions of feeble penetration. 
Even at maximum penetration 5-mev a-particles produced oscillograph pulses in such 
crystals which were only about three times as large as those from 1-mMev f-particles. 
Strong correlation was established between the wavelength of the continuous ultra-violet 
transmission limit and the magnitude of the electrical counting response of the specimen. 


Electrical Counting Response of Two Large Diamonds under Beta-Irradiation, by 
K. StraTTON and F. C. CHaMPION. 


ABSTRACT. 'Two large diamonds §, and S,, both of good optical quality, were subjected 
to B-irradiation and the electrical conduction responses were studied under a variety of 
conditions. With S, special attention was paid to the formation and removal of space 
charge. There is evidence that the mechanisms of space charge removal by /-irradiation 
‘and by illumination respectively are not quite the same. The efficiency of this diamond 
as a B-particle counter was not appreciably different from that of a standard Geiger counter. 

Small volume elements of S, and S,; were examined individually for S-counting 
efficiency. While S, did not show a variation greater than 20°{, variations as large as 
ten to one were found in different parts of S;. ‘The results were not inconsistent with the 
existence of a correlation between the ultra-violet transmission and the counting efficiency 
of any selected volume element. 

The f-response of S,; was more than three times the best response from S,; under 
comparable conditions, while the ultra-violet transmission limits were 2250 and 
2850 A respectively. : 


Rectification Phenomena Exhibited by Natural and Sulphurized Galena, by 
A. L. REmmann and J. V. SULLIVAN. 


ABSTRACT. Specimens of n-type galena were sulphurized, either by a wet chemical 
process or by exposure to sulphur vapour at various temperatures, and the resulting 
changes in rectification characteristics were observed, both with d.c. instruments and 
at 3000 Mc/s, using various cats-whiskers. According to treatment given, either an 
improved n-type rectification was obtained or the polarity was changed to p-type. The 
results are briefly discussed in the light of present-day rectification theory. 


Experiments on the Mechanics of Rubber : II—The Torsion, Inflation and Extension 
of a Tube, by A. N. GENT and R. S. RIVLIN, 


ABSTRACT. 'The formulae derived by Rivlin (1949) for the forces necessary to produce 
simple extension, radial inflation and torsion in a tube of incompressible highly elastic 
material, which is isotropic in its undeformed state, are specialized to the case when the 
amount of torsion is small. The resulting formulae are compared with measurements on 
tubes of vulcanized natural rubber and conclusions as to the form of the stored-energy 
function are drawn, in general agreement with those reached from earlier experiments. 
The experiments indicate that hysteresis between the loading and unloading curves for 
the vulcanizate is associated with certain terms in the expression for its stored-energy 
function and, presumably, with the mechanism which accounts for the presence of these terms. 


Abstracts for Section B 575 


Thermal and Electrical Conductivity of Porous Metals made by Powder Metallurgy 
Methods, by P. Groorenuuts, R. W. Powe t and R. P. Tye. 


ABSTRACT. Measurements are reported of the thermal conductivity and electrical 
resistivity over the range of 20° c to 200° c on several specimens of a sintered porous 
bronze material (89%, Cu, 11% Sn). The constituent powder particles were substantially 
spherical in shape and the specimens studied covered particle diameters ranging from 
0-00133 to 0-040 cm and had densities from 5:27 to 7:01 g/cm*. The modified Lorenz 
relationship between the thermal conductivity K the electrical resistivity p and the 
absolute temperature T was found to apply, the experimental results satisfying the straight 
line K=0-58 x 10-8T/p+0-005. Existing data for solid bronzes of similar composition 
also conformed to this line. It is shown that, to a first approximation, the conductivities 
are independent of the particle size of the powder from which the specimen has been made, 
but can be correlated to. a density basis, the values falling along a straight line drawn from 
the point for the solid material to cut the abscissa at a density corresponding to the 
maximum porosity which can be attained on packing equal sized spheres. This method 
of correlation is also confirmed by data from the literature for the electrical conductivity 
of copper, nickel and iron. Suggestions are made for estimating the thermal conductivities 
of similar porous materials. 


Effects Associated with the Flow of Vacancies in Intermetallic Diffusion, by 
R. S. Barnes. 


ABSTRACT. Experiments on the interdiffusion of copper—alpha-brass and copper—nickel 
couples and sandwiches, have shown that there is an increase of volume in the diffusion 
zone as a result of diffusion. This increase in volume is associated with the formation, 
in the diffusion zone, of voids with crystallographic faces. "These voids appear near to 
the original interface and on that side from which there is a net loss of atoms, or a net gain of 
vacancies, as shown by the movement of inert markers. Grain boundaries influence the 
movement of the markers in regions where there has been preferential grain-boundary 
diffusion. A close micrographical examination of the diffusion zone shows it to have an 
unusual structure, various observations pointing to the presence of strain, and x-ray back 
reflection photographs of a copper—nickel couple show that the diffusion zone is polygonized. 
Subsidiary experiments suggest how strain in the diffusion zone arises. 

A vacancy diffusion mechanism where there is a preferential flow of vacancies in the 
one direction, is used to explain the above phenomena. This flow will produce strain in 
the diffusion zone if vacancies are generated on one side of the original interface and 
collapse on the other. Mechanisms are discussed for the generation and condensation, 
within the diffusion zone, of the large numbers of vacancies necessary for the process. 
To explain the volume increase and formation of voids it is postulated that not all the 
vacancies on the one side of the interface are eliminated by the lattice collapsing, but that 
some vacancies coalesce to form microscopic voids which enlarge by assimilating further 
vacancies. 


New Observations of Crystal Overgrowth on Silicon Carbide, by A. R. VERMa. 


ABSTRACT. On the faces of some silicon carbide crystals showing ‘ growth spirals ’ 
superimposed oriented overgrowths can sometimes be observed. ‘These oriented over- 
growths which took place during the process of manufacture consist of stepped triangular 
pyramids with craters at their peaks, the pyramids having heights up to a few wavelengths 
of light, and built sometimes on platforms of similar heights. ‘These conclusions have been 
arrived at by the application of the techniques of multiple-beam interferometry and 
light-profile microscopy. In addition to the pyramids some trigonal, hexagonal and their 
compound figures have also been observed, which by the application of multiple-beam 
interference fringes are shown to be molecular oriented layers. These oriented molecular 
layers, the existence of which is revealed from the observation of the above overgrowths 
on silicon carbide crystals, may be silica deposit or possibly the growth of silicon carbide 


on itself. 
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Influence of the Apparatus Function on Crystalhte Size Deen with Geiger 
Counter Spectrometers, by F. R. L. ScHoENING, J. N. vaN NIEKERK and 


R. A. W. Hawt. 


ABSTRACT. A brief survey of the various methods used for determining crystallite 
sizes from x-ray diffraction line broadening is given. It is shown that the diffraction 
profiles determined experimentally with a Geiger counter x-ray spectrometer can be 
represented by functions of the type y=c/(1+-’x*)?. Using these functions a correction 
curve for instrumental influence on line breadth measurements is derived by solving the 


++ © 
integral equation B(x) = b(v) B(x —v)dv by means of the Fourier transform method. The 


apparatus function b(x) is determined by using a standard substance. Experiments are 
described which illustrate the effect of sample preparation, absorption, etc. on measurements 
-of line breadth. 


The Spectral Emissivities of Iron, Nickel and Cobalt, by the late H. LuNpD and 
L. Warp. 


ABSTRACT. ‘The spectral emissivities of iron, nickel and cobalt have been measured by 
the direct comparison of surface and black body radiation at temperatures between 1000° c 
and 1 300° c, over the wavelength range 1-0: to 2:6, using a lead sulphide cell and amplifier 
as the detecting unit. Values of the temperature coefficients of emissivity referred to 
1000° c indicate an X point for iron at 1:6, but none for-nickel in this spectral range. 
The results for cobalt show an anomaly at a temperature close to the Curie point. 


A Study of Circulation Patterns within Liquid Drops moving through a Liquid, 
by K. E. SPELLs. 


ABSTRACT. The technique used to reveal the circulation patterns was developed by 
Hagerty, and depends on the striae which appear in glycerine and glycerine—water mixtures 
when sheared. 

Instantaneous flash photographs were taken of circulation patterns in drops of glycerine 
falling slowly (about 1 cm/sec) in castor oil, and in drops of various glycerine-water 
mixtures falling at about 10 cm/sec in heavy white oil (viscosity 0°37 poise at 15°c). In 
the former case the drops were spherical and the circulation pattern agreed reasonably 
well with that predicted theoretically by Hadamard; in the latter case the drops were no 
longer spherical and the centre of circulation was displaced below the equatorial plane. 


The Electrification of Liquid Drops, by E. W. B. G1Lu and G. F. ALFREY. 


ABSTRACT. Charges on drops breaking up in electric fields were investigated and 
it was found they could all be explained by the ordinary laws of electrostatic induction, 
with allowance for contact potential differences. Electrification due to splashing was 
also studied and again the effects could be explained by contact potentials alone so that 
theories put forward of double layers to explain previous experiments were unnecessary. 
A warning is given of the necessity of making due allowance for possible electrification 
by splashing in interpreting experimental results and a critical examination of an 
experiment of Zeleny is reported. 

A few experiments of splashing water drops off ice gave results which may have a 
bearing on the production of electricity in thunderstorms. 


